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ABOUT THE BOOK 


This book is primarily inten- 
ded as a text on queueing theory 
for students of graduate pro- 
grammes in management 
business administration, mathe- 
matics, statistics, computer 
science and the like. It may also 
be found useful by advanced 
undergraduate students and O.R. 
practitioners. The theory of 
queues comprising description of 
various queueing models has 
been given in the first four chap- 
ters, pointing out applications to 
practical problems form time to 
time. The last, ie. the fifth, 
chapter is exclusively devoted to 
applications of queueing theory 
to telephone traffic, inventories, 
machine interference, scheduling 


and transportation. 


Some special features of the 
book are : 

1, A problem set is provided 
at the end of each chapter, with 
answers given for all problems 
and hints for solution or complete 
solutions given in several cases. 

2. Worked-out examples are 
provided in the text, wherever 
found necessary, 

3. For most measures like 
queue length, waiting time, etec., 
variances are evaluated in addi- 
tion to means, often with proce- 
dures appropriate for evaluation 
of higher moments. 

4. Almost throughout the 
text, explicit results have been 
obtained after inversion of 
Laplace transforms or generat- 
ing functions. 

5. Preliminaries and other rele- 
vant data are gi 
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PREFACE 


The present book is primarily intended to serve as a text for graduate 
students studying operations research in institutes/schools of management] 
business administration or of technology/engineering or in allied departments 
like mathematics, statistics and computer science of universities/colleges. It 
may be found useful by advanced undergraduate students as well as by O.R. 
practitioners, who need a greater exposure to queueing models than that 
available in texts on operations research containing a chapter or two on the 
subject. 


The book has grown out of the lecture notes used by the authors for about 
two decades, modified from time to time and enriched by some of their recent 
research findings as well as those of contemporary authors. Our aim has 
been to present the subject matter in a simple, lucid style, devoid of technical 
jargon as far as practicable. Our thoughts on the subject and our notes 
have drawn considerably from the existing literature on queueing theory, 
especially from the books written by authors like Cox and Smith, Morse, 
Prabhu, Saaty, and Sasieni, Yaspan and Friedman, and some more recent 
books including those by Chaudhry and Templeton, Gross and Harris, 
Kleinrock, and Neuts, and to them, among others, we acknowledge our sincere 
gratitude. Of course, such a list can never have the semblance of being 


exhaustive. 


The book consists of five chapters. Chapter 1 contains a brief introduc- 
tion to the subject and goes on to the study of steady-state Poisson queues. 
Starting with the M/M/1! queueing system, it includes various other Poisson 
queueing models like M/M/1/N, M/M/c, M/M/oo, studied through the use of 
state-dependent parameters. Chapter 2 includes the transient solution of 
the M/M/1 queue and the analogous derivation of the busy period distribution. 
In Chapter 3 we study steady-state non-Poisson queues like M/G/1, M/D/e, 
GI/M/1, and GI/M/e, using the imbedded Markov chain method and append 
a continuous-time analysis of the M/G/1 system, using the supplementary- 
variable technique. Also included are the continuous-time study of the 
M/E, /1 system as as well as results for M/HE/1 as a special case of the M/G/1 
queue. The chapter concludes with a matrix-geometric analysis of the 


GI/E_/1 system. In Chapter 4 we include : 
al /M/1, the bulk-service model M/M?/1, queues i 


miscellaneous topics like the bulk- 


arrival queueing system Me 


x 


series and priority queues. For simplicity, only the case of non-preemptive 


priorities is analyzed. Chapter 5 includes a discussion of practical applica- 
tions of queueing theory to areas like telephone traffic, inventories, machine 
interference, scheduling and transportation to serve as an illustrative set of 
areas of queueing applications, though, really speaking, queueing and con- 
gestion problems arise almost in every walk of life. 
Finally, as complement to the main text, six appendices have been added. 
As, for various reasons, many queueing problems are often solved through 
simulation, we thought it proper to provide a concise introduction to simula- 
tion and its application to queues and this has been done in Appendix A. 
Necessary information on probability as well as on the Laplace and Z- 
transforms is given in Appendix B, whereas Appendix C contains miscellaneous 
topics like Rouché’s theorem and some not-so-common topics from differential 
and difference equations, and inversion of transforms. Appendices D and 
E provide a list of the symbols used and a table of Laplace transforms, 


respectively. Appendix F is a table of random digits for handy use in 
simulation. 


Some words regarding the special features and the raison d’étre of this 
book will not be out of place here, especially when a rich literature on the 


subject is already available. Some of the special characteristics of the book 
are the following. 


1. All chapters contain Problem Sets at the end. Answers to all pro- 


blems are given and, in some cases, the problems are completely worked out. 
2. Wherever considered necessary, suitable solved examples have been 
included in the text to give the reader an insight into problem solution. 


3. Often the mean is found to be a crude indicator of an efficiency 
measure of a system. For most of the measures like the queue length, waiting 
time, etc., we have evaluated the variance as well, sometimes with a procedure 
given for evaluating higher moments. 


4. For almost all topics included in the text, complete solutions are 
given, including the inversion of the Laplace transforms or generating func- 
tions, so that the results are available for use in practical applications. 

6. To sum up, though the book is expected to make an interesting 
reading for any one who has gone through undergraduate mathematics and/or 
statistics courses, we believe that, with the simple language used and with 
inclusion of topics in Appendices B and C, it can also be used fruitfully by 


the not-so-mathematically-oriented students of business, and practising 
managers, 


Sr alall 


ig 


TS 


eet 


— <r + -Seie 
on wot 


xi 
To quote Francis Bacon, ‘Studies serve for delight, for ornament, and 
for ability.” We believe this book will be found useful on the three sotiitta 
We are very grateful to Professor J. G. G. Templeton of the 
of Toronto for having kindly agreed to go through the entire text a 
suggestions for improvement. 


University 
nd offering 
Members of our families have helped us in many ways during the writing 
of this book, and to them—Raj, Rakesh, Meenakshi and Rajesh Kashyap, 
and Parphulat, Vandana, Amit and Ajit Chaudhry—we owe many thanks. 
Last, but not least, we thank the printers, Eka Press of the Indian 
Statistical Institute, Calcutta for their excellent work. 


We will welcome the readers’ comments and suggestions for improvement 
of the book. 


B. R. K. Kashyap 
M. L, Chaudhry 
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CHAPTER 1 
POISSON QUEUES | 


1.1 Introduction 


eae 


In queueing theory we study situ 


at a service facility for receiving sey 
having to wait for service, 


ations where units of some kind arrive 


vice of some description, some of the units 
and depart after service. A queue or waiting line 
service facility cannot cope with the number of units 
The units calling for are often called customers in 


fined as something hz 


develops whenever the 
requiring service, 


service 
a generic sense. 


A system is generally de wing an input, 


in between, which changes the input 
Thus, the arrangement 


regarded as a queueing system, 
the input, the serviced ¢ 


an output and a transformation process 


into the output. described above qualifies to be 


where the customers requiring service 
ustomers the output and the service 


form 


rendered is the 
transformation process. 


~ a51AIIS o1qNd uolun fe | spaepue} 


customers 

0::-00 O i 
—— 

input } output 


service 
facility 


The subject of queueing theory had its origin in the pioneering work 
done by Agner Krarup Erlang, an engineer at the Copenhagen Telephone 
Exchange around the beginning of this century on the application of probability 
theory to telephone traffic problems. It soon drew attention of many other 
probability theorists and has remained a popular field of research —— 
throughout the eighty years since then. There are many situations in - 
life where queues of the type described above develop. Eh we can have a 
queue of broken-down machines waiting for repair at a repair shop, a queue 


. ee 


ene) La ie 
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of customers at a departmental store cash counter, 


we have cases where a physical queue is absent, such as the Ww 
passengers for a railway or airline ticke 
say, for the purchase of a car. 
supplied from future production. 
1.2 Characterization 

A queueing system is specified by stating the following details about it, 
(a) The input or the arrival pattern, 


System occur at times ty, ty, ...; then u, = bri;—ty (r = 1, 2 
arrival times. We assume that the » 

distributed (i.i.d.) random variables (r.y.’s) 
by stating their probability distr 
called the inter-arrival time distribution or, 
the input distribution. Such information as whethe 
in groups may be added here, 


(b) The queue discipline. 
customers for service out of those Waiting, and may include rules regarding 
formation of queues. 

(c) The service mechanism. It 


is stated by Specifying the number of 
servers ¢ > 1 and the 


probability distribution function B(v) 
customers »,, Usper ns 
.8 and also independent of the inter 
may also state whether the se 


of fixed or varying size.* 


of the service 
times of the successive 


which are assumed to be iid. 
r.v 


-arrival times uw, (r = 1, 2 +). We 


’ Ma J 
rvice is rendered to customers singly or in groups 


1.3 Topics of Study 


In any analysis of 


a queueing system, one or more of the 
of it are investigated, 


following aspects 

(a) The queue length. This means the number 
before service, Sometimes the term i 
customers present in the sy 


of customers waiting 
8 used to mean the total number of 


stem, including the one or more 
are being served at the moment. Generally, the precise 


the context. However, some author 
term queue length in the for. 


customers who 
meaning is clear from 
8 prefer to avoid confusion by using the 
mer sense, i.e, to mean the queue proper (before 


* Some recent studies on queues with correlated arrivalg or service times have been made, 


[ Sec. ifs?) 
or such dangerous queues 
as the one formed by planes circling above an airport waiting to land. Often 
aiting list of 
t or of persons who register their names, 


which is not readily available and is to be 


Let the successive arrivals to the 

, 2, ...) are the inter- 

r are independent and identically 
and the input process is specified 
ibution function (P.D.F.), A(u). A(u) is 
simply, the arrival distribution or 


r arrivals occur singly or 


It is the statement of the rule of selection of 


— 


f 


y 
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“2 
service), and by calling it the 


number in the System or the 
those being served 


system size when 
are intended to be included, W 


e shall be following this 
arrangement, 


(b) The wailing time, 
waiting before service. 
authors 


This means the time 
This term, too. 
to mean the time spent 
call the time spent iy 


spent by a customer in 
is used in a wider sense by some 


in waiting or being served, and they will 
i waiting hefore 


service the queueing time. We shall 
prefer the former meaning, which seems to be more popular, and we shall 
call the waiting time in the w ider sense the total time spent in the system, 
which is clearly the waiting time (before service) plus service time. 


(c) The busy period, By a busy period is meant & period of time over 
which the server ig continuously busy. Clearly, a busy period starts with 
an empty service facility and ends when the 
in the system drops to zero for the first time. 


Starts, which ends at 


an arrival at number of customers 
At this instant an idle period 
the commencement of the next busy period. 

The queue length, the waiting time and the busy period are studied 
through their probability distributions, from which 


moments like mean, 
variance, ete. can be obtained. 


1.4 Notation 


We shall use Kendall’s (1953) notation, explained below, to denote the 
various queueing systems to be discussed. Occasionally, when necessary 


we shall enhance the notation following the extension of Kendall’s notation 
suggested by Lee (1966). 


A queueing system will be denoted by a triad -/-/-, in which the first 
two members will be letters which stand for the forms of the input and the 
service time distributions, respectively, and the third member is a angie 
signifying the number of servers. The forms of the arrival and service dis- 
tributions commonly used are described below. 

(a) M (Markovian or the neyative exponential diairebulion). This is the 
most commonly used distribution in queueing literature for reaeoca tobe 
discussed later. When the letter M is used in the first piace in the tried 
mentioned above, it will be taken to mean that the input distribution is 
given by 


(1) A(u) = 1—e"**, u>O 


= 0, otherwise 


we 1 By 
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or, that the corresponding inter-arrival timo density function is given by 


(2) a(u) = Ae*¥, u > O 
- 0, otherwise 


The mean of the above distribution is 1/A, which is, therefore, the mean inter- 
arrival time. It follows that A is the mean arrival rate, ¢.e. the average 
number of arrivals per unit time. From the well-known relation between 
the (negative) exponential and the Poisson distributions, it can be seen that 
the number of arrivals in a given time ¢ will have the following Poisson 


distribution : 
(3) P{n arrivals in time ¢} = (At)” e-4#/n!, n = 0, 1, 2, ... 
Because of the above, the case of exponential inter-arrival times is equi- 
valently stated as Poisson arrivals or Poisson input. 
In the same way, when M stands at the second place in the triad, it will 
mean that the service time distribution is 
(4) Biv) = 1—e-*”, v > 0 
= 0, otherwise 
the corresponding density function being given as before, with A, w in (2) 
replaced by mw, v. As before, we see that 1/u is the mean service time and 
fe the mean service rate, and 
(5) P{n services in time t} = (ut)" e~*#/n!, n = 0, 1, 2, ... 


Note: Special features of the negative exponential distribution. The 
exponential distribution possesses the following two properties, which make 


it very convenient for use in the analysis of queueing systems. 

(i) The Markovian or the “forgetfulness’’ property. The conditional 
probability of an arrival during the time interval (u, u-+-du), given that an 
arrival has not taken place in [0, uo], is equal to 
(6) Plu < Uy < u-+-du| up > ug} = a(u)du/{1—A(uy)] 

= Ae*utty) dy 


which shows that the entire inter-arrival time and the residual inter-arrival 
time after time uw) have the same exponential distribution. Taking Up =U, 
we see that the probability of an arrival during (u, u-|-du) does not depend wpon 


the time since the last arrival. That is why Poisson arrivals are also called 


random arrivals, 


(4 
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Similarly, in 


time and the 

naa 1K ; have the same distril 
ty of a service completion is inde 

to this fact, 


the case of ex i 
he case of ©xponential service times the 


be residual service 
ntire service time 


ution and the proba- 
pendent of the clapsed service time Due 
distribution is not suitable for 
the ©xponential distribution is often a 
especially where 
machino failure, or 


the exponential vi 
On the other hand, yee 


arrival times good fit for inter- 
; arivals are due to chance phenomena like 
come from an unstructured 


system like bus arrivals at ; 
bus stop far from the ; ween 


Starting point, 


In many instances ‘ / 
: any Btancoes oustomers’ arriving 
at a business shop exhibit | 


such random arrival pattern. There is a saving 


in Hindi hie Nf ; 

indi which means that a customer and death have no fixed time 
ii) In the case of nopat; i istri 
(il) In the eas of ne gative-exponentially distributed inter-arrival times 


since the arrivals in a given interval follow the Poisson distribution, we have 
(7) P{one arrival during (¢, t- At)} = (AAtje>44/1! 
AAHI1—AAL+...) 
A At--o(At) 
where o(At) stands for terms which are negligible 
sense that 


in comparison to Af, in the 


(8) lim aad) 0 
Sst—>0 At 
Similarly, 
(9) P{no arrival during (#, t-+-At)} = (AAL)% AAtIO| 
L—AAt--o(At), 
(10) P{two or more arrivals during (f, (--Af)} <= o(Ad) 


In the same way, if service times are alao exponentially distributed, 


wo have 

(11) P{one service completion during (1, t-+Af)} = wAt--o(At) 

(12) P{no service completion during (¢, t-|-At)} = 1— aAt-+-o(Af) 

(13) P{two or more service completions during (t, ¢--Af)} = ef At) 

(14) Pino arrival and no service during (f, f-}- Af)! 1 —(A-+- 2) At--o( At) 
(15) Pfone arrival and one service during (f, t-+-Atl)} = o( Af) 


(b) D (deterministic or constant inter-arrival/service times). The letter D 
in the first (second) place in the triad in Kendall’s notation indicates that 


the inter-arrival time (service time) is constant. In the case of constant 


uoIun Pty 


NNR ee i 
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a > Z i é 2, é ut, W hen each 
Ai LI calle Peowlar ary vals or regular inp ut, 
ive es é 1 30 all ed reg ular arTl 
In r-arT1v al times, £ 


7 2 i istribution 
inter-arrival time is a, say, we have the input distribu 


(16) A{u)=0, waa 

= 1 sa Sia 2(0: <= od) 
and the density function 
(17) a(u) = d(u—a), 
where d(-) is the Dirac delta function. When service times are constant, 


we have similar expressions for B(v) and b(v). . . 
(c) By, (k-Erlang distribution). We say eas the Sai a 
have a k-Erlang distribution, or an Erlang Cipro, of i es fs eu. 
inter-arrival time is made up of k phases, the time spent in each p be g fe 
a negative exponential distribution with a un i a pe 4 ae 
independently of other phases, so that the mean inter-arrival tu s 1, 


Thus, the inter-arrival time density a(u) is given by 
(18) a(u)du = Phu < uy < u+tdu} 


— P{(k—1) phases are completed in time uv and 
. one phase in (wu, w+du)} 


co) oom _ kAdu 


(kA)Fuk-1 e- kau 
- (E11 


This is clearly the k-fold convolution of the exponential distribution ee 
If and Also, it 


du. 


itself and, therefore, gives the exponential distribution for k = 1. ; 
. , i i re at the Laplace 
ean be seen from (18), or by using the convolution theorem, that the Lapla 


transform of a(w) is 


a kA \# Oa as Glee 
(19) J emawdu = (qy55) =(+—) 78 


as k-»o, Now, since the Laplace transform of the density function of a 


i i (e7 re see that e~*/* corresponds to the case of 
r.v. X is the expectation H(e-8*), we see that e correspo 


constant inter-arrival times = 1/A, the mean value. So, for ko, the 
Erlangian input becomes the regular input, given by (16) or (17). 

The case of k-Erlang service time distribution can be dealt with similarly, 
the density function b(v) being given by (18) with A, uw changed to yp, », 


respectively. 


() 
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(d) G (general distribution). The letter G 


ace 1s used for a general or 
arbitrary distribution. However, . 


rs . following the convention, we shall write 
an t for *foanan, . ri . : 
a1 (standing for general independent ) for general inter-arrival times and 
‘ ~ _ . . . . : ; 
G for general Service time distribution, 


From the notation explained above, we see that 


(i) M/M/1 stands for Poisson arriy 
one server, 


als, exponential Service times and 


(ii) M/D/e stands for Poisson arrivals, 
c servers, 


constant service times and 


(ui) M/E,;/2 stands for Poisson arrivals, k-Erlang service 
two servers, 


times and 


wT Tift], > 
(iv) GI/G/e stands for the most general case of general input, general 


service times and ¢ servers, 


Lee (1966) has extended the above notation to -/+/- :(- | ‘), where the 


triad is as explained earlier, and within the parentheses the first member is a 
number, showing the maximum number of units that can be accommodated 
in the system (queue plus service) and the second member is an abbreviation 
for the queue discipline, like FIFO. LIFO, SIRO, ete., explained below. 


We shall now describe some of the queue disciplines encountered, 


Queue discipline. The most commonly practised queue discipline is 
“first come, first served” or, as it is sometimes called. “first in, first out”, 
abbreviated, respectively, as FCFS or FIFO. In this case, the waiting 
customers are selected for service in the order of their arrivals, Various 
other possibilities are there. Sometimes we have just the opposite rule of 
“last come, first served” or “last in, first out’, in short LCEFS or LIFO. In 
Some cases we can have a random selection for service, or “service in random 
order”, briefly, SIRO. And there are also a number of priority disciplines, 
according to which the customers are divided into two or more priority 
classes and a higher priority unit is chosen for service before a lower priority 
unit. Further, under a “preemptive priority” rule, a lower priority unit is 
preempted or taken out of service whenever a higher priority unit arrives 
during its service, the preempted lower priority unit being attended to only 
after serving all the units of higher priority classes, Under “non-preemptive 
priority’ or the “head-of-the-line priority” rule, a service once started is 


allowed to continue up to completion, irrespective of arrivals of higher 
priority units, which are attended to only after completion of the service in 
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ef i jority” disciplines, where the priority 
progress. There are also “dynamic priority disciplines, wher I Vv 


i * »riori isciplines in 
increases due to waiting. For a detailed account of priority discipli 
queues, one may refer to Jaiswal (1968). 


It may be noted that out of the three queue characteristics of the queue 


a ‘ ati ‘me is affected 
length, the busy period and the waiting time, only the waiting time is affect 
by the queue discipline. 


Impatient customers. In addition to the queue ciaeiphe some bag 
of behaviour of waiting customers have often to be taken into account. Thus, 
a customer is said to “balk” if, seeing the size of the queue, he Sead nob 
to join the queue and is thus lost to the system. If after wailing ior corny 
time, he gets impatient and leaves the system before his service, he is said 
to “renege’’. 


Systems of queues. In a multi-server queueing system like M/M/c, the 
servers are operating in parallel and are fed by a single queue. Alternatively, 
we may have a separate queue for each server, @.c. a number of quente 
parallel. In such a case, the customers exhibit the phenomenon of shifting 
between queues time and again, for example, to be always in the ahorterl 
queue. This is called jockeying. We have an arrangement of gieue in series 

or queues in tandem when the output of one queue becomes the input to the 


i in the series. We have a case of 
next queue and so on until the last queue in the series. We have a ce 


cyclic queues if the output of the last queue again becomes input to the first 


queue. Sometimes we have queueing networks consisting of queues arranged 
both in series and in parallel. 


1.5 The Queueing System M/M/1 
(a) The queue length 
Let us consider the single-server queue with Poisson arrivals and expo- 
nential service times. Let us write 


(1) 


P(t) = Pf{there are n units in the system at time #} 


So, using eq. 1.4 (7)-(15), we have, for n > 1, 

(2) P,(t+-At) = P,,()P{no arrival and no service in (¢, (4-At)} 
+P,,_,(t)Pf{one arrival in (t, t-+-At)} 
-+-P.,,.,(t)P{one service in (¢, t-+At)}+-o(At) 


= P,(){1—(A+m)AB-+P,, a(t) AAL+P, (0) wAt-+0(At) 
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whence 

; P,(t+At)—P,,(d) A; 
(3) AE = ALM O+AR, (0+ AP gas(t)-+ 200 


vad 
AL 

Pacaeeal jes 

roceeding to limits as At > 0, we have 


(4) 


1 
ai P(t) = —A+p)P,()-+AP,()-+HP, A(t), n= 1,2, ... 


For n = 0, we shall have 


(5) Po(t+-At) = P,(t)P{no arrival in (¢, + At)} 


+ P,(t)P{one service in (E, t--At)}+-o(At) 


= Py(t){1—AAL} + P,(t):wAt--o( At) 
whence, as before, 


d 
(6) dpe!) = —AP ot) +HP,(t). 


Kq. (4) and (6) are differential-difference equations, which can be solved 
to evaluate P,(+), thus giving the time-dependent probability distribution 
of the number in the system. This solution is quite involved and will be 
discussed in the next chapter. 


However, for most of practical applications 
of the queueing model, a “steady state” solution suffices. 


We say, in general, 
that a system has reached a steady state or a state of “statistical equilibrium” 
when fluctuations with respect to time can be disregarded. Hence, to obtain 
a steady state solution we shall regard the P,(t), n = 0, 1, 2, ... as independent 
of time. Consequently we have dP,(t)/dt = 0, and we shall drop the argu- 
ment ¢ in P,(t) and write simply P,,. 

(7) 


(8) 


Thus, (4) and (6) yield 
0 = —(A+p)P,+AP,iteP aap n= 1,2,... 
0 = —AP)+yP. 


Kq. (7) and (8) are easy to solve recursively, though there are other methods 
to solve them. (For details see Section 1.7.) It is customary to write 


(9) (A/z) = p 
and. call this ratio the “traffic intensity” of the system. 


From (8), we have 
(10) 


P, = pPo 


Now, writing eq. (7) for n = 1, substituting for P,, and solving the resulting 
equation for P,, we have 


(11) 
2 


2 = p'Py 


e sand YoOluNy 


<3 
— 


~ 


am adil Ad 


7. 
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Similarly, taking eq. (7) for » = 2 and substituting fer P, and P,, we have 
(12) P; = p*°Po 
and so on for P,, ete. 
It is easy to show by mathematical induction that 


(13) P, = p"Py, n= 1, 2, 2% 


Hence, the steady-state probability distribution P,, is known in terms of By: 


— i ~ 
te 


To evaluate P,, we use the fact that the total probability over the entire 
distribution is one. Thus, if we have no upper limit on the number of customers 
= whe 
: that can join the queue, we have 
a oo 
(14) =Pi=—1 


n=0 


ay which is called the normalizing equation. 
From (13) and (14), we have 
(15) Po(l+pt+p +...) =1 


Now, the infinite geometric series 1+p+p?+... is convergent only when 


p< 1. We conclude that the steady state exists only whenp <1. Assum- 
ing that p < 1 we have, from (15), 


(16) P, = 1—p 

so that, from (13), 

(17) P, = (—p)p*, (p a1, =, 1):2)..) 
which is a geometric probability distribution. 


The mean number in the system, which we shall denote by ZL, is given by 


(18) EAS ap: 


n=0 


= (1—p)(p+2p?+-3p?-+ ...) 


= A1—phl—p)s 
Thus, 


(19) PPE 
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r 
The mean r i 
: mM number in the queue (wait; for servi : 
| (waiting for service), denoted by La, is 


similarly given by 


(20) Liq = Ps (n—1)Pp 
=1 


n= 


Or a A* 


l—p = w{w—A) 


(b) The waiting time 


To study the waiting time, let us assume that the queue 


: discipline is 
first come, first served’. 


‘ “,* . . 

Then, the mean waiting time in the queue and 
in the system, denoted by W, and W, respectively can he fo 
arguments. 


und from simple 


Let us first take Wy. The customers found in the queue at any time 
are those that arrived during the waiting time of the customer at the head 
of the line. Therefore, the mean number of customers in queue is the mean 
number of arrivals during a mean waiting time, that is, 


(21) Lig = AWs, 
whence 
(22) Wa = LglA = Alfu(p—a).. 


By similar reasoning, we have 


(23) Lia, 
giving 
(24) W = L/JA = 1/(u—A). 


Now, to find the distribution of the time spent in system (w,), let a typical 
customer arrive to find n customers in the system. Then, for w, to lie in 
(w, w+dw) all the (n+-1) customers including the new arrival himself should 
complete their service within time w+dw. Thus, the waiting time density 


f(w) is given by 
(25) f(w)dw = Plw < ws < w+dw} 


@ 
= 2 P,:P{n services in time w}-P{1 service in dw} 


a tw) ree 
— 3 (1 —p)p"- (ut 3 
n= 


Ae : pdw 
> l 


4 
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oO 
= (1—p)we"Vdw & (npw)?/n ! 


n=0 
= (1—p)ywe"” dw.e?” 
= (1—p)ue—"1-9)#" dw 
3 Ut 
i i : as an exponential distribution with mean 
Hence, the time spent in system has an expon 
= 1/{(1—p)u} = 1/(u—A), which agrees with (24) 
istri i raiti ime in the queue, wg, we see that 
To find the distribution of the waiting time in the q ; 
(26) P(wg = 0) = Py = 1—p 
i ity fi i - the non-zero values of w 
and, arguing as before, the density function for the n uy 


is given by 


(27) folw)dw = Plw < wy < w+du} 


kad . . = = 5 . 
= & P,.P{(n—1) services in w}.P{1 service in dw} 


n=1 
co (uw)nle-Hw 
= 3 (apie » pdw 
caf ( p) (n—1)! 


= (l—p)pye-*” dw.e4?” 
= pu(l1—p)e-A-9) 4" dw 


From (26), (27) we see that the distribution of W, consists of a discrete atom 
at zero and a continuous portion which is a truncated exponential distribution 
over (0+, 00). The mean waiting time in the queue is 


(28) Wq = 0.11—p)+ f e.pu(1—p)e-—P)mo dw 
0+ 


p A 


~ WI= 2) ~ HA) 


(c) The busy period 


The distribution of busy period, which involves time-dependent proba- 
bilities, will be discussed in the next chapter. However, the mean duration 
of a busy period can be obtained from simple arguments as follows. 


The server's entire time is made up of alternating busy and idle periods; 
so the number of busy and idle periods in a long period of time 7 can be 
taken to be equal. Now, an idle period is the time from the instant the server 
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becomes idle until the next arrival which, according to the 
of the exponential inter-arr 
inter-arrival time. 


Markovian property 
ival time distribution, is distributed as ti 


Hence the mean duration of an idle period is 1/A. Now, 
the probabilit d » server is j : r ; 
if lity that the server is idle is Pj = 1—p, Therefore, over the 


time 7’ the server is expected to be idle for a duration T(1 
equivalent to TA(1 —p) idle periods. 


up of T'A(1—p) busy periods, and the 
(29) 


he entire 


—p), which is 
Thus the server’s busy time Tp is made 
mean duration of a busy period is 
Tpl[TAL—p)] = 1](~—A) 

It follows th 
(30) 


(d) Higher Moments 


at the mean number of customers served in a busy period is 


b(u~—aA) = (1—p)3 


In the foregoing discussion of the queue length, waiting time, etc., we 
have kept the analysis up to the derivation of the mean. which is the first 
moment, about the origin, of the distribution in question. However, it is 
often easy to evaluate higher moments like the second, third, fourth moments 
about the origin or other values, and so on, by using their definitions or by 
successive differentiation of the Laplace transforms or generating functions 
as the case may be. We give below derivations of the variances. 


(i) System size. Let N be the rv. denoting the number in system. 
Then, the variance of N is 


(31) Var[N] = ELN?]—(E(N}) 


where H[X] stands for expected value of X, and E[N] = L, given by (18), (19) 
Also, by definition, 


(32) E(W2] = = ntP, 
n=0 
= (1—p) 3 np" (on substituting for P,,) 
n=0 
= Aes p 2p?) 
a Laem er 
SU ae ie 
1—p _ (1l—p)? 
Thus, (31) simplifies to 
: p* 
fant =e a 
(33) Var[N] EPR ET. 
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is i i actorié ont, of order 2, of the system 
It is interesting to note that the factorial moment, of order 2, b 


size distribution is 


(34) ELN(N—1)] = ELN?]—EW] 
ee 
is 


ii J rj > and number 
(ii) Queue length. If Ny, Negery denote the number in queue 


in service, respectively, we have 
(35) N= Ne =e Nerv 

and from the description of the system it is clear that Ng and Ngery are inde- 
pendent. Hence, 

(36) Var[N] = Var[Nq]+ Var[ Necro]. 

Now, the r.v. Nsery has the probability distribution 


f P(N serv = 0)= Py =1—p 


(37) 
L P(N serv =1)=p 


whence it is easy to show that 
(38) Var[N sero] = p(1—p) 
From (35) and (38), we have 


pega ite Ce et ee 
(39) Var[Ng] = 1—p' (1—p)? pl p) 


(iii) Time in system. Since ws, the waiting time in system, has a negative 
exponential distribution with mean 1/(~—A), (see (25)), its variance is given by 
(40) Var[ws] = 1/(u¢—A)? 

Incidentally, it is easy to see that the r-th moment of ws about the origin is 
oo 
(41) f= J wtf(w)dw 
0 
= r!/(u—Ay 
f(w) in (41) being given by (25). 
(iv) Waiting time in queue. Proceeding as in (35), (36), we have for 


variance of wy, the waiting time in queue, 


(42) Var[w,| = Var[ws]—Var[v], v being the service time 
ee 1 
los ee a 


Note that (42) can also be easily derived from (27). 


Sec. 1.6 | 


POISSON QUEUES I L5 
1.6 State-Dependent Parameters 


We now consider the gs i 
consider the same queueing model as in Section 1.5 except that 


the mean arrival g service 7 i 
on wrival and service rates A and # will now be taken to depend upon 
16 state 4 > system ¢ ; j 
aay of the system and, therefore, will be denoted by A. and pi 

respectively } i iff ‘ “* : 

pectively. By choosing different forms of the functions A, and p,, we 
. £ in ; 1 ; ‘ 4 
shall obtain various interesting : 


exponential service times, 


models of queues with Poisson arrivals and 


Proceeding as before, eq. 1.5 (4), (6), now become 


1 d P 
( ) di n(f) 7 —(An+e,)P(Q+A,1P,,. Wb) +A nik p i (4), (n a 1, 2, ea 


‘ d 
(2) qe Poll) = —AgPold)+-psP y(t) 


Hence, the steady-state equations are 
(3) 0 = —(A,-+4,)P +A aan ata avs (n= 1,2, vee) 
(4) 0 = —AyPo+ MP, 
From (4), we have 
. A 
(5) Py = — Py 
fy 

Writing eq. (3) for n = 1, substituting for P, from (5) and solving for P,, 
we have 

AA 
= ee & P 

Pifle 


(6) P 


and so on. 


By mathematical induction, we can prove that 
(7) Po Stat p. 


In the general case when there is no upper limit on n, we have the normalizing 
equation 
Ay 


(8) Pol 1+—2+ 


AoA Te es AoAy Ana 2 a = ] 
By Male, a fge bn : 


For a steady-state solution to be possible, the infinite series inside the brackets 


in (8) should be convergent. Let it have a sum S; then 
(9) P, = 8 
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beady -sté tion 
From (7) and (9) we have the steady-state solutio 


AoA Ay -j S7 
Pyflg: fen 


a a Pap dy Byron 
(10) Par 


where an empty product is taken as unity. | e 
We proceed to consider a few particular cases of the model. 
(a) The queueing system M/M/1/N | | ite 
To obtain the solution for the queueing edi ie abs hier - ee 
oy rates A and yp, respectively, where no more than £ 
to enter the system at any time, we set 
(11) An =A(0Kn< WN) 
=0 (n> N) 

faa (@— 1h, 2h.) 
Then, from (8), we see that S is the finite sum 
| 5 1—pN11 
(12) S=liptet..te¥ = Fo, Sg 

i ‘ince S is finite, we conclude that the 
mere! fe Pr ee i ‘3 Substituting for S in 
stea f 2X ISTE 
(10), we have 
(13) Py, = (1—p)p2/(1—p* 4), p #1, n= 9, 1, 2, ..., N 
When p = 1, (12) gives S = N+1 and, therefore, P 
(14) Py, = 1f(N-+1), p=1, n ‘ 40s Ee eae ee 
which can also be obtained from (13) by letting p — 1. 


in the 3m is 
The mean number in the system 


N 
= ¥ nP. 
(15) L an a 
= PIL—(N + 1)p¥ + NpX*}/[(L—p)1—pX™)], (p # 1) 
= Ni2, (9 =) 
The mean number in the queue is 
; 1)P 
(16) Iq = 2 ee ) n 


= p{1—Np¥4(N—1)p%[(1—p\(1—p")], (p # 1) 
= M(N—1){°N+1)}, (p = 1) 
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b) The queueing system M /M|elao 
Let there be ¢ severs (channels), the ser 
©xponentially distributed with me 


with mean rate A, 
(17) 


vice times at each of which are 
an. service rate / and let arrivals be P: 


oisson 
» solution for this case, 


To obtain the we get 


A, aN for all n 
fr=ne forn<e 
= forn>c 
The justification for the 
the probability of 
one of the 


above substitutions is that 
a departure during (t, t+ Af) i 
nm Servers c 
probability, ig nyt. 


When # servers are busy, 
8 the probability that any 
ompletes service which, by the addition theorem of 


In this case, we have, from (8), 


A A2 Ach Ae A A2 
18 S = _ a - i 
(18) eo } 2 2 boat { 


+ {_ eek See | \ 
Cw cle teat t...| 


For a c-channel System it is customary to write the traffic intensity 
(19) Pp = Alc) 
Thus (18) gives 
cp)” cp) (epyje ; 
(20) S = 1+¢ep+ at Font GO (+) {1+p+p*+...} 


c~1 
= & (cp) fj | +{(ep)¢/e }1—p), provided p < } 

4=0 
The series in ( 


20) is convergent only when p < 1, 


conclude that the steady state exists only when p < 1, 


Using (17) in (7), we have 


(21) P,, = Py.c%p"/n ! for n <¢ 
= Py.c’p"/c! forn > c 
where 
c-1 
(22) z 


Po = 89 = [2% {coli Y-+{epyee ! (1—p))}1-. 


It is easy to sce that the probability of having to wait is 


(23) P(n > c) = Pete!) & pn 


nec 


= Pylep)ei{e | (1—p)} = Pel(1—p) 


1.é., when Ajew <1. We 


a swAoid w 3 


puejs ue!pul J 


ang uolunr ] nie 
sy QaDIAIZGS o1qNd as fe SPIE 
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The expected number waiting in the queue 16 


(24) Liq x (n—o)P, 


Pot(op)*/e I x (n—cjp™" 
= Pyplep)'e | —py"s 
AWA/p)oPo 
(c—1) ! (¢“—A)’ 


; Vaee ; , an number of busy 
The mean number of customers being served O1 the mean 8) 


channels is 


e~J bd 
Py 37 ) > >» 6 
(25) Ls = a. nP 4 Bye cP, 
C 7) 
. ’ ) 
> nF, | >» cl n 
nel not-1 


c (p)* Fp of 
Li WED ra x ¢.P, ol 
an NM! geet #3 


0 
cp = Alp 
Hence, the mean number in the system is 
(26) L = Let Ly = Lg t(Alp) 


The mean waiting time in queue is 

(27) Wa = Lq/A = MAlp)PPo/{(c—1) 1 (cu—A)?*} 

The mean time spent in the system is 

(28) W = Li/A = Wg+(1/p) 
(ec) The system M/M oa 
This is the system where every arrival is able to find a server and is some- 


times called the ample-server system. Possible applications are self-service 


systems where the customer carries the server with himself. In such systems 


ea 
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no suecing ja i , = 

19 queueing is involved and some authors like Conolly (1975) have pointed 
out that they should not be called queueing systems. Here we have 


6 
(29) A,, A for all n 
fy nh, n==%,1,2,.,. 
Therefore, 
Py, Lata Alo? 
(30) S 1-+-(A/y)+ ~ 1 = g/t 


Note that the series is convergent and, therefore, the steady state exists for all 
values of A/y. It is easily seen that we have 


(31) (A/p)re~*”* 


n ' 


P 


= 


so that the number in the system has a Poisson distribution with mean (Ajy). 
Thus 


(32) LD mAfh 


(33) W=LIA=1 fp 


which is otherwise obvious: since there is no queueing, the mean time in 


system equals the mean service time. 


(d) Queue with balling 


Let us consider the M/M/1 queueing system with mean arrival and service 
rates A and yw, such that when there are n units in the system, an incoming 
customer joins the queue with probability 1/(1+-n). So we take 


(34) A, = Ajf(a+1), » = 0,1, 2, ... 
Li for all n 
Then, 
r r\ 2 , 
(35) 8 =14(rjpyt MAY 4... = etm 


as before. Hence, in this case also, the number in the system has the Poisson 
distribution, given by (31). It follows that the mean number in the system 
is A/u, as before. But in the present case queueing is allowed, and we have 
the mean number in the queue 


oO 


(36) Lg = = (n—1)P,, = (Alp) t+e*"—1 


nol 


feos 


ne ueIpUl se ! 


me weve ie Ve ee 


20 


AN INTRODUCTION TO QUEUEING THEORY 


(e) The system M/M/c/c (loss system) 


This is the system with ¢ channels, where no queueing 18 allowe 
that the maximum number allowed in the system Is ©. 
this model is a telephone exchange W ith ¢ trunk 


ing subscribers who fail to get a free line. 


d, so 
An application of 
lines with no facility for hold- 


Taking A as the mean arrival rate 
and yw as the mean service rate for cach channel, we have 


(37) ieee forn <¢ 


= 0 forn > c 
fn =me fornge 


Therefore, S is the finite sum 


A/ 2 Alnye 
(38) 8 = iia Aju? lw) 


>] Toeee TF 


c! 
Thus, the steady state exists for all values of A/u, and we have 
(39) 


p, = WEE gO <n <0) 
‘ nN! 


where JS is given by (38). The probability that all the channels are busy and, 
therefore, the incoming calls are lost is given by 


(40) 


faye) ped Sa 


(40) is called Brlang’s loss formula and is still in use in telephone traffic work 


(f) The machine interference model 


Let the customers come from a finite source of k members, such that the 
probability that any member will call for service, independent of other 
members, in the interval (¢, t+At) is AAt+-o(At), and let there be a single server 
with mean service rate #, the service times being exponentially distributed. 
Then, we have 


(41) 


A, =(k—n)A forn<k 
sally forn > k 
i, =p Ot i= Ney 
Therefore, 
(42) 


S = 14 MAlp) + Mb--1)(Alp)?+...-+h V(Alu)E 
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and 
4 
(43) Py = 1/9 
44 oe HG) ; 
(44) P= bk—1)...(k—n-+ 1\(AJu)P,, (w= 1,2, ..:,b) 
Possible applications of this model 


are situations where there are k ati 
e ore p automatic 
machines under the charge of a si 


ngle operator and the machines fail and 
ai gy ac an 
require the operator’s attention from time to time 
industrial workers re 


or situations where | 
crib, ete. 


quire service at a facility like a grinding machine, tool 
(g) The brand-share model 

s} 7 ol Pa Jatt 

Consider a market consisting of N consumers, out of which at any time 
” consumers are using a particular brand, say brand A, of a product. Let 
the consumers switch into and out of brand A in a Poisson fashion with rates 
A and 1, respectively, so that when the brand-share of A is n, we can write 
(45) A, =(N—n)A forn < N 

—0 


for n > N 


Ly = MM forn < N 
Thus, we have 


(46) 


N \N 
S=1+(', )Am+(, )Almt+...+ (5, ale 
= 4AM} 
whence 
(47) Py = 1/9 = {p]A+m 


(48) - 


Clearly P,, is the probability that A has a brand share a, while 
(49) Py = {A/(A+p)}” 


is the probability of market saturation by brand A. It may be noted that 
P,, follows the binomial probability distribution : 


(50) 


where 


(51) p=Al(A+z), ¥=1—-p 


See AE LS 
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Hence, the expected brand share of A is 


(52) E(n) = Np = NAf(A+ /) 
Ex. 1. 


A television repairman finds that the time spent on his jobs 


has an exponential distribution with mean 30 minutes. If he repairs sets 


in the order in which they come in, and if the arrivals of sets are approximately 
Poisson with an average rate of 10 per 8-hour day, what is the repairman’s 
expected idle time each day ? How many jobs are ahead of the average set 
just brought in ? [By permission from Sasieni, Yaspan and Friedman (1959)] 


We have an M/M/1 queueing system with the mean arrival rate 


0 ‘s 
al = = i sets per hour 


and the mean service rate 


ft = 2 sets per hour 


Therefore, 


Hence the repairman’s (server’s) expected idle time per day 


3 : 3 ; 
= = x duration of a day (i.e., working hours) 


— x 8 hours = 3 hours 


Also, expected number of units in the system 


= 4 
ee ee p 5/8 5 


Pye 
ip ea i 3 jobs (sets) 


Ex. 2, A bank has two tellers working on savings accounts. The 
first teller handles withdrawals only and the second teller handles deposits 
only. It has been observed that the service time distributions for both 
deposits and withdrawals are exponential with mean service time 3 minutes 
per customer. Depositors are found to arrive in a Poisson fashion throughout 


the day with mean arrival rate 16 per hour. Withdrawers also arrive in a 


Poisson fashion with mean arrival rate 14 per hour, What would be the 
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effect on the average waiting 


ba) 


time for depositors and withdrawers if each 
teller could handle both withdrawals and deposits 4 What would be the 
effect if this could only be accomplished by increasing the mean service time 


Qe we eat qe Bae ees * 4 . P r 7 6 
to 3-5 minutes? [ By permission from Sasieni, Yaspan and Friedman (1959)| 


Originally we have for depositors, A = 16 arrivals per hour, w = 20 


services per hour, so that p = Alp = 4/5 < 1 (steady state exists). 


g uelpu Jo ne 


nat 


: . A 16 1 
Therefore, i : Sg: eee = 12 mi 
qa alia) 20(20—16) 5 hour = 12 minutes 


puke} 


And, for withdrawers, A = 14 per hour, ~ = 20 per hour, p = 7/10 <1 
Therefore, WV, == is hour = 7 minutes 
a 20x 6 


In the second case, we have an M/M/2 system with A = 16414 = 30 
per hour, = 20 per hour for each server, so that 


P—Alzn— 2 <1 (so steady state exists) 


Therefore, from 1.6 (22), taking c = 2, we have 
p= [3S (5) ea) 


Hence, from 1.6(27), 


9 

20x 

gene BALE? Bs 
(i= (2u—A)? 0 100 


AN adlAlIS anand uolun ty Sple 


Ore. 
x 7 hour = 3-— min. 


In the third case, when the mean service time increases to 3°5 min., we have 


i services per hour, and A= 30 per hour, as before. Therefore, 
ere 


=A/2z = 7/8<1. Thus, we have 


ee ae ee 
Py = peer | = (14744 = 1/15 
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} whence 
120 9 
ai (7/4)? a 343, Ll _ min 
= sgt * iG 10ur = 30 OS 80 : 
(80) 


1.7 Alternative Solution Methods 
We give below two alternative methods for the solution of the set 
of equations 1.5(7)-(8) for the steady-state queue length probabilities P,, 
in M/M/1. 
} (a) Difference-equation method 
: The equations can be written 
(1) PP ayo—(At Pasi tAP, = 0, n= 0,1, 2,... 
(2) pP,—AP, = 0 
Now, (1) is a difference equation with constant coefficients. Using the 
difference operator H defined by 
CP = be aigh i Lee 
we can write (1) in the form 
[uB*—(A+p)E+A]P, = 0 
so that the characteristic equation is 
pB*—(A+p)B-+A = 0 
or, 
(4H—A)(H—1) = 0 
which gives 
EH = 1, p, where p = A/u 
Hence, the solution of (1) is 
(3) P, = A+Bp" 
where A, B are arbitrary constants. 


Now, substituting from (3) in (2), we have 
A+Bp = p(4+B) 
which gives A=0 and hence B=P,. Thus we have, from (3), 
(4) P, = Pop” 


which is the same as 1.5(13), and we can proceed further as before to show 
that P, = 1—p(p <1), which completes the solution. 
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(b) Generating-function method 


Let us define the generating function 


(5) P@) = E Pawn, jz] <1 
n=0 


Multiplying eq. 1.5(7) by 2, summing over n = 1, 2,..., and adding 1.5(8) 
we have st 


(6) O= —A+m) E Part yP tra Py mtn S Pog 
n=0 ~<a aol atom 

or, 

(7) 0= —(A +H) P(2) + WP y+-AzP(z) +(u lz) P(z)— P,} 

or, 

(8) {Az2— (A+ p)z-+ po} P(z) + (z—1)P, = 0 


This yields 
(9) Pe) = ._Mi—2)Po 
AA—(A+-p)z+p 
The denominator of the right-hand member of (9) can be factorized as 
Ae —(A+-p)z-+- pe = (Az—p)(z—1) 
so that (9) can be written 
iP ma 
(10) P¢) = zk = P,(1—pz)1 = P, = (pz), since pz < 1 
, - n=0 


Recalling the definition (5) and equating the coefficients of z” on the two sides, 
we have 

(11) P, = Pyp 

Py can be evaluated as in 1.5. Alternatively, we may note that P(1) = 1, 
and then (10) gives P, = 1—p, as before. 

Ex. 3. Find the steady-state probability P,, for the finite-waiting-space 
multichannel queueing system M/M/c/N, where the number of channels is 
¢ and the maximum number allowed in the system is N (N > c). Also deduce 
expressions for L, Lg, W, Wg and the expected number of busy channels. 

In the notation of Section 1.6, we have 

A= A when OC n< WN 
Wh) when n > N 
fin = ye ~=0Owwhen n << 


=cu whenn>c 
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Therefore, writing A/cu = p, we have as in 1.6(21) 
P, = Poctpt/n!, (0 n< 0c) 
= Pyrprie!, (ec<n< WN) 
where (as in 1.6(18)) 


: e1 (cpt, (cp) 8° 
g = Py =a OY OF E 


SF leper)’ VPN eP when p # 1 
> | 


© geo 9) VTS ese 
cK Osi his een ean 
j-0)! c! 


Note that the steady state exists for all values of p. 


Now, the mean number in the queue (waiting) is 


N 
Iqg= © (n—c)P, = Pylepy'le}} E (n—c)pr-< 
1 n=ct1 


n=ct 
= {P,(cp)¢/c }}{p-+2p?+ 3p%+ ...+(N—c)p%—*} 


Py (cp)p a ae 
=— 7, [1—p¥-*—(1—p)(N—e ] 
Sia ame (1—p)(W—e)p 
The average number being served or the mean number of busy channels is 


ci N 
L, == = nP,+ = cP, 
n=0 


n=C 
= opPo [ 3 (ep) in—1) H+ pe] 
= cp 0 | 2 ko) n : c! eS 
- C—1 } cp\e N-e-1 
= pP, | = (ool! 4A zp 


t n=0 


= opPy [ Py 1—SPK . pare | = ep(1—Py) = (Aln)(1—Pw) 


Therefore, the mean number in the system is 


L = Lg+Lg = Ly-+(Au)(1—Py) 
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The effective arrival rate at the system is A’ = A(1—P 


Nn). Therefore, from 
Little’s formula, 


W = LIN = Wet+(1/u) 
Wg = Lg|r' 
Ex. 4. At a port there are six unloading berths and four unloading 
crews. When all the berths are full, arriving ships are diverted to an over- 
flow facility 20 miles down the river. Tankers arrive according to a Poisson 


process with a mean of one every 2 hours, It takes an unloading crew, on 


the average, ten hours to unload a tanker, the unloading time following an 
exponential distribution. 


(a) On the average, how many tankers are at the port 2 
(b) On the average, how long does a tanker spend at the port 2 


(c) What is the average arrival rate at the overflow facility ? 


= 1 
Here we have an M/M/c/N system with c = 4 and N — 6. Also, A = z 


arrival per hour, and p = = service per hour. Thus, cp =A/u=5, and 
p = 5/4. Therefore, as in Ex. 3, 
UPa= rt a Toy 7 te ot ae +e x 
= 138-62 
Therefore, 
P, = 1/138-62 ~ 0-0072 


6 
Py = Pyetp/c | = 0-0072 x 256 x ( ) [24 — 0-293 


we | or 


Hence, 


(a) The mean number of tankers at the port is 


L = Lg+(Ajp)(l—Py) 


Now, 
Py(cp)*p ” is 
= —0\PTP 11 pN-e_(1—p)(N—c)pN-e 
0.0072 x 54x (5/4) Bi J gee Aiea 
ee NIL pO. 2 egg 2O |) airgap 
24 x (1/16) [) iwia 16 


Bie 


“yeaa 


r ih eae ] ; eC. 1 
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Thus, = 0-820--5(1—0-293) = 4-36 tankers approx. 
(b) The average time spent at the port is 
W = Li{A(1—Py)} = 4:36/(4 x 0-707) = 12-33 hours 
: Thus the average time spent in waiting before unloading is 
i Wq = W—10 hours = 2:33 hours 
} (c) The average arrival rate at the overflow facility is 


j APy = $x 0-293 = 0-15 tanker per hour 


PROBLEM SET 1 


: 1. Arri vals at @ el phone booth are considered to be Poisson with an 
: V t t e 
average time of 10 minutes between one arrival and the next. The length 


of a phone call is assumed to be distributed exponentially with mean 
3 minutes. 


(a) What is the probability that a person arriving at the booth will have 
to wait ? 

(b) The telephone department will install a second booth when convinced 
that an arrival would expect to wait for at least 3 minutes for the phone. 
By how much should the flow of arrivals increase in order to justify a second 
booth ? [By permission from Sasieni, Yaspan and Friedman (1959)] 


[(a) 0-3; (b) new arrival rate one every six minutes] 


2. A supermarket has two girls ringing up sales 


at the counters. If 
the service time for each customer is exponential with 


mean 4 minutes, and 
if people arrive in a Poisson fashion at the rate of 10 an hour, 
(a) What is the probability of having to wait for service 2 


(b) What is the expected percentage of idle time for each girl ? 


(c) Find the average queue length and the average number of units in 
the system. [By permission from Sasieni, Yaspan and Friedman (1959)] 


[(a) 1/6; (b) 6639/; (c) 1/12; 2/3] 
3. In a railway marshalling yard, goods 


trains arrive at the rate of 
30 trains per day. 


Assuming that the inter-arrival time 


follows an exponen- 
tial distribution and the service-time di 


stribution is also exponential with an 
average of 36 minutes, calculate the following : 


(i) the mean queue size (line length). 


(ii) the probability that the queue size exceeds 10, 


[P/(1—p) = 8 trains; p — 0-056] 
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4. At what average rate 

to ensure a probability of 0. 
than 12 minutes ( 


must a clerk at 
90 that the custome 
a) betore service. (I>) 
is only one counter, 
average rate of ] 
tial distribution. 


# supermarket work in order 
r will not have to wait longer 
including service ? Assume that there 
at which customers arrive oisson fashion at an 
5 per hour. The length of 


[By permission from S 


ina P 
service by the clerk has an exponen- 
asieni, Yaspan and Friedman (1959)| 
[(a) 2-48 minutes per service: (b) 2-26 minutes per service) 
5. Problems arrive 
average rate of five per d 
man waiting to get 
being solved. 


at a computer centre 


in a Poisson fashion at an 
ay. 


The rules of the computer centre are that 
his problem solved must 
[f the time to solve 
distribution with me 


any 
aid the man whose problem is 
a problem with one man has an ex 
an time of 1/3 day, and if the 
inversely proportional to the number of people 
approximate the expected time 
[By permission from 


ponential 
average solving time is 
working on the problem, 
in the centre for a person entering the line. 
Sasieni, Yaspan and Friedman (1959)} 
{8 hours) 


planning to install telephone booths in a 
It has established the policy that 


wait more than 10 per cent of the time 
for use is estimated to be Po 
average phone 


6. A telephone company is 
new airport. ® person should not have to 


8 he tries to use a phone. The demand 
isson with an average of 30 per hour. The 
call has an exponential distribution with 


@ mean time of 5 
minutes. How many phone booths 


should be installed 2 [6] 
7. Show that for the k-Erlang distribution given by 1-4(18), the me 
variance and mode are, respectively, 1/A, 1/kA®, and (1— 1/k)(1JA). 
8. For an M/M/2 queueing system in steady state 
ALg = pPy 
(A+y)P; = AP, +2uP, 
(AF SBP p= AP. 2hP aay (nm = 2, 3, ...) 


Hence prove that P, = 2p"P, (n= 1, 2, ...), 


an, 


, Show that 


where p=A/2u, and 
Py) = (1—p)/(1+p), and deduce that L = 2p/(1—p?), L 


1q = 2p /(1—p?), and 
that the probability of having to wait is 2p" /(1+-p). 


9. Find the steady-state probability P,, given that 


((N—n)A, O0<n<N 
A, ae 4 


L 0, n 


Nt, 1 
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Describe a practical situation represented by the above model. 


AT 


[P.=(_) Alp, 0<n<e 
nm’ 


T 


= ( ) (mle ten-ey(AjuymP,, oS n< N 


10. In Problem 3, assume that the yard can admit 9 trains at a time 
(there being 10 lines, one of which is reserved for shunting purposes). Cal- 
culate the probability that the yard is (i) empty, (ii) completely occupied, 
and find the mean number of trains in the yard, {$67 7% 

[0-26; 0-02; 2-4 trains] 


11. For an M/M/1 system, show by simple arguments that Wa = Lp, 
and W ==(£+1)/u. Hence use the expression for Z to deduce those for 
Wa and W. 


12. At a railway station, only one train is handled at atime. The 
railway yard is sufficient only for two trains to wait while another is given 
signal to leave the station. Trains arrive at the station at an average of 
6 per hour, and the railway station can handle them on an average of 12 per 
hour. Assuming Poisson arrivals and exponential service distribution, find 
the steady-state probabilities for the various number of trains in the system. 
Also find the average waiting time of a new train coming into the yard. 


[P, = 0-533, W, = 32 minutes] 


Cy 


_ SEE! oo at ca ROTOR imiReindsrmaat toe 1 By 


CHAPTER 2 


POISSON QUEUES II 


2.1 Introduction 

In Chapter 1, we obtain the steady-state solutions of various queueing 
systems with Poisson arrivals and exponential service times like M/M/\, 
M/M/c, ete. In this chapter we shall discuss the time-dependent or transient 
solutions of similar queueing systems. 


2.2 The Queueing System M/M/1—Transient Solution 


As we saw in Section 1.5, the probabilities P(t) that there are n units 
in the system at time ¢ are given by (Cf. 1.5(4), (6)) 


d 
(1) dt P.( a —(A+P)P a(t) FAP, 1(t)+ 2B a sr(t), = IL, 2, ee 
‘ d 
(2) di P(t) = —AP,(t)+-~P,(t) 


Define the p.g.f. (probability generating function) 
o 
(3) G(z, t) = X P,(t)z 
n=0 
which clearly converges inside and on the unit circle |z| =1. Multiplying 
(1) by 2", summing over n and adding (2), we have, on simplification, 


(4) : Gz, t) = —(A+p—Az- L2G, t)+-w(1—1)z)P,(t) 


Let the system start at time ¢ = 0 with i units in the system, so that 
(5) P,(0) = din 
f I, whenn =i 
0, whenn 4i 
Jin is called the Kronecker delta. Thus, we have 
(6) Giz, 0) = 2! 
Let g(z,s) be the Laplace transform (L.T.) of G(z, t), defined by 


(7) g(@, 8) = f e#tare, t)dt, Re(s) > 0 
0 
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and let the L.T.’s of other functions be similarly denoted by the corresponding 
lower case letters. From (4), on taking L.T.’s, wo have 


(s-+-A-+-—-Az—p/z)g(z, 8) = 2+ (1 —1/2)pol8) 


(8) 
or, 
zitl— 4(1—z)po(8) 
(9) He. 8) = a ge Aa 
Since g(z, s) converges inside and on the unit circle for Re(s) > 9, the 
zeros of the denominator inside and on |z| = 1 must coincide with the 


. . gar , ’ , inator are 
corresponding zeros of the numerator. ‘The zeros of the denominator ar 


a,(8) = (2A)-[s+A+p—V(s+-A+p)?—4Az] 


(10) 
(11) G9(8) = (2A)“Ys+-A+pu+ V(s+A+u)?—4Az] 
Now, let 
(12) fe) = (8+A+p)z 
g(z) = —(Az*-+-) 
Then, on the unit circle |z| = 1, 
(13) If@| = |st+A+p| > [A+z4| = 9@) 


Hence, by Rouché’s theorem, f(z) and f(z)+-g(z) have the same number of zeros 
But f(z) has only one zero, and so f(z)+g(z) also has 


inside the unit circle. 
This zero must be #,. Thus we have 


only one zero inside the unit circle. 


git! 
were Pde) = Tay 
Therefore, 
ziti_(1—z ay * 
(15) Ie a= 1—, 
A(z—a1)(2— Gy) 
or 


{| 


241(1 —a,)—(1—z)a,!*1 
J 


g(Z, 8) Aate(z—,)(1—z/,\(1—or) 


(2643. i ait) —2y(2' — at) 
Aoto(7—O)(1 —2/ 4y)(1—ar,) 


_ tae! t+... +04) tay (2! + o2'*-+... +. ,-1) 
Nato 1 —2/G%q)(1 —o) 


tae... +0f)(1—ay) +aft? 


Note 1—z]arg)(1—ary) 
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or, using binomial expansion, 
(16) gfe, a) = Glas bs es 
Ue, ¢€ vA Zz J ” ale et a," $3 
: Ads (2'-| a ee ies | aj) & (2/a@_)* 4 2X (z/a,)¥ 
kum Ad,< I 44) k= ro 


Note that |2z/c,| < U3 


| No, P,{8), the L.T. of P,(t) is the coofficient of zn in g{z, 6). The contri 
ution to the coefficient of 2” from the second term on the right-hand side 


of (16) is 


ay t 
m4 aft 


(17 re (chs 
AaA(I— oy) ~ Jag +e +ai+...) 
i otis wo 
= ee? Ce — 
A (@ya,)0*3 - a 
1 A n-+4 tvs) 
: Ss al 
A pL k=nti+a | 
] ( A perk. Ss ( ey” J : 
A\ p k weit \ A a” (nen, ee er eee 
And, the contribution from the first term is 
ut J ay f [\\m 
(18) >» : s __ (KIA) 
m= (i~n)* Aa, ag tm m rh n)+ Ans 2m+1 
where 
(19) (x)*+ = max (2, 0) 
Therefore, 
t 4 rs 
(20) P,A8) = a (u/Aym » : )" >» (m/A)* 
A lime (i—nye EE ] petine od 


Now, from Appendix E, eq. (18), we have, for the inverse Laplace transform, 


(21) "(8+ 4/s*9—a? )-” ] = va-” tf] (at) 
where J, (z) is the modified Bessel function of the first kind and order v, given 


by 


22 ae eC 
ie He) = BEIT e+e) 


ld SWINTRODUOTION TO QUANT HIN PHMORY 
Pharntin 
(any : } , yp (AE (UA)A Ata 4 Vo) MY 1/y(B 4 Wit) 

AVE Cy ALAet Wipe wApe A) 
flanwe, fron (80), tolling fhe inverse t Aplin Wvanslormnea, we have 

poke iA tii : from bel Dpiael ' 
(B4) D(H) [ (¥ , | Fyody pny (Oy Aft #) 
\ ji f i f 
teeth fh 
A , wed ' | we yf h 
| os iY VApt) 
7 homers \ / 


On using the well loiowi reourrenee relation 
(Uh) (Sv/eys le) Fy g(@) md yy g(@) 
aid ather some olyvioia ain plifivation, we have 


(Aju (4 hi] A)! i ij | (\ BIA)! WING day 


(U0) I’ (f) 
(1 Vy Alaye 7 {\ HIAYE LY 
hkeowypepe 
where /, Lie VAS) \valn, alioe for iatepral v, 
(#7) / te) /i(e) 


aa can be meen from (89), the lin TLL (YH) oni byes repilneod hy ie 


The Tuey Poviod Distribution 
Vor finding the diateiliition of Ditey period, we fmeninis fhat ah tine ¢ i) 


the ayatem aharte with the arrival of # customer whieh makes the timber in 


y Pr] 


the ayatem wv f and that there ia an absorbing barrier at the abate w if) 


Then, oloarly, 


(1) f(t) ’ (litany period « 


id the probability dietoibution finetion (D9) of the biey period duration 


f) M(t), win 


Vidor these adsunphiona, the 77(/) ave given by 


7] 
(3) di L(t) (Ad MP MAAR, GD AP yh), 1 1, B, 
: d ; 
(il) di M(t) (Ad LO) wll) 


it 


4) 
(4) i 


i Pil) phy 


Pilar Qui ma tt is 


May | | 


| 
mbiie thie onerahinnge Canebion 


(hh 4 
(et) amt 


we have fain (4) mel ( ') 


(A) il 
ft 


file, tp) LA H \ fil jiflm, ty HP 


OP(O), wed ete niobiii the teaneforne te 
y 


) 


Now, taking o/T a on howl aides 


CoOrPanpondiiig @inadl loibeva. we have 


(7 | 
} pA fA flee, a) jip la) 


or, 


(Mt) g(a, a) 
i | A | pt \ Hl 


TT 
‘he two gerow of tha donomninaoy ava given Wy ,OC1O), (11) A : ) 
t ith Hohn 
7 ' } ] 
7; Li ; 
F ¥ Thonioha'a theorem the hiineraho of the righthand member of (4 
} 


Vivtiialion fal 4, ‘White 


(1) fipiylo) tj 


Wrom (1), (4) and (0), the hiimy period denaliy funetlon 4 


/ 
(10) f(A) ome, WU) ome Y(t) wm ten, 


yp ees bp) lad As pte dA 


fi 
AA 
/ | t 
(A, A4 ft) | lay) At yy? bAyi 
; I ! / 
ii (% ity (AVHYE 
Way j ey 
h gelatine i? 
VIIA / (leap f) 


whore I, ) Ie (hia molified Hepaal finietion of the fret bind ail opley 4. ane 


we hive ised 8,9 (41) for evaluating the inyern Laplac tranaform 
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The Initial Busy Period 


The initial busy period may start with any initial number » = 7 in the 


system instead of x | as we assumed above, So, proceeding as above and 
taking G(z, 0) = 2!, we see that tho busy period density is 

, @—1/ yi 
(11) Silt) = &-(a;) 


| 
= e-(AtHe . i(4/fe]A)! r Ti(2/Ap t). 


Alternatively, (11) may be proved as follows, 

We know that the busy period duration does not depend on the queue 
discipline, The initial busy period under considoration starts with ¢ customers 
in the system; let these customers be numebered 1, 2,...,%. Let us first serve 
the customer ‘1’ followed by service of all customers that arrive during his 
service time, then followed by all customers who arrive during their service 
times, and so on until all such arrivals have been serviced, Then, we troat 
similarly customers 2, 3,...,7. Thus it is clear that the distribution of a 
busy period starting with ¢ customers is tho t-fold convolution of that of an 
‘ordinary’ busy period starting with one customer, Let tho corresponding 
probability density functions (p.d.f.) of these distributions bo, respectively, 
Ji(t) and f(t). Now, we have shown in (10) that 


(12) L(t) = ay 
Therefore, by the convolution theorem of Laplace transforms, 
(13) Lilt) = af 


from which (11) follows, 


2.4 The System M/M/co 

The steady-state solution of the service system M/M/co was discussed in 
Section 1.6, The probabilities, P(t), for n units in the systom at time ¢ for 
this system are given by 


d 
(1) di P(t) —(A-+-np)P,(t)4 AP,,_,(t)+-(n-4 uP 4(t), (n = 1, 2 Sep. 


es | 


i | 
(2) ap Poll) = —AP\(t)-+nP,(t) 


Define the generating function 


(3) Hz, t) = 5 P(t)", |z| <1 
=0 


n 
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Multiplying (1) by 2”, summing over n and adding (2) as usual, wo have 


OG © 


a Atm P OAL Py ety F (ntl) P. (t)zm 
ie nol Pom nt4 

where @ = Gz, 1) 
or, 

he aq og 

= = AG Zz oi eT y \_, : 7 

ot @ N—me Ze +AaGe, t+ Bz 
or, 
4) 0G : i ae: 
| op ~1l—z) — = —A(1—2)G4, t) 


my . . “ff . ‘ a . 
lo solve this partial differential equation of Lagrange type, we write the corres- 
ponding simultaneous equations 


dt dz dG 


(5) l  =MI—s) ~” —A ae 


The solution of the equation obtained from the first two members of the above 


equation is given by 

(6) (l—z)e~*? = ¢, 

and the solution of the equation obtained from the last two members is 
(7) Ge-Aimez — fe 

where ¢,, cy are arbitrary constants. Hence, the general solution of (4) is 
(8) G(z, t) = el4/Med{(1—z)e-*#} 

where ¢ is an arbitrary function, 


Now, if the initial number in the system is i, we have the boundary 
condition 


(9) Gz, 0) —— A 


From (8) and (9), we have 


at = ef4/M) (1 —z) 


whence 

(10) $(@) = (1-2)! exp [—(AJj)(1—2)} 

Using (10) in (8), we have 

(11) Ge, #) = oxp{(A/m)z]{1—(1—2)e-*8}! + expl—(A/u)(1—(1—2)e*4] 
= (q-+p2)' exp[—(Aq/)(1—z)] 
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where 
(12) 


Taking the coefficient of z” in (11), we have 


= enHt, q= 1—p 


n t (A/u)* —Ag/tpn—ki—nt+2k 
eo es AE * e7*d n q 
(13) P,() = 2 tigen hie # 


The mean number in the system at time ¢ is given by 


ag | MOE! eax, 2 fet ee 
(14) L => ae = p-- it == ni | (i hh Je 


And, the variance V of the number in the system at time ¢ is given by 


é 0G ‘ 
(15) Visa | “pee 
z=1 
which simplifies to 
A 
(16) V = tpg+ 2 


= (A/u)+-(i—A]p)e—"*# —te- Ht 


It may be noted from (11) that G(z, t) is the product of the generating 


functions of the binomial distribution 


(17) BAR cs BY @ gi-kpk, k = 0,1,..., 4 


‘ 


and the Poisson distribution 


(18) P(X = k) = e-*9/*(Ag/u)E/k !, & = 0, 1, 2, ... 


Thus, the distribution of the number in the system at time ¢ is the convolu- 
tion of these two distributions. This gives us an alternative way of writing 
(13) by taking the convolution of the two distributions given by (17) and 
(18). Also, the mean and variance of the number in the system at time ¢ can 
then be obtained by adding the means and variances of the above binomial 


and Poisson distributions, which directly gives the results given in (14) 
and (16). 
For the steady-state distribution of the number in system, we let t— 


in (11) and obtain 


(19) Ge) = lim Gee, t) = exp [—/p\(1—z)] 


oh 
~ 
i 
KE 
Be 
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which is the generating function of the Poisson distribution with m« n Aj 
y ICANN [iL 
IPs 


giving the steady-state probabilities 
20 ies 
(20) Pp. e MAI)" In | (n = 0, 1, 2: ree) 


which are as derived in 1,6(31). 


PROBLEM SET 2 

is 1. From 2.2(26) deduce the steady-state probability P,, by making t+ o 
[From the asymptotic behaviour of Iz), we have | ‘ : : 
(1) Iz) ~ exp(z)/4/2nz 
which is independent of k. Hence, for large t, 
(2) eCMtT (24/Ap t) ~ expl—(4/A— Vp) t/t 

+ O0ast4o 
Thus the first two terms on the right-hand side of “.2(26) vanish. Also. it 


is well-known that 


(3) expl(a/2)\(y+1/y))= & ytiye) = "S ytiucrs = : 
: C= oe YL (2)+- XL yk yx) t+ YD yk le 
k=—co k=0 a) a K\) k Rid Ifa) 


since I_y(x) = x(a) for integral k. 


Taking « = 2V Aut, y = Va, and using (1) and the fact that Alp <1, 
we see that the first and the third sums on the right-hand side of (3) vanish 


as ¢—>oo. Thus (3) reduces to 
@ alte = 
mal VEIT 2A t) = oxpl(/Ape t)(4/p]JA+/Aln)) = exp[(A+y)4] 
Taking m = n+i+2, we see from 2.2(26) that 
P(t) > (1—A/n)(Aly)™ as t+ co,] 


2. Obtain the distribution of an ordinary busy period in M [M/1 system 
and use it to show that the mean duration of a busy period is 1/(u—A). 


[Mean busy period = J pode, where f(t) is the busy period density fune- 
tion given by 2.3(10). Use Erdélyi et al. (1954a, p. 195, eq. (1)), viz. 
L& [1 (at)] = a*(s*—a?)-1/(54 38 @2)-¥ 
Alternatively, mean busy period = — f'(0), where f(s) = LG (t) = a, 
from 2.3(10).] 
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40 AN INTRODUCTION TO QuE linary) busy period is A+) |(H— Ay. 
e (ordinary 8} 
oO > that the variance of the (ore hosts pt a 
3. Show that *s the busy period, Var [B] = BLB")— (B(BI 
: stands for i F , eae , 23 10). oo 
(if B star E(B?) = i ep(ede, with f(t) given by ( N 
E(B] = 1(u—A), and EE : 


iq. (2), viz. 
e Erdélyi et al. (19548), P- 195, Eq. (2) a 
use w J = 5 9\—3/2 a/32—a? : 
Lit I,(at)] = a(s+v \/ 32—a?)(s?@—a@*) (s+ 
Z vy 
iid 7 f(s) given by 
Alternatively, E[B*| = f’(0), with f(s) g 
BB? = 2n/(n—A)?*.] ‘ sy period 
Ree cat a busy perio 
ay M/c queueing system, show that the dur a GF 2 . te a 
a va ae white all the c channels are continuously DUS) 
defined as the time 


the density function 


2.3(10), which simplifies to 


a. 


e-(atemte-1y /op JA T,(2y/cAp t). 
F riod is 1/(c“—A). 
how that the mean duration of such a perio / ; 
Hence show . ee eee 
2.4(17), (18) give, respectively, the a en mber 7 
5. Show that Eq. 2. , at time ¢ out of the initial nu 


‘ill in the system ager 
the customers who are still > Reeanen Oi 


b m zer a wh 


x 


CHAPTER 3 


NON-POISSON QUEUES 


3.1 Introduction 

In this chapter we shall consider 
distribution of inter-arrival times 
arbitrary distribution. 


queueing systems in which either the 
or that of service times is a general or 
As was pointed out in Section 1.4, the 
exponential distribution possesses the Markovi 
80 that in the case of an M/M/ 


the time since the last 


negative- 
an or the forgetfulness proy 
1 queue we do not have to take in 
arrival or the elapsed service time of the unit 
This is not the case with queues 
time distributions are differe 
next section. 


erty, 
to account 
in service, 
where one or both of the arrival and service 


nt from exponential. It is explained in the 


3.2 Imbedded Markov Chain 


To be specific, let us first take the case of an M/G/ 
Here the inter-arrival time distribution is expone 
property. However, the service 


1 queueing system. 
ntial, which has the Markovian 
time distribution being general, the 
bility of a service completion in time (, t+-Af) 
the service has already bee 


proba- 
depends upon the time for which 
n given. To overcome this complication, we study 
the system not in continuous time but at the 


discrete time points at which 
customers depart after service. Between any 


two consecutive departures, 
we have to consider only arrivals which are Markovian. 


Suppose we are 
considering the queue length process or 


the number in the system. The 
queue length at the departure epochs is said to constitute a M 


arkov chain 
imbedded in the non-Markovian queue length 


process in continuous time. 
The departure instants are called regeneration points and are 
knowledge of the state of the process at any rege 
previous history of the process irrelevant. 


such that a 
neration point renders 


In the case of the GI/M/1 queue with general arrivals and exponential 
service times, the arrival epochs are the regeneration points. Clearly for an 
M/M/1 system, all time points are regeneration points and the whole process 


in continuous time is Markovian. 


6 
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3.3 The System M/G/1 
Poisson arrivals at a mean 


m with 
and one 


with density function 6(v) 
from the system, ” = Lio) avs 
the system at time f, 


queueing syste 
eneral service time distribution 
of the n-th departure 
number of customers in 


Let us consider the 
rate A, ag 
server. Let #, be the time 
t,= 0. Let N(t) denote the 
and let 
(1) N,= N(t,+°) 

i.e., N,, denotes the number in the system left behind by the n-th departing 
customer. As explained in Section 3.2, the discrete-time process {N,} consti- 
tues a Markov chain imbedded in the continuous-time process {N(t)}- Further, 

s of this Markov chain, viz. 


it can be shown that the transition probabilitie 


(2) py PLN ns =x j\Nn = i}, (i, 2 0, 2 1) 


Such a Markov chain is said to be time-homogeneous. 


are independent of n. 
it follows that 


In fact, from simple probability arguments 
(3) po = ky (9 > 
pis = Mytay G>1j2t-) 


= 0, otherwise 


(4) 


where k, is the probability of r arrivals during 4 service time of arbitrary 


duration and is given by 


(5) i= if [e-2e(Av)r/r {Ib(v)dv 


The probability generating function (p.g.f.) of the distribution {k,} is given by 


(6) Ke) = © bear = fe b(v)dv 
r=0 0 
= had), tie <0 
where 
(7) B(s) = f e-svb(u)dv = L{b(0)] 
0 


of the service time probability density func- 
the Laplace-Stieltjes transform (L.S.T.) of 
function (P.D.F.) Biv): 


is the Laplace transform (L.T.) 
tion (p.d.f.) b(v) or, equivalently, 
the corresponding probability distribution 


(8) 


(15) (ni 
mm — 7_in-l)p — 7) Pn 
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Clearly, from (6) we have 
(9) 
K"(1) = —AB(0) 


vh sre th i 5) stead | 
rr ee pr imes len O iff 
v € € aS ( fo] differenti: Z 
Ove srentiati yn A 3 
t t ( Iso from (7) 


(10) B’ f 
—6'(0) = J vb(v)\dv = I/u, say 


0 
which i 
hich is the mean service time 
(11) 


Thus, (9) and (10) yield 
K'(1) = Alu = 

. ‘ . ; ie 3 
which is the traffic intensity ‘ 


From (3) and ( 


4) it follows that iti 
a 8 that the transition probability 


P = (py) 


matrix 


is given by 


mana ky ky ks 


(13) Ban cae ae 


Ch e { Live ¢ e-ste I 5 { ( I } b es e 
1e € le ments 
m of 
P 21 one-s *p t ansition wropna ] 
: ; % } « 1 iti 8 1n th 


distribution of NV 


fede ie. sense that if 
) he probability vector giving the probability 
m, the number in the system just af | as as 
ee ra just after the departure of any 
probability distribution of N i 

, on of Nm, will be given by 


14 
(14) Wm) — (0p 


Similarly, if a” j : 
arly, is the probability vec 
a y ctor for N m (n =]. 2 
in = ), we shall 


have 


7?) — 7p — 7_op2 


I} jO ! 


S ueipul 


Sor 


Mena 


spuepue} 


neo OnaNd uoluny” i] 


a 
ss ¢ 


f% 
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Thus, it is clear that P” is the matrix giving the n-step transition probabilities, 

whose (i, j)-th element, Say ps, is 

eS) aig = P{Nm ae jl Nm = i} 

which is independent of m as discussed earlier. 

for example, Prabhu (1965a), p: 65, problem 25, 
- 1, the Markov chain {N,,} 


shown (see; 
that when Pp <— 


on Markov chains) 
a steady-state vector 
nm = (Mo; M1> 72 a) 


Tt can be 
or a text book 


is ergodic, te. 
(17) = 
exists, where 
(18) T= in 
Then, using (14), we S°¢ that the steady-state probability vector 7 is given by 


(19) 


From (19), we have 


n= 1P 


j+1 
my = Tokyr 2 mij t (j 2 1) 
= 


(20) 
Let us define the p.g-f. 


(21) II(2) == mgd (121 <) 


Multiplying (20) by 2 and summing over j, We have 


ao co §f+1 
= Mo = kyzi+ a: = mgky_ass2) 


II (z) 


(22) 
=0 j 


wo a 
= 1,K(z)+ hoes miky_tss 
ial j=i-l 


changing the order of summation in the second term. 


the second term 4s 


(23) Eng ( 3S byte 


= {K@)/A{1@)—ml 


Thus, (22) simplies to 


(34) I(z) = mK@)+ {K(z) {0 @)—70] 


Now, we can write 
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which gives 


(25) 
I(z) = To(1—z)K (z) 
K(z)—z 


J fi d > e ta { that lim Iz = ] B 7 1 He it a] e. we he ve 
’ =] e t , 
Lo n 77 we us th f UC Zz 0 
0 . Dy yp ] rul 1a 


(26) 
oe 
lim T(z) = m)K(1—0) lim —'—)) 


zn i 253 


Now, K(1—0) = 
( 0) 1 and from (11), K’(1) = 
os )= p; therefore 


and, hence, ail 
(28) 
Ie) = L=P)\—2) Ke) 
Ke)—2 


(28) gives th 
s the generating functi 
so that 7; i g function II(z) of steady-s 
4 18 the coefficie ; (z) of steady-state probabiliti : 
of (28). lent of 2/ in the expansion of the A ahrye {ny}, 
of the right-hand member 


Waitin gq Tim v2 


Let N deno 
4 te the steady 
departing ct ady-state number ir : 
sustome Gk n the system ie 

the ee : al and w his total time spent in tl his left behind by a 

of arrivals during time w. Thus 1e system, Then N is 
(29) } 

ty = P{N = ae . Len re 
IS J [e*w(Aw)4/j Y] f(w)dw 


where f(w) i 

s the probability F 
be ae P ility density function of w. ¥F i 
umming over j, we haye w. From (29), on multiplying 


(30) a 
II(z) = Al 
(2) J e~A-zhf(w)dw = f(A—Az) 
where 
(31) ~ 
Fis) = J ef (w)dw 


is the L.T. of f(w) 
eae WwW), Now, let fF (- b 
spent in queue , let fo(:) be the p.df. of wg, the queueing time (ti 
be hi ia ey of the typical departing custome eam queueing time (time 
is service time; then r in steady state and let v 


(32) 
wo = Wgtv 
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and hence 
(33) {w= Sav) « b(w) 


i convoluti si 2 ¢ ion theorem of 
where the asterisk denotes conv olution. Using the convolut 


Laplace transforms, we have 


Fis) = Fols)b(s) 


and b(w), respectively. From (30) 


(34) 


where f,(s) and b(s) are the L.T.’s of fa(w) 
and (34) it follows that 


T(z) = fg(A—Az)b(A—A2) 


(35) 
or, using (6), 
(36) (2) = fo(A—A2)K@) 
Now, writing A—Az = s in (36), we have 
7 ze) — J (l—s/A) 
en fal) = “Rafa 
whence, on using (28) and (6), 
i ee 
(38) Sa(s) ai s—A+AB(s) 


From (28) and (38) we can find the mean number in the system and the 


mean queueing time. The mean number in the system is 


a ek, MEH (v?) _ : 

(39) E(N) = 1'(1) = p+ 3q— py = DL, say 
and the mean queueing time is 

NR re a Deo : 

(40) E(wq) = —J,(0) = ln). Wa, say 


where E(v?) is the expectation of v? or the second moment, about the origin, 
of the service time distribution. 


If W be the mean time spent in the system, we have 


(41) W = W,+(1/z4) 
From (39)-(41), it follows that the two relations 
(42) L=AW, La=AwW, 


hold for this system also. 
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In (39) it may be noted that 
(43) LE 
U(v?) — oF+1/p2 
where o? is the vari j 
2 ariance of the service time y Hence (39) k 
p e (99) may be written 


(44) 


L= p+ Mot pr 
mp} 
and, therefore, 
“ie Yi ks ee 
2(1—p) 


Deductions for M/M)\ 


Considering the special case of M/M/1, we ha 
AVL} 2 ., WE Nave 


(46) 
b(v) = jee-#” 


and, therefore, 


(47) i pe 


(48) = 67 
K(@) = b(A—Az) = p[(y+A—2z) = (l+-p—pz)} 


Hence (28) reduces to 


(49) ej 2 
1—pz 
and (38) simplifies to 
(50) Jae) = Seer) 
s+yu—A 


(49) gives the well-known M/M/ 


1 system-size distributi } 7 
- e distribution [Cf. 1.5(17)] : 


my = (1—p)pi 
Also, (50) can be written 
52 
a) Fale) = (1—p)+-(1—p). —A_ 
s+ p—A 

whence, on taking the inverse Laplace transform, we get 
(53 

) fol) = (1—p)8(w)-++(1—p)Ae@—mne 


which i 5 7 
ch agrees with 1,5(26), (27). [d(w) is the Dirac delta function.] 
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pone 


Polity 


m 1| enjppuels Ue!PUl ht 


,o |91and uoluf) ; 
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Higher moments es < 
V - in the system, 8 ven by 
] e of T, the number in shi y 
i T > arL_ance AV 
(1) The 2 g 


o WY 7 2 
Var[N] = BiN?|—(BW)) 


E | 8 i >] oted that 
| L is given by (39). Als Oo, it may be no 
there E# NI =. y 
where | j 
Il "(1—0) = iim 7a m) 


= BLN(N—1)] 
— B{N?)]—2IN] 


imi curring in 
ate the limit occurri 
f N, of order 2. To evaluate the 
i roatorial moment Of JV, 
is the factorlé a 
(55), we write (28) in the form 
55), 


Ne{K@—4 = (1—p)(1—2) K(@) 


e ee sniligving, | We 
io 6) w.r.t. 2 thrice, proceeding to limits and simplifying 
Differentiating (56) w.r.. @ 

oh A2B(v?) , AMH(v*) aba wil) 

(57) I*(1—0) = i-p | 3(1—p) 2(1—p) 


Using (54), (55), we have 

2 
(58) Var[NV] = ’(1—0)4+- L—L 
5 
(1—0) are given by (39) and (57). 


” 
where L and Il o 
1 with constant service time 1/p, 


Considering the particular case of M/D| 


e have 
: B(v?) = 1/p?, B®) = 1/3 
(59) 
so that ‘ ‘ : i 


p*(p?—4p+ 8) 
~ 6p)? 


Again, for M/M/I, since 
Blo?) = 2/2, Blo®) = 6|p? 


(61) 


Sec. 3.8 | 


NON-POISSON QUEUKS 49 
we have 
(62) (a Oj ee “Ae ae 
i p i pP (1 py 
2p* 
(1—p)? 


which tallies with 1,5 (34). 
(ti) The variance of wq, the waiting time in the queue, is given by 


(63) Varlwe] = f;(0)—(f, (0))2 


where f(s) is the Laplace transform of the probability density function of Wa 
given by (38) and - ~f,(0) = Ewq| = Wg, given by (40). 


To evaluate fi(0), we write (38) as 
(64) {8—A-+AB(s) }fa(s) = (1—p)s 


and differentiate w.r.t. 8 thrice and proceed to limits as s — 0. Thus, on 
simplification, we have 


Feo) = MB Wal(v*) + B(v4)} 


65 
” 3(1—p) 


LqE(v*) +-AEB(v)/3 
| —O 


Ex. 1. [f Q(x) = N(x)/D(x) be indeterminate of the form 0/0 at « =a 
t.e., N(a) = D(a) = 0, show that 


; 


a D'N’—N'p" 
: i ‘ ” — j é 
(1) i a ()' (x) jim si 


where the primes (dashes) stand for differentiation w.r.t. x, and the arguments 
have been suppressed on the right-hand side for simplicity. 
We have 
DN’—ND' 
0’ a= 
d (x) Di 


which is also of the form 0/0 at x = a. Applying the Hopital’s rule twice 
we have the result. 
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oe 7 ; 1 
is very important and often ust d in qu oe 
rth, mean waiting time, etc. For examp/e, in 
1 using Eq. (1) above, 


Eq. (1) above 
for finding mean queue leng 
z)K(z), D= K(z)—2 anc 


Eq. 3.3(28), taking N = (l— 
we have 
USO 


Now, K’(1) =p, and fiom 916); ee) — A26"(0) = A2H(v2). Hence 3.3(39) 


follows. 
and similar results in 


The same method can be used to obtain 3.3(40), 


the following sections. 


3.4 The System GI/M/1 


r the single-server qu 
and exponential service time with mean 1/p. 


he n-th customer (” = 1, 2,...), and let N(é) 
Further, we write 


eueing system with general inter- 


Let us conside 
arrival time distribution A(u) 
Let #,, be the arrival time of t 
denote the number in the system at time t. 


(1) N,, = Ntr—9) 


so that VV, 
since the service process is 
chain imbedded in the process Nii). 
uring an inter arrival time oO 


is the number in the system just before the n-th arrival. Then, 
{N,,} constitutes a Markov 


Markovian, the set 
Let k, be the probability of r service 


completions d f arbitrary duration, given by 

(2) ly = f [eminent |r  aA(a) 
0 

Then, the transition probabilities of the Markov chain are 


(3) Mm =han CP 0<76 tT) 


(4) 


ll 


Pad Spy (ult 
be qw [1- Be r ! 


i 
= 1l— & kh = %, say 
r=0 


It may be noted that 


(5) 


m= Kinyt hye t tee 
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3) 


That Dio = si Vv | 8 
7 = a 
0 % May also be arrived at as follow 


| } - s * . 
: i Ls ju e 15) ne © ary 9ecome 
in t 16 8 stem st ft ok t 16 ival b 2COMeSs 4 
J SU e ¥] Cxt arr if he ) Vv v 
ust bef re bl 16 ne t 1V al t 2 Service f 


, 
“4 The number of customers 
+1 and this will drop to zero 


least . woe 
2+-1 service : acility is cans : 

Fl services during Qn iitkawcatees y 1s capable of cor a 
g inter-arrival time mpleting at 


The transiti 
ansition probability matri 
atrix P = (py) is gi 
4) is given by 


( 
Gaia! a 8 


0 
Hy ky 0 0 


1 


(8) a ae le i oa: 
P — 2 Xe k, ky le 
3 


1 pec 3.3 t e steady-state I | bilities are given 
As Ir Sex tion 3 
’ h CAC Sté > probabilit a ¢ 
4 & by 


(7) : 
7 aes Thy 
i=0 
(8) eh 
ae ee Wki_jus (j > 1) 


As a trial solution of (7), (8) we take 
(9) 

mj = abi (a>O0, 0<b<1) 
Substitution of (9) in (7) gives 


(10) a= Sabi( ¥ t) 


Lies | 


ie] 
a ak,bi ° 
eiikior 3 changing the order of summation 
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or 
, We have NO 
? N-POISSO’ Ee 
(11) K(b) =6 i ISSON QUEUES 53 
id fi+-0) =a'¢ 
whert Also, )=00, f(1) , 
oo 
(12) K() = z kat (jz| <2) (20) 
s, "(2 6 « 
oe I) = Degz-9-4 Dhey+ Cheye-t 5 
opatiks PLA oe 
ore, f’(z) is a monotonic incr 
re 


(9) in ( ) Vv 
2BIN fu Hi é 
£ th n und i h me { a] 


And, substitution of 
<= | if and or re 
il / 
yYuUsl)=e—1 > 0, since we have f’(0 
a ? ) HD. 


is the generating function of Kr. 
indO<z 
oo F F 
(13) abi= & abiky_gu1 
i=j-1 
f(z) 


oo 

— abi X k,bt = abJ-1K (b) 
r=0 

which also implies that K(b) = 56. Thus we see that (9) is a solution of (7) 


8) if b is a root of the equation 


and ( 

(14) K(z) =2 

jc. b is a root, less than 1, of the equation 
(15) O(p— pz) = % 


T, of inter-arrival time 


and (12), where a(s) is the L. 


as can be seen from (2) 
a from the normalizing condition 


density function a(v). Evaluating 
oo Let {'@) = ( 
: eae J(@)=0atz= - 
Lm (0,2) and me BPs (0 << 4% < 1); then f() is n 
motonic increasing in (z 1] oN nonotonic decreasing in 
© “1) . since fi) =) = 


(16) 
j-0 
Se fee a Bay’ ' a 
1 und hence there exists exactly one val it follows that 
ven > value Zz (0 < ; 


f@q)=1, or EK 
. , A (2, — 4 T r4 yA 1 
(@) =%. We can a fortiori assert that 0 <%) such that 
that O<2<1. This 


trie distribution 
completes tl : 
a Hats) mi 
proof of the theorem. =o 


(17) my = (1—b)bs 


we obtain the geome 


where b is as given above. Thus 
. 8, Wwe see that (17) : 
at . determine 
uniquely nes the steady-st 
. y-state probabilitic 
€ nt 714 


It remains to be shown that the equation K(z) = 2 has a root b such that 


Lich. we do in the followin : 
: 8 Queueing Time. 


Using (17 

: est 7) and proceedi Heke 

obtain the distribution of a proceeding as in Section 1.5, w 

of the queueing time w, of 5, we can 
) Wq Of an arriving 


and assume that hk, > 9, steady state as follows 


gO b = 1, wh 
customer in the 


Theorem. Let Kz) = 5 Iya be a pee, 
rea 0 


o< K(ljp=ae< Then the equation K(z) = 2 has a root z such that: (21) 
0<%< lif and only if a > 1. [See Feller (1968); Prabhu (1965a), p. 63.] P(wg = 0) = 2, = 1—b 
(22) Pott 
Pw < wy < w+dw) = z (1—b)bd- (yw)i-le-#e 
j=l Gj—1)1 + pedw 


Let us consider the function 
= (1—b) bye-bye dw 


(18) fe) = K@e = Mole +t... (0<% <4) 
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(‘17 d (21), ( 
The relations (17) anc ( 


ae 
; + nortant 
54 AN INTR 92) are imp 


tric and the que 


the interarrival time 
form of 
of the 4 


f M/M/1 queue, with, of course, 
ne he unique root in0<2z2<1lof 


. 6 ) Ss eo e 
1 5 

Z jistri 1101 

en A 

he sy stem 1 i oO 1e 


spective 
exponen 


bial irre 

is truncated : 

distribution, just = 

different parameters, ee 

the equation a(pe— HA ooh 
For the M/M/1 system, 


i(s) = AJ(A+8) 


in the 
od by b, t 


hey were 


being repla 


(23) 
Therefore, 3 


(24) 


is the root of 


AyAt+H—HA) = # 


or, 


: — ¢ : s b e. 
. : and p: H n b > 
in s are be 


3.5 GI/M/c ‘ ‘3 applic 
r thod of Section 3.4 18 apr "et 
nee GI/M/ec. As before, w 

system Gj"; tive arrivals as 


ble to the more general case of the 
~ % . . 
hall have the one-step transition 
multichannel sc 

abilities between two cons 


prob cs 0, 0 <j < i+1) 


y — k 4 x . 
(1) ae are completed during an inter- 


. » gervices 
ability that 7 SeTv- 10 as all the ¢ servers are busy, 
* 1 fo 


zis b 
where K, is the pro Now, as lo 


arTrlva time 8) AT tI ry ] ration 
] J f a *b a Cc . 
a a 


'>c—1, we have 
which requires that 1 2 ¢ 1, 


@ eee (Huy GA(u) 


(2) kp = 7 | 


HE ge 
The steady-state probabilities 7 are given y 
ne stead) 
co 
ae x MP 
(3) 2, 


| x ve olve where e Vv by sa trial 
e bt i 7 (1 and (2). A e 

I ha to al iG (3), h Te the Pij are gi ren. A ( ) 

1us, W' , ' AY)» 


8 on (sug; i GI/N s write 
lution (suggested by the solution for /M/1), let u 

8 lon (sugs , . : 
‘ : Tj = OW, 0<0<—1,3 > c—l, C a constant 
( ) ? 


Substitution from (4) in (3) yields 


co 
S Ai-jH = + Ok 
(5) 6= = 6 Kegs ra . 


imj—1 
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oO 
ao 


Now, using 


(6) 


(2) and simplifying, we have 


O = G(cu—cpl)) 
where, as before, G(s) 
in Section 3.4, it ¢ 
let us call it b, 


is the Laplace- 
an be shown that ( 
Then, for 9 = b, 


Thus, we have 


Stieltjes transform of A(u). Now, ag 


que reot in 0 < 9 = i; 
the values of 7 (j > c—1) are given by (4). 


6) possess a uni 


(7) m™ = Cbi(j> c—1) 
The evaluation of OC and th 


e determina 
complicated 


and will be discussed | 
below some interesting results which 


tion of 7 for j < c—1 are rather 
ater only for the case ¢ — 2. 


We give 
can be obtained without e 


valuation of C. 
The conditional probability 


that there are j units in the system, 
that the customers have to que 


given 
ue for Service, is 


eo 
(8) m™| % mm = (1—b)bt-, (5 > c) 
r=c 
Therefore, the conditional 


queueing time distribution of the customers who 
have to queue up has a de 


nsity function fa (w) given by 


(9) falvo)dw = ¥ (1—byps-c, (oprw)ee-emw endw 
j=c (j—c) ! 
= (Il—b)cu exp[—(1—b)cuw)dw 
(8) and (9) prove the interesting results 
of the number in the system for customer, 
queueing time are, respectively, the geom 
tions for the GI/M/ec system as they were 


that the conditional distributions 
8 who have to queue and of their 
etric and the exponential distribu- 
for M/M/1. 

Note that in writing equation (9) we have used the fact that if there 
j units in the system, i.e. j—e in the queue, an airiving customer has 


to wait 
in queue until j—c+1 services are completed. 


are 


3.6 The Supplementary Variable Technique 


Besides the imbedded Markoy chain method 


, another way to analyze 
the non-Poisson queues 


is to use one or more extra variables, called supple- 
mentary variables in the definition of states of the system, so that the process 
with the new definition be For the M/G/1 or GI/M/1 


ary variable, being the time since the 
ectively. 


comes Markovian. 
system we need only one supplement 


last departure or last arrival, resp 


The most general queueing 
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y using c+] supplementary variah 
rrival and c for the elapsed service tim ee 
the method below for the system ee 
1 has the advantage Over the imine r 
cide 


y the system in conti 
) INUOUS tj 
Ime 


56 

¢ can thus be studied by 

ince the last @ 

We illustrate 
methoc 

nere we cad stud 


stem GL/G 
e for the time 5 
¢ servers: 
variable 


sy 

on 

each of the 

The supplementary 

Markov chain method that | 

instead of at discrete time points. 

£/G/1 

joint probability that at time ¢ ther 

here are 


0), be the 
and the elapsed service time of the unit in seryj 
ry. 
let Po(t) be the probability that ia 
6 


-state, We consider the limit 
ing 


The system I 

Let P,(*; t)dx, (” 2 

rits waiting in the queue 
ad). Also, 


terval (#, 
at time ¢. For the steady 


n ut 
lies in the in 
system is empty 
probability density 


P,,(#) = lim P,,(a, #), (” > 0) 


» i> @ 
ting probability 
Po = lim Pol) 


(2) ae 
n distributed with mean rate A. Also, let n(x) b 

: ; er e 

a service completion, given that the 


and the limi 


be Poisso 
obability density of 
so that 


Let the arrivals 
the conditional pr 
elapsed service time is *, 
: b(x 

(3) (x) = oe 
are, respectively, the distribution function and the density 


where B(v) and b(v) 
From (3), we have, on integration, 


function of the service time. 


(4) Biv) = le. le oi 
and, therefore, 
(5) b(v) = atee 6 i 
Arguing as in Section 1.5, we have the following transition equations. 
(6) P,,(a+Az) = P (x)[1—(A+q(a)}At]+ Pal) AAe, (n > 1) 
(7) P,(v+Az) = P,(x)[1—fA+7(@)}A7] 
8 ao 
(8) Pop = Pop(t—AA*) + f n(x)Ac. Py(x)du 
From (6)-(8), we have : 
d 
9 Leas od 
(9) 4 Pla) = [A+n(e)Pale}+APaate) (m2 Y) 


EFI os — 


aye? 


Sec. 3.6] 
br NON- 
F POISSON QUEUES 57 
ns Ae) = — 
e rs [A+-9(x)|P,(z) 
| AP oy = f P 
1ese equations are to b ; Be 
' 0 be solved wi 
“s solved with the followin I 
nO) = ~ £ boundary conditi 
2 I i Pvila)y(ax)der, (n > 1) oye 
P,(0) = f P 
and the normalizi : sah a 
rmalizing conditic : 
(14) a sy P ge 
| ek «a @ 
head the probabili ry 6 Ses 
of the valu a , 
the value of x, is gi c a — en 
Gs ies sons in the queue, ir 
"ue, irrespective 
P,= {Pid | 
P,= P41 f 
efine the Lae 
Define th generating functi a 7 
tio 
Z ictions 
Paz, x) = » 
a, 7%) = X P(x 
¥ -, n(X)z", jz] <1 
, » 
i. P j Paz, a)dx 
Multiplying (9) by 2", s a ce 
, Summing over 
‘ n, and adding (10), 5 
g (10), we have 
aE a(@, ©)+[A(1—z)+ P. == 
whose solution is a oan: 
(21) “= P| —AUL—z)a— 
Paz 
q\*, 4 
) a(z, 0) exp [—A(1—z)x fa )d 
{oe C v)dv 
ot = F(z, 0) exp [—A(1— yarl{ £ : 
F / z)x},1—B(x)} 
(22) = CD eS) 
Se = 
Paz) Pz, 0) - 1~8(A—2z) 
whe v 
re we have used the fact that Sal 
1a 
(23) 


J e-8* Biw)da = € b(v)d dx = 6b(s)/s 


0 
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papcae! ne a & 
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58 Al from (12 4 
Plt: 0), we have: 
yaluate dding, 
poe powers of z and ac % endu-+-AP 
priate Pp l t s P,,()9(@)2 axe - 00 
Pg(z, 9) = Zz 0 n=! 
(24) a ; 
1 ) + 1 { Pye; x)n(x)da 
= AP oo \ l Z ee 
ee . 21) in 
Y P,(z,2) from (2 (24), 
ert substituting Nee : BE by 
- using (11). Now. ee w that Py lz, 0) isg NM 


x or al, 
evaluating the integ AP oo (1 —?) 
Pal, 9) = 5(A—Az)—2 

(25) 


7 (95). we have 

From (22) and (25), W Poo [1—2( b(A—Az)] 
- Pye) = — ¢a—Aa)—2 
(26 


1e (19) gives 
Hence (19) Pelt 


Pe) = Fa—A2)— 


(27) 
Hépital’s rule 
7 sing ’Hépital’s rule, 
)) = 1, we have, using 
Since P(1—9) ; 
Poo = 1—-P 
== L=AB'(0). Therefore, 


(28) 


where p is the traffic intensity 
= P@) = 5a—da)—2 


lation between the random point ihe, Lites. 
M 7 i robabilities 

‘ties P. studied above and the imbedded Markov chain pro ae 

ee Let P+ be the steady-state probability that there ar 

‘ (11), we have 


We proceed to find the re 


Section 3.3. 
' D 
the system just after a departure; then, as i 


(30) APY = f P, (e)nlade, (mB 0) 


Using a probability generating function 
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AP *(z) =f ¥ 2™P,,(x)n(x\dax 


0 ned 


@ 


= J Pal, x)n(x)dx 
0 
= Pelz, 0)b(A—Az) 
on using (21) and integrating. Now, using (25) and (28) in (31), we have 
(32) P+(z) = (1—p)(1—z le 
b(A—Ag)— 
It may be noted that P+ is 


; the same as the 
Eq. (32) above is the s 


7, Of Section 3.3 and that 
same as Kq. 3.3(28 


Cle early, the generating function of P 
vl by (2) +P oo = *(z). Hence, 


(33) 


n, the random epoch probabilities, js 


P+ — Pan = Oe 
From (29), by differentiating P(z) 


wr.t.z and letting z4 10. we have 


Lg = P(i—o) — A%8"(0) _ A*E(v*) 
2(1—p) 2(1—p) 

which is the mean number in the queue at a random epoch, If we write L+ 

a 

for the mean number in the queue at de *parture inst 


ants, we can similarly see 
from (32), or 


as already proved in 3.3(45) that 


2h p) 2(1—p) 


(34) Be i A2B (vy?) p*4 A2e? 


Particular cases 


(a) M/D/1. Here o, — 0 and henee 


(35) Po ee 
“~~ 2(1=p) 
which is half the mean queue length for M/M/1. 


(b) M/E,/1. Here we have, using 1.4(18), 


(36) Eo) = f v®{(ku)tok-e-¥n0/(k1) 1) dy = = 
0 fe 
Therefore, 
_ ae a? k+1 p> k+1- A? 
Pe ae a i jiu) 


It can be seen that (35) can be deduced from (37) by letting k + o, (see 
1.4(19)). Also, setting k= 1 in (37), we obtain DL, for M/M/1 


9c o1GN, UOIUNF 
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, 
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HU 


H i opba € we ean ob ain L I} 
y o 1 DL ae, 5(21) ( 3)) : 
st que i€ ing systems (sec oer}, \=* 

mo : 


which hold for We ‘ Lala 
(35) 4 i 
(39) fs ; : : 
o L= Lg r Pp 
( { 


_ : ‘ | 
a. Poisson arrivals at a mean 


=x — 7. with 

| et us conside r the e-channe ] SV stem is . : 

a d hs ant serv ice time equal to 0. Le t ] n t be th pl | yabilit , 

rate / and cons ti 8 . >, fe re have 

he t the re are 7 uni S in the sys te m a time ; he nh we e 
tha are t > : > t t t } 


P.(t--b) = S Pilt)e*[(Ab)"/n 1) 


(1) " e 
Ss P,(t)e*0[ (Ab)! (n—i-+c) !] 
j=c+1 


rotate . ahili ies 
Let us assume that the steady-state probabilitie 


Pe ces lim: 2 fe) 


(2) ae 


exist; then from (1) we have 

nee 
as (P Py. Pejky, D> Piky_ire 
ct 0! 1 (Ban 


and simplifies to the Poisson probability 


n 


(3) 2 ad 


where k, is as defined in 3.3(5) 
function with mean Ab. 


Using the generating functions 


(4) i ae n=O 
. y »—AB(1—z) 
(5) K(z) = 2 Fn h=€ 


we have, from (3), 


xX Py(z'—z°) 
i=0 2 


1—2/K(2) 


(6) 


which can be written in the form 


c—1 ‘ 
(l—z) & Qi! 


i=0 
(7) P(z) — 1 —2%eAb(l—2) 
where 


Qr = Po +P,+...+Pr (r > 0) 


nage SYREN 


Cy 


z= 
: ae 
Sec. 8.7 | 
NON-POISSON QUEUES 6 
Now, we have 
¢=j 
I lim P(e) = 0 
'—>1-—0 O—Ab 
80 that 
(10) s 
p> Or = © -Ab 
teoO 
: 
Thus we should have 
11) 
( P 7 Able mc) 
Which is therefore, t} iti ‘ 
‘ ; 16 6 ; . s eX ; 
© condition for the existence of the steady state, 
Now, for iciently s iti 
sufficiently small positive 3, we have, at jz] 1-+-6, 
(12 >—AB(1—2) -AD+AD (2 
) |é | < enAd+Ad/(2) — enh = |] +d)¢ = |zJe 
since Ab < ¢. fr Tene My 
A . ¢, from (1 1). Hence, by Rouché’s theorem, the denominator of (7) 
1a8 exactly ¢ zeros inside the cin. > |z i 
: 8 inside the cirel |2| 1-8, which we denote by 2, Z» 
®e1,2e = 1, Since the orator of (7 i de 
*e—1) 2 Since the humerator of (7) must also vanish at these values of z, 
it follows that the equation 
: e—1 
(13) x Qiz' = 0 
t=0 
8 ae Ee Pe 
has the c—1 roots “Is “a +++) 2.4. It can be shown that % (r= 1, ...,¢—1) 
are a isti 3 3 re nave ec i 
: e all distinct, Thus, we have c—1 equations and (10) to determine the ¢ un- 
cnowns Q; (i = 3 30 t i 
= Gs (§= 0,1, ... 6 1). So the P(z), given by (7) or (6), is completely 
determined, 
As a particular case, for the system M/D/1, we have on setting ¢ = | 
in (6), 
> a 
(14) P(e) = —*0o(l—z) 
| —zeAb(1-z) 
where 
K 
(15) Po = islbe l—p 
Thus, 
BI ‘> coe ~ F D ‘a> 
(16) P(e) = (1—p)(1—2)/[1—-zer-2)] 
Fgh 
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It may be noted that (16) agrees with 3.6 (32), as could be i d an 
(16) the mean queue length and waiting time can be found as befor et 


the mean number in the system 1s 


Spare | Pp 
(17) L=P(i—0) = Sip)” p+ 2(1—p) 
Be 
(18) lg Oe Sp 
a9) Wa LX 
(20) Wa = Lyla 


3.8 The System M/E,/1 
The single-server queue with Poisson arrivals and k-Erlang service time 
distribution admits of simple analysis as follows. Let A and y« be the mean 
arrival and service rates, respectively. Then (Cf. Section 1.4(c)), each service 
time may be considered as made up of & phases, where the time spent in each 
phase is distributed according to an identical negative exponential distribu- 
tion with mean 1/kyu, so that the mean service time, being the time spent in 
all the & phases is, 1/z. The phases are numbered in the reverse order, so 
that the service starts with phase &, after finishing which the phase k—1 
begins, and so on until the service ends at the completion of phase 1. 
[Equivalently the phases could be numbered in direct order : 1, 2, ..., k.] Some 
authors further generalize the model by introducing a variable number of 
phases, so that an arriving customer may require r phases with probability C, 
(r = 1,2, ...). In such a case, numbering the phases in reverse order is parti- 
cularly suitable, so that the service starts at phase r and stops on completion 
of phase 1. 


Let Py;r(t) be the probability that at time ¢ there are n units in the queue 
and the unit being served is in phase r of service, and let P(t) be the corres- 
ponding probability that the system is empty. Then the transition equations 
are (ignoring the o(At) terms) 


(1) Par(t+At) = Par(t){1—(A-+hp)AR+P,,_,, r(t)-AAt 
+P raat) uA, (n> 1; 1<¢7<k) 
(2) Pra(t+-At) = Pra(t){1—(A+ku)AR+P,,_,, e(t)-AAt 
+P (t)kwAt, (n > 1) 
(3) Pop(t-+At) = Pop(t){1—(A+-ku)Ab-+ Py, 74 (t)-kuAt, (1 <r < k) 
(4) Po,(t-+-At) = Porlt){(1—(A+ kn)At} + Pog(t)-AAt+-P,,(t)- hut 
(5) Poo(t-+-At) = Poo(t)(1—AAt) + Py, (t)-kuAt 


a ee 
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From the above equations 


we have the following 
babilities P,,,, ete., ; 


(6) O= —(A+kp)Py FAP, rkpP,, 41 
(7) O= —(A+kw)P ap AP,-1, k4 kpeP nit, 1 
(8) O = —(A+kp)Poet kpPy, v4 

(9) 0= —(A+ky)P oe AP oo kepP,, 

(10) 0 = —AP y+ kpP oy 


Let us define the generating functions 


(11) F(z) = EPs, je] 1 
n=0 


(12) k 
n=0 
Multiplying eq. (6)-(9) 
(11), (12), we have 


(13) 


TAY Poot (kuy*/a)[F (2) —Po] 


or, simplifying and using (10), 


(14) (A+ ku—Au—kply)G(a, y) = kyly*/2—1)F (er) 4 Ay*(1—1/a)E 


To evaluate F(x), we put 


ges Ap ee 
A+ku—Ax Yar SAY; 
in (14), which gives 
A(l—2x)P A(L—«)P 
(15) F. (x) = " 00_ = =a. 00 — 
; ku(1—a/y*) kul j—2( Ath Ax \*) 
L \ kyu Pay 


Substituting for F',(x) from (15) in (14), setting y = 1 and simplifying, we have 


[( A+kp—As \# 

Pro| ( “tina pax 

(16) G(x, 1) = ri Wm Ta. oY 
i) 


for the steady-state pro- 


k ie) C 
Ole, y) = X Pray = EE Pryaryr, ty) <1 
r=l Ci : Sia 
by appropriate powers of « and y, adding 


0= —(A+kn)G(a, Y)+ArcG (a, y)- (kuly)(Gla, y¥)—yF (x)] 


OP 


<i RC ih ake) 


moe ueipUu JO ft 


and using 


| 
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Clearly, 
os 
(17) G(a,1)= = P,xm, |x| <1 
i n=0 


i i ion of ‘obabilities f finding n in the queue 
is the generating function of the probabilities P,, of finding n 1 , 


service ase. Fr j aking limits and using 
irrespective of the service phase. From (16), taking limit g 


l’Hopital’s rule, we have 


we 
lim G(a, 1) = Poop 


(18) cpl 1—p 

where p = A/u. Also, from (17), we have 

(19) lim G(x, 1) = 1—P,, 
ctl 


From (18) and (19) we have, as usual, 
(20) Py = 1—p 


Substituting for P,, in (16), we have the generating function of the queue 
length probabilities. From (16), differentiating w.r.t. x and taking limits as 


x — 1—0, we have the mean number in the queue 
(21) b OMmee tt 66 (aed 0 | Jem) ee 


Using eq. 3.6(38)-(40), we can find Wa, W and L. 


It may be noted that (16) agrees with 3.6(26), since in the present model 
the Laplace transform of the service time density is 


(22) b(s) = = 


Special Cases. Taking k = 1 in (21), we have the well-known result for 
M/M/1 


2 
(23) Ig = 
as es 

And, letting t + 00, we have for M/D/1 


(24) Ts =5 4 


which agrees with 3.7(18). 
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Direct Derivation. We gi 3 : 
dttorste give below a method of obtaining Wa, W and L 


We know that P,, is 


the ili ee 
and, therefore, n+1 Probability that there are ” units in the queue 


in the system (n = 0, Li 2, a) ences 


(25) L= DH (n+1)P, 
0 


n= 


I 


2 nP,,-- z P. 
r= 0 


I 


Iq+(1—P,,) = L,+p, using (20) 


ie k+-1 p? ; 
eee Mar "Top? Using (21) 


Also, Pry is the probability that there are (nk-+-r) phases of service in 
the system (yet to be completed). Hence, the expected number of phases 
in the system, say Ly, is given by 


o ok 
(29) Ip= XX (nk-+-r)P.,, 


nmr] 


=kKXUEZ aP,,+ DU rP,, 
EN nr 


i 2A 
=k ae G(x, 1) bar se G(1, y) 


aa y=1 


From (21), it is clear that the first term on the right-hand side of (26) is kL,, 
For evaluating the second term, we have, taking « = 1 in (14) 


> 


(27) a1, 9) = ew) 

and, letting «1 in (15), we have 

(28) Fy(1) = Afku 

Therefore, we have 

(29) G1, y) = (yE+ yk... +-y)(Afku) 

Hence, 

(30) ouw| = MY Ae, 
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Thus, (26) simplifies to 


k+l pp? k+1 
(31) ly = ap, Sap a 
Be 


Since each phase occupies an average time 1/ky, the mean time 
K € eae 2 
queue is 


Spent in the 
(32) 


Wo = Ly/kp 
k+1 es Pot. a k+1 A 
~ Be  w(l—p) 2h M(u—A) 


gain, since the incoming customer adds & phases to the sys 

Again, since the ee us I ystem, the inbae 
time in the system is 

(33) W = (Lp+h)/ku 


= Wg+(1/z) 
3.9 Hyper-Exponential Service Times 


It may be noted that the variance of a k-Erlang distribution W 
1/u is equal to 1/ky? and is thus, for k > 1, less than the variance 1 


ith mean 
corresponding exponential distribution with the same mean. 


/ of the 
Thus an Brla; 
distribution is not useful to approximate a distribution w 
mean has a variance greater than that of the corres 


hich for a given 
distribution. For studying such distributions, consider 


ponding exponentia] 
: : two channels arranged 
in parallel (instead of the series-arrangement of Erlang distributions), Let 
the service times at the two channels be exponentially distributed with means 
1/(20p) and 1/{2(1—o)}, and let the customers be served at the two channels 
with probabilities o and I—g, respectively, (Q<¢o < 1). Without loss of 
generality we may assume o < 1/2. Then, the overall mean service time is 
(1) 
and the overall variance is 
(2 P 2°[or.2orpe-2ene 
) fe [o.2on€ + (1—0).2(1—o7 e202 ner dey (1 Jy2) 


een 


aes OL J 
~ Lemay] oe 

which has its least value 1/4? when ¢ = 1/2 

than 1/y? for all oth | ” 


and is, therefore, always greater 
er values of ¢ in (0, 1), Nay 


—— 
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Thus, the above distribution 
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suitable for approxim 


67 
» ealled hyper-exponential distribution, is 
ating distributions which have 

the exponential with the 


& greater variance than 
sonar: same mean, Similarly, we can 
exponential inter-arrival time distribution 


It is easy to see that for the s 
and hyper-exponential 


have a hyper- 


ingle-server queue with Poisson arrivals 
service 
we have, using 3.3 (45), 


times, (which may be denoted by M/HE/1), 

(3) Lg = ME(*)/[21—p)] = p*/[4o(1—o)(1—p)] 
which is greater than the value p?|(1—p) for the M/M/1 queue with the same 
parameters A and « expect when ¢ = 1/2. 

The values of L, W, Wa can be deduced from (3) by using Eq. 3.3(41), (42). 
Ex. 2. A hospital clinic has a doctor examining every patient brought 
in for a general check-up. The doctor averages 4 minutes on each phase of 
the check-up although the distribution of time spent on each phase is approxi- 
mately exponential. If each patient goes through four phases in the check- 
up and if the arrivals of the p 


atients to the doctor’s office are approximately 
Poisson at the average rate of three per hour, wh 


at is the average time spent 
by a patient waiting in the doctor's office 2 What is the average time spent 
in the examination? What is the most probable time spent in the 
examination ? [By permission from Sasieni. Yaspan and Friedman (1959)] 


This is an M/E,/1 system with A = 3 per hour, and 1/u = 


total service 
time = 4X4 minutes = 16 min. = 4/15 hour, so that = 


Hence, using 3.8(32), we have 


15/4 per hour. 


average time spent waiting in doctor’s office 


rene Fs 3 5 16 fe se 
= W,= id (s/Daa) 5 x cc hour = 2/3 hr. = 40 min. 


Average time spent in examination (service) = 1/ = 16 min. 
Most probable time for examination 


= modal value = (1—1/k)(1/u) = (3/4)(16) = 12 min. 
(See Problem Set 1, Prob. 7.) 


Ex. 3. At a certain airport it takes exactly 5 minutes to land an 
aeroplane, once it is given the signal to land. 


Although incoming planes 
have scheduled arrival times, the wide variability in arrival times produces 


an effect which makes the incoming planes appear to arrive in a Poisson 
fashion at an average rate of 6 per hour. This produces occasional stack-ups 
at the airport which can be dangerous and costly. 


Under these circumstances, 
how much time will a pilot expect to spend circling the field waiting to land ? 


[By permission from Sasieni, Yaspan and Friedman (1959)] 


age onand uociuf) 
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Here we have an M/D/1 queue with A = 6 per hour, and w = 12 per 
hour, so that p= A/u = 0.5. Therefore, 
Lq = p*[{2(1—p)} = 1/4. 
Hence, Wg = Lg/A = (1/24) hour = 2-5 minutes. 


Ex. 4. In a certain bank, customers arrive according to a Poisson 


distribution with a mean of 10 per hour. From observations on the teller’s 


performance, the mean service time is estimated to be 4 minutes, with a 


variance of 8 (minutes)*. It is felt that the Erlang would be a reasonable 


assumption for the distribution of the teller’s service time. Also it is assumed 


that there is no limit on the number of customers joining the queue. Bank 
officials wish to know, on the average, how long a customer must wait before 
service, and how many customers are waiting for service. 


Here A = 10 per hour, 1/4 = 4 min. = 


ig hr. 


Therefore, 2 = 15 per hour. 
For an Erlang distribution, we have variance — 1/(ku?) = 8( 


min.)?, 
where 1/z = 4 min., so that k = 2. Hence we have an M/E,/1 queueing 
system. Therefore, 


average time in queue = W, = Ary iui. =e pe ae: 
2k yu(u—A) 4 °° 155 
= (1/10) hr. = 6 min. 
and average number in queue — Lg =AaW, = 10x = 1 customer. 


Ex. 5. In Ex. 4, show that the answers are independent of the fact that 
the service time distribution is Erlang. 


This follows from the fact that for the system M/G/1, Lg, and therefore 


W,, depend only on the mean and variance of the general service time dis- 
tribution. Thus, if the mean and variance are given, as in the case in Ex. 3, 
the form of the distribution will not matter. [See 3.6(34) 


‘| In fact using 
the results for M/G/1, we have 


p?+1o2 
2p) 
Here p=Au= . , A = 10 per hour 
and of = 8(min.a — _.8 Rind hack 
(min.) 60 x 60 (hr.) 450° (hr.)? 


[ Sec. 3.9 
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4 j 
ronnie Lit b eee 
So, | Ne 450 
a 2(1/3) = 1 customer 
and W = glk (1/10) hr. = 6 min. 


3.10 The System GI/E,, |1 — The Matrix-Geometric Method 
The matrix-geometric method [Neuts (1981)| provides a technique for 
analyzing queueing systems with fairly general assumptions. In what follows 
we use this method to study the System size behaviour of the system 
GI/E,/1. . 


erg eee ys : 
Starting, for simplicity, with the case m = 2, we consider the queueing 
1 = * ~] 
system GI/E,/1, ie. a single-server System with general inter-arrival times 
of two exponential phases, marked in 


Thus, each customer, entering service, 
will begin service at phase 2, at the end of which he goes to phase 1, the service 
being completed with the completion of phase 1, 


and with each service time consisting 
the reverse order, say (Cf. Section 3.8). 


We shall denote the state of the system by (n, r), which stands for ‘n in 
system and service in phase 7’, n = 1,2,3,... and r= 1,2. However, when 


‘0’ to denote the empty state of 


n = 0, r has no significance, so we shall use 
the system. Using a notation similar to that of Section 3.4, let k, denote 


the probability that r phases of service are completed in an interarrival time 


of arbitrary duration, and let a = 1— 5 


f= 


k,. Then, corresponding to 3.4 (6) 
0 
we have the transition probability matrix P, given by 


oO FAL 4 ia 8) (8 Ty ee 


_ 


a, at 0 0 0 0 
| 

Oy ky ks k, ky 0 0 ae 

Le ky ke 0 ky 0 0 


k, 


_— 
bo 
i) 
— 
—_— $e , 
R 
i bi 
= 
= 
= 
o 
19 
>= 
ro] 
~~ 
a 
~~ 
_ 
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By an extension of the notation, wherein we use 


Koy Kor 41 
(2) adie ( ) 
Kor_1 Koy 


with 4_, = 0 in Ko, (1) can be simplified to 


(Gy Koy oO oO O ) 
| | 
kG Kio 30 O | 
; 
(3) P= | boo ER, Rie . 
| 
| K5o K, K, K, K, | 
| ; J 
e 
Evidently, in (3), we have 
(4) k,, = (4 &,), a row vector 
and 
Csity : 
(5) i ( wa ) , 2 column vector, 7 = 1, 2,... | 
Aoi | 


To find the steady-state solution for the GI/E,/1 system, we have to determine 


a matrix 
(6) P = (Pp Pi Po ---) 
which satisfies the equations 
(7) pP =P | 
and | 
(8) tie pent 
n=1 


where py, a scalar, is the probability for empty state of the system, the py are 


row vectors 


(9) P, = (Png Pni): n= 1, 2, tre 
with 
(10) Pnr = P(n in system and service in phase r), 7 = 1, 2 


(4 
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and e is the column vector 


an “= (1) 
} 


Now, Eq. (7), on using (3) and (6), implies 


(12) Po = Pot + ¥ Pik, 

i=1 
(13) Py = Poko, + z PK; 

=] 
(14) Pp= 2 pK, 5 = 2,3,... 
i=j-1 

Analogously with 3.4(9), we take a trial solution 
(15) Dy = AB), j= 1, 2, ... 


where A is a 12 vector and B a 2x2 matrix. Note that the scalar Py being 
a degenerate case has to be treated separately. 


Substituting from (15) in (14), we have 


(16) AB/= © ABK, ,,, 


i=j—l 


= AB! 5 BrK, 


r=0 
Thus, Eq. (15) will satisfy (14) if 
(17) = BK,—B 
r=0 
Again, substitution from (15) in (13) yields 
(18) AB = p,k,,+ = ABIK, 
i=l 


= Poko, +A(B—K,) 


using (17). Thus, (15) with B given by (17), satisfies (13) if 
(19) PK = AK, 
Finally, substitution from (15) in (12) gives 


(20) po(l—a,) = ba AB'‘k;, 


i= A ) T Q THING TH HORRY 
NTE ODUCTION oO Qt BUELL 
N INTS z 
= , i Te h e 
‘ the de fini t1o. 1 C f 4 we € 
using ie Vv 
Now ; from (5) y B > I ) 5 1ve 


«© 
— >») K,e 
r=itl 


i act that 
the last simplification following from the fact th 


ic rs 2), since 1 ky = 1 
which follows from (2), siz Re) 5 


Thus (20) simplifies to 


(23) p(l—x%)= 2 = AB'Ke 
( 2 


r=i+l 


= >» x AB'K,e, changing the order of summation 


2 i=l 


=5 A(I—B)“(B—B’)K,e, adding the zero term corres- 
=] 


r ponding to r = 1 
= A(I—B)-1 (B = Ke = BrK,e ) 
1 1 


= A(I—B)-'(B(e—K,e)—(B—K,)e), 
= A(I—B)-(I—B)K,e 
= AK,e 
= Poke, using (19) 
= Polko+ky), using (4) 
which is seen to be true from the 


definition of a. Thus, (15) satisfies (12) 
also, if we use (17) and (19). 


Yec. 3.1¢ P 
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The foregoing shows that 


(15) satisfies the 
through (14), if B ig 


transition equations (12) 
® root of the m 


atrix equation 


(24) = WE, =z 


r=( 


provided further that (19) holds, 


Now, (24) is the matrix equation correspondir 
equation 3.4(14), In fact, analo 
negative solution of (24) 


ig to the well-known algebraic 
Section 3.4, B is the minimal non- 
P< 1(CE£, Gross and Harris (1985), 
computation for B. by successive approxi- 
fe have to find A and Bs, 
alizing condition that t 
in matrix notation can be writte 


gously to 
and is unique for 
pp. 321-322), wherein a method of 
mations is given. Further, w 


To do this, we can 
use (19), and the norm 


he total probability is 1, which 


nh as 
Ls =) 
(25) Pot X& pe =1 
j=l 
or, using (15), 
co 
(26) l—p, = + ABle 
j=l 
= A(I—B)""Be 


From (19), 


(27) A = pk K51 
Eliminating A between (26) and (27) and solving for p,, we have 
(28) Po = [1+ky,K>(I—B)“"Be}1 
Using the value of p, in (27) we obtain A. Thus the steady-state probabilities 
Po Pj, J =1,2,... are all evaluated. This completes the solution for 
GI/E,/1. 
For the system GI/E,/1, the states of the system are denoted by (n, r) 
with n = 1, 2,.., as before, but with r — 1,2,...,m, And, as before, we 
have the empty state n = 0. In this case a service starts at phase m and 
proceeds through phases m-1, m-2,..., 1. The column headings in the matrix 
P in (1) will now read 0, (1, m), (1, m-1), ..., (1, 1), (2,m), (2, -1)5 (2. 1 
(3, m),.... The rest of the analysis is similar to the foregoing theory for the 
case m = 2, For example, now we shall use the matrix K, given by 


heme kins vee Ker s.m_—1 


(29) K, = kms kme ms . ker m—2 


kemp_maa Kmrmig ++ kg 
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Further, we shall have 
ky = (ko ky .» Km) 


(30) 
Ami.m—1 
(31) ky = Amim—2 
Ami 


and the (Ist row, Ist column) element of the matrix P in (3) will be ap _, 


We now have to evaluate 7», Pj (j =e iy Bie es)) WLU Po & 8S alar 


instead of &;. 
xm). Proceeding as for the case m = 2, 


as before and the p, each a vector (1 
we shall see that the solution is still given by (15), but with A a (1 xm) vector 
and B an (mXm) matrix, which is still the minimal non-negative solution 


of (24), where Z is now (mxm) and not (2X2) as before. Thus, formally, 


the results will remain the same as before with the dimensions of the matrices 
and vectors changed. Though A and py will still be given by (27), (28), res- 
pectively, I will be the identity matrix of order m and e a column vector 
consisting of m 1’s. 

Mean Number in System, i.e., L, is obviously given by 


@ 


(32) L = & jpye 
j=l 
= % jABle 
1 


=A > jBI-'Be 
1 
= A(I+2B+3B?+...)Be 
= A(I—B)-*Be 
Also, the second moment, about the origin, of the system size distribution 
is 
(33) E[X?] = = j*ABle 
j=l 
= A[(I—B)“B-+ 2(I—B)-*B’le 


as can easily be seen by summing up the series y j°B/. Hence we obtain 
the variance of the number in the system as : 

(34) Var[X] = E[X*|—L* 

where L is given by (32). 


oR ERE 
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, 4 


The Case of GI/M/\ 


It is interesti : 
’ ; ng as well as : : 
4 well as instructive to obtain results for GI/M/1 as a 


analysis by s etting = ivi O1 

‘. J y zm Li; Kv ide ntly fo the 

case m= 1 we sh: ave — 7 a | 
hall have ¥ e€=t1. K, k,, = ky nd B=hb, the root 


special case from the above 


ae 
of & brky = b lying between 0 and 1 


r - Thus, (28) will reduce to 
(35) Po = [1+(1—b)1b}7 = 1—b 
as is well-known. Similarly (27) will become 
(36) A = p, = 1—b 
and thus (15) will become 
(37) pj = (1—b)bi 


Note. A further generalization of the system GI/Em/1 discussed in this 
section 1s to GI/PH/1, where ‘PH’ stands for a ‘phase-type distribution’. 


Interested readers may refer to Neuts (1981) 


PROBLEM SET 3 


1. Repairing a certain type of machine which breaks down in a given 
factory consists of five basic steps that must be performed sequentially. The 
time taken to perform each of the five steps is found to have an exponential 
distribution with mean 5 minutes and is independent of the other steps. If 
these machines break down in a Poisson fashion at an average rate of two 
per hour, and if there is only one repdirman, what is the average idle time 


for each machine that has broken down ? 
[l hr. 40 min.]| 


2. Inacar manufacturing plant, a loading crane takes exactly 10 minutes 
to load a car into a wagon and again come back to the position to load another 
car. If the arrival of cars is a Poisson stream at an average of one every 
20 minutes, caleulate the average waiting time of a car. 

[5 minutes] 


3. The queueing system of machines needing repair at a repair shop 
behaves like M/E,/1 with A = 6/5 per hour, and # = 3/2 per hour. Find the 


average duration of time a machine has to wait for repair and also the average 


time spent in the shop. ; 
[2 hours; 2 hours 40 minutes} 


anne 
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4. A warehouse in a small state receives orders for a certain item and 
sends them by a truck as soon as possible to the customer. The orders arrive 
in a Poisson fashiin at a mean rate of 0-9 per day. Only one item at a time 
can be shipped by truck from the warehouse, which is located in the central 
part of the state. Because the customers are located in various places in 
the state, the distribution of service time in days is one with probability 
density 4te-. What is the expected delay between the arrival of an order 
and the arrival of the item to the customer ? Service time here implies 
the time the truck takes to load, get to the customer, unload and return to 
the warehouse. Loading and unloading times are small compared with 
travel time. [By permission from Sasieni, Yaspan and Friedman (1959)] 

[74 days] 


5. A colliery working one shift per day uses a large number of locomotives 
which break down at random intervals; on an average one fails per 8 hour 
shift. The fitter carries out a standard maintenance schedule on each faulty 
loco. Each of the five main parts of this schedule takes an average 1/2 hour, 
but the time varies widely. How much time will the fitter have available 
for other tasks per day and what is the average time a loco is out of service 2 


[5:5 hours; 3-18 hours] 
6. In an M/G/1 queueing system, let Y,, be the number of customers 


left behind in the system by the n-th departing customer, and X, the number 
of arrivals during the service time of the n-th customer. Show that 


Y,4, = max(Y,—1, 0)+-X,., = Y,—U(Y,)+Xasa 
where U(x) is Heaviside’s unit function, given by 
U@)=1 (c>90) 
= 0 (x < 0) 
Using the fact that in the steady-state the expectation AY 44] = B[Y,) 


and E[Y;,,]= E[Y2], deduce that the steady-state expectation E[Y] is 
given by | 


ELY] = p+ p1+c?2 


2(1—p) 


where C, is the coefficient of variation of service times, and p is the 
traffic intensity. 


(Hint. Equating expectations of Y,., and Y,, gives ALU(Y,,)] = E[X 
\ n bi oa 


ntl: 


TE aie 
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Now, 


Dm oH 


nal = f X r[(Avjree/r 1) dB 


BUX 
r=1 (v) = p 1 


yrhner w) Is » oC + | * . 
where B(v) is the service time distribution function 


Similarly, 


MXnal = ME +p = "(1 /p*) +-02|+-p 
= Pu +02) p 
Note that xU(x) = x, so that [U(z)? = U(x) 


dln ta aeee : Now the result follows by 
ating the expectations ‘~ F sel y 
Meanie ies I of Y?., and Y?. Assume that ¥ n and X,,, are 


7. In M/G/1, show that the busy period 
given by 


distribution function G(x) is 


~ 


oO 
a {(Av)"e-A0/n 1} G,(a— vd B(v) 


G(x) aa 


oo 8 


where B(v) is the service time distribution and G(x) denotes the n-fold con- 
volution of G(a) with itself. Hence deduce the functional equation 


g(s) = b(s+-A—Ag(s)) 
where 


g(a) = f e*dG(x) 
0 


is the Laplace-Stieltjes transform (L.S.T.) of G(x), and b(s) is the L.S.T. of 

Biv). Verify the expression for g(s) obtained in 2.3. (9) for the M/M/1 system. 
[Hint. Regard the busy period as made up of the service time of the 

first customer plus a busy period initiated by n customers in the system, 

provided that there are n arrivals during the service time. Now use the 

argument of Section 2.3 for the initial busy period. ! 
For M/M/1, b(s) = w/(~+s); thus we have a quadratic equation in g(s), 

giving g(s) = [(s+A+p)—~V/(s+A+p)2—4Ap]/2A. Note that minus sign is 

taken before the radical, since we shall have g(coo) = 0.} 


8. For the system GI/M/2, show that the steady-state probabilities 7, 
are given by 
my = Cb (j > 1) 
where b is the root, in 0 << @ <1, of the equation 
§ = a(2u—2p0) 
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and ‘ 
C=(1 —b)(1—2b)a(#)/[b1—b—a())] 


Tm = (1—b)(1—24())/[1——4(H)] 
and d(s) is the L.S.T. of the arrival distribution A(u). 


[Hint. For evaluation of C and 7, proceed as follows. Considering 


one-step transitions, we have 


(1) Mo = Mo Poo tM Pio + (1—M— 7) Pao 
where 
fe) 
(2) Poo = JS (1—e-*4)dA(u) = 1—a(p) 
0 
te) Pro = J (L—e-M)*dA(u) = 1—2a(n)-+4(2p) 
0 


and Py is the probability that a customer, who arrives to find two or more 
units in the system and, therefore, has to wait, is followed by the next arrival 
who finds the system empty. Let the interval between these two arrivals 
be wu. Then, if the queueing time of the first customer is w, we should require 
both the servers to complete service within a duration u—w. Also w follows 
the distribution given by 3.5(9). Therefore, we have 


i 


oo # 
q =f f 2(1— »—2(1—B) My p—M(U—w) 2754.7 / 
Pu J ! 2(1 b)ne ( ip [l — e—Mu—u | du (LA (u) 


which, on using the relation a(24—2b) = b, simplifies to 


4 ON Ea oe 1—2b ; 
(4) Pao [_ | 2—4da(w)— 5 2p) | 
The results follow on using (1) through (4) and the facts that 7 = Cb) (j > 1) 
and » y= 12) es 
j=0 


9. For the system GI/E,,/1, 


issn by show that ky as defined in Section 3.10 is 


ky = (—I){(mp) [re amp) 


where » is the mean service rate 
=f) 


and 
3 Ir 
amp) = > as) | 
} s=Mp 
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is the r-th derivative of the 
the Laplace transform, 4(8), of the inter-arrival time 


density function, evaluated at 
y ’ aluated at » — mp. Deduce the expressions for k, for 
ssions for k, 


the systems : 
(i) Ey/Em/1, (ii) D/Em/1, 
[ By definition, 


wo 


a f \ 7 j 
ky i (mpn)re mit irVa(u)du 


whence the result follows on using the L.T. relation (9) of Appendix E 
(i) For E,/Em/l, as = ( Lidice sas. 
) | 1/En/ i(s) fees! , whence on performing the successive 
differentiation and simplifying, we have 


I+-r—1, 


ky = ( ) p'gr 


fi 


with 
p=V|(lA+mp), gq = 1—p 
which is a negative-binomial distribution. 
(ui) For D/Ep/1, a(s) = e-s/, Proceeding as before, we have 
ky = (mp]A)re-m4l4 Ir! 
which is a Poisson distribution.] 


10. For Ex. 2 of this chapter (p. 67), evaluate the average number of 
patients waiting in the doctor’s office. Also obtain the average number 
waiting, when : 

(i) the inter-arrival times have a 2-Erlang distribution with the same 

mean as before, 


(ii) the inter-arrival times are constant and equal to 20 minutes, Note 
that this is the case when the doctor makes firm appointments and 
the patients adhere to the same. 

[2; 1-6636; 0-5362] 
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CHAPTER 4 


BULK. TANDEM AND PRIORITY QUEUES 


4.1 Introduction 
In this chapter we discuss queueing systems 1n which either the alrivals 


take place in groups or the service ip: 
a random variable given by a probability distriby- 


is rendered in groups. The size of the 


groups will be regarded as 
Hon. Tie queueing systems described in earlier chapters with group size 
always equal to one become particular cases of the models discussed here. 


Also included in the chapter are queues in tandem and priority queues. 


4.2 The Queueing System M*/M/1 

This is a single-server bulk-arrival queueing system where the times 
between consecutive arrival epochs and the service times are both exponentially 
distributed with means 1/A and 1/u, and, at each arrival instant, customers 
arrive in groups of variable size X, given by the probability distribution 


(1) PX =r) =a, P=1,2) 8, oe 


It is easy to see that in the steady-state, corresponding to Eq. 1.5(7) and (8), 


we shall have 


(2) 0= —(A+p)P,-+A y P, 1dr +P 41, n= 1,2,... 
r=1 


(3) 0 = —AP, +P, 


Define the following generating function : 


wo 
(4) P(z) a bar Pz", |z] < | 


n=0 


From Eq. (2) and (3) we have, ag usual 


(5) 0= —A+mPe)+uPrte Ed Petras FP yaen 
n=0 n=1 r=1 P 
The double summation in the last term, on changing the order of summation, 
gives 
(6) 3 5 
BS (Pn )ayar = P(z) A(z) 


r=1 ner 


(4 
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7 ia ia 
fel { . 
is the generating function of ’ ) 
n of the group size, ‘Thus (5) finally yields 
4 ot 
(8) P@) = __ “(1—z)P, \u 
NeA(2)—(A+ pep } 
i ; 
Now, since P(l—0) = 1, on using l’Hpital’s rule we have from (8) . 
where @ = 2 ray is the me roti ais 
waitin ! mean group size, and we use the following obvious i 
(10) A(1) = ] : | 
(11) Atiyes See Sg 1 
dz 2=1 
It is clear that the traffic intensity for this system is given by | 
So (9) shows that we have the usual relation : 
: 


(13) P15 
Thus, for the steady state to exist, p < 1. 


Substituting from (9) in (8) and then evaluating the derivative at z = 1—0 
with the help of l’Hopital’s rule, we see that the mean number in the system 
is given by 


(14) L = P(1—0) = pop (1+ 4") 203 
where 

A") =F 4@)| 
(15) — z (r?—r)ay 


o® being the variance of the group size. 


ll 
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From (13) it follows that E [Nserv] = P, and hence 


(16) L, = L—p 
In the case of single arrivals we have 
Af) =z, @=1, A*(1)=0 


and, therefore, (14) reduces to 
A p 


(17) be ee 


which is a well-known result for M/M/1. 


Waiting Time 

We now consider the distribution of the time spent in system by an 
It may be noted that in the system M*/M/1 arrivals 
We shall assume that within an arriving 
Let us first take the 


arbitrary customer. 
occur in groups of variable size. 


group the customers are served in random order. 


first-served customer of an arriving group. If on his arrival he finds 
customers in the system, he will stay in the system until the completion of 
(n+1) services. Therefore, the waiting time (in the system) of the first- 
served customer of a group has the density function f,(w), given by 


(uw)eme 


(18) f(w)dw = = P,: ——— . pdw 
n=0 nm: 
Hence, the Laplace transform of f,(w) is given by 
(19) Fils) = f e-*’f,(w)dw 
0 
oe oo / ph cae be a 
Fi a Pa | p+s pbs ie ( pits 


where P(z) is defined by (4) and is given by (8), (9). 

Now, if a customer is served as the j-th customer in his group, he will 
have to wait for an additional time during which (j—1) services are completed. 
Also, it can be shown that in random selection for service the probability of a 


customer being served j-th is 
j=l 


(20) 2 = (1— 2 a)/a 


r=] 


Therefore, as before, the L.T., fils), of the p.df. of the additional waiting 
time of an arbitrary customer over the first customer is 
j-1 


3 He = Sof | 


1 1 
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On substituting for Pp; from (20) in (21), we have 


oni 1 ~ ee 
(22) he= 4S (ia Se, ake ie 
@ 4-1 \ a pets } 
ee 1 © j-1 
"Ty ai ae {1 | £ = 
ee esd eS j=-2 r=1 fe+s } 
: hs a 
i ral oe p> = a ( : fe dr 
¢ + r=i jorsa’ \ pte pts } | 
] pets + 
ly: Be » pb . 
a F | Mt )1 
Hence, using the convolution theorem of Laplace transforms, we have for | 
the L.T. of the p.d.f. of the waiting time in the system for an arbitrary ij 
customer : 4] 
" ae 
(23) Fis) = fuerfe) il 
i} 
| 
ae By EE me 
as E ( fis )] 
: : 
Differentiating f(s) w.r.t. s and taking limits as s—»0 with the help of ; 
PH6pital’s rule, we have for the mean waiting time in the system 
(24) W = —fo) = —[L414()a"a)] 
=a | (7) ed 


If we eliminate A”(1) between (14) and (24), we obtain 
(25) L = \aw 
which shows that the well-known Little’s formula holds in this case also. 


Note that here the mean number of customers arriving per unit time is Ad. 


Queueing Time 
Proceeding as above, we can show that the L.T. of the p.df. of the time 


spent in the queue by an arbitrary customer is 

(26) fals) = Farls)fals) 

where f,(s) is as given in (21), (22), and fals) is the L.T. of p.df. of the 
queueing time of the first customer of a group, given by 


n 


(27) farls) = = P,| rie an | es 
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Also, we have for the mean number in the queue 


© = e 
(28) Lg — » (n—1)P,, = » nP,,— Dy ie 
n=l n=0 n=1 


= L—(1—P,) = L—p 
And the mean queueing time is given by 
(29) W, = —f; (0) 

From Eq. (26) through (29) and (14), we can show that 


(30) Ly = AGW, 


which verifies the other relation in Little’s formula. 


4.3 The Queueing System M/M?/1 


This is the single-server bulk service queueing system with Poisson 
arrivals and exponential service times, in which whenever the server finds 
a number of customers waiting, he begins service on b of them or the whole 
queue, whichever is less. Let P,, be the steady-state probability that n units 
are waiting in the queue (n = 0, 1, 2, ...) and the server is busy, and let P,, 
be the probability that the server is idle (the service station is empty). 


The steady-state transition equations are 


(1) (A+p)P,, ow AP, a t+HP rtp (n 2 1) 
b 
(2) (A+p)Po = APog tp & Py 
j=l 
(3) 


AP oo = "Po 


Let us introduce the generating function 
22) 
(4) Pe) = = P,z", |z| <1 


From (1) and (2), using (3) we have 


b 
(5) (A+p—Az—pz)P(z) = pw X P,(1—zn— 
n=0 
whence 
b-1 
# x P,,(e*—z?) 
(6) 


le n=0 
0 = Opp 
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Now, let us take 
(7) f@) = —(A+-p)zb, q{(z) = Agbt14 p 
and take a closed contour 0 given by \z| 

we have, on C, 


= 14-6, where 8 is small. Then, 


(8) \f()| = |(A+-p)z? | = |A+p|+|2]e Es (A+ y)(1+42)P 
and, = (A+-)(1+-b8)-+-0(6?) 
(9) \g@)| = [Ab | < IAl-fe}+24. [ut 


= A{L4+-(64+1)8}4 p+0(82) 


= A+p+Ab +-1)8+-0(6?) 
Therefore, we can write 


(10) \9@)| < \ff| ond 
if 
(11) Ad < phé 
ie., if 
A 
(12) ere Pies 1 


which is the condition for the steady state to exist. When this condition is 
satisfied, using Rouché’s theorem we can write that the denominator of the 
right-hand member of (6), which is f(z)-+g(z), has the same number of zeros 
inside the contour C as f(z). But f(z) has b zeros inside C: so the denominator 
f(z)+g(@) also has b zeros inside C, one of which is z = 1 and therefore the 
remaining b—1 zeros are inside the unit circle \z| =1. Hence the denomi- 
nator has only one zero outside |z| = 1; let us call it zp. 


Now, corresponding to each zero z; of an expression, the expression should 
have a factor (z—zi). Therefore b of the factors of the denominator on the 
right-hand side of (6) should cancel with the b similar factors of the numerator. 

= . . 
l i oe : inate of degree 6+1. 
Also the numerator is of degree b and the denominator is of deg 


So. on cancelling the like factors, (6) can be written 


A rs) gh 
\= = A D 
(13) Oo ge 


where A is a constant. The expansion given in (13) is valid since clearly 
\z/z)| <1. Noting that the coefficient of z® in P(z) is P,, we have, on setting 
0 : 


eu, 


(14) A= Pio 
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Thus, (13) becomes 

i (15) 

i Therefore, we have 
(16) 

Also, from (3), 

(17) 

Thus, (16) yields 

(18) 


Po = (A/u)Pop = bpP oq 


P= bpP oz", m= 0; 1,2, ... 

Finally, P,, can be evaluated from the normalizing condition 
@ 

(19) Po+ 2 P, = 1 
n=0 


which, after simplification, gives 


—1 
20) P= Sea 
( © (iat 
This completes the steady-state probability distribution of the number in 
/ the queue. 
a The mean number in the queue is 
| oe I e's 
: 21) lg = E nP. = Pp {| * za t-.) 
| ( . par = tg 2% 
Ete 
(@—1)? 
Using (17) and (20) we have 
(22) j ed bpzy 


(2o— 11+ bp)z.— 1p 
Note that (21) could also be obtained by differentiating (13) atz = 1, 
It may be observed that z, 
of the equation 


(23) 


is the root, outside the unit circle Jz] = 1, 


APH (A+ n)ab4 y = 0 
which can be written 


(24) ie A 


2 Mpa 
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1 
Lotti z= —, We see ; ; a 
Setting ae aus I/z is the root, within the unit circle, of the 
equation 
A 
(25) t= 


A+ p— pnd a G(— pa) 
where @(-) is the L.T. of the inter-arrival time density function, since we have 
(26) G(s) = Al(A-4¢) 
Particular case: M/M/1, Setting b = 1, it is clear that z= p/A = 1/p 
and we can verify the 


Various results of M/M/1 from the foregoing. In 
particular, from (20) and (18) 


(27) 


, 


Poy = 1—p 


(28) P,, = P(n+1 in system) = (1—p)pn+1, y 


ah, 2... 
and, from (22), 
(29) Lq = p*|(1—p) 
Queueing Time 


Since the waiting customers are taken in for service in batches of b 
customers or the whole queue, whichever is less. a customer who finds (ib+-j) 
customers waiting in the queue (0 < j < b—1) has to wait until (i-+1) service 
times are over—i for the ib customers preceded by one for the residual 

of the batch being served. Therefore, the queueing time density 
fa(w) is given by 


service 


, function 


(30) falwydw =% = Pars MOD e-woudr 
(31) P(w = 0) = Py 

Substituting from (16) in (30), we have 

(32) falw)dw = 2 Pept OY e-sudo 
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Thus, the mean queueing time is 


eo 


(33) Wq= J w.fg(w)dw 
Wee ee at 
= Po gost * a—25") 
pl —252)(1—25°) 
Now, 2) is a root of (23); therefore, we have 
(34) a? = (A+ —Azy)|fe = 14+bp(1—%) 
Hence, (33) simplifies to 
Pe 
io 0~0 i: 
85) «= “ue —1ybpG—1) 
Poo _ Lala 


~ bp (&—1)? 
using (12) and (21). Thus we see that Little’s formula holds for this system. 
4.4 Queues in Series (or in Tandem) 


Let us consider two service facilities arranged in series as shown in the 
figure. Let customers arrive in a Poisson fashion with mean rate A for 


0...00]0|— 0...00]0]> 
ie A 


service at service station 1. After completing service at service station 1 
the units join another queue for service at service station 2. Let the service 
times at stations 1 and 2 be both exponentially distributed with mean service 
rates i, and 4, respectively. This system thus consists of two distinct queue- 
ing systems which we shall call subsystems 1 and 2. Let ey no!) be the 
probability that at time ¢ there are n, units in the first subsystem (in queue 
or in service) and m, units in the second subsystem. Then, we have 


(1) Pa, ng tAt) = Pai nolt{1—(A+oa tHe) Ab} 
+P at, ng) AAP, 4, no-ilé) Hat 
TP notilt)'MgAt+-o(At), (n, > 1, m > 1) 
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. DE 
whence we have, for the M AND PRioR; 
steady- TY QUEUES 


State Probabilities Pp : 
(2) 0= Atm a)P, n +AP pi 
, 1g “4, at 
ny 1 Ng MP act ny~1t bP, Rgtl 
5 *S 

rye is 1 > f i 
Similarly, we have the following Special cases of thi er : : 
a 18 equation : 
Y (A+ AP, oFAR,, a gee 
(4) 0=—A rere.) (my > 1) 
= —(A+p,)P 

Hs) 0 ng tyP, mathe, (m, > 1 : 

; Neri a | 

(5) 0= AP 9+ MP o, | : 
We take as a trial solution : 
) 
(6) 
: 


Pat my = Cpypr? 
where C is a constant and Pi = Alp, p AJ 
Ly Fe = Aj py. 
Dividing (2) throughout b 
y A and substituti 
suituting from (6), we have 
: 0 ree 1 ay a n n,— 
(7) OTe POPP Cpl pte ptt. 6 mate 
a | : 2 TrUpr"po 
Her is true. Hence the trial solution (6) satisfies (2). Simi i 
also be seen to satisfy (3), (4) and (5). Also using the A lg ais 
; zing condition 


So i 


n1=0 Mo=0 M1, Ny 
we have 
() C = (1—p,)(1—p,) 


and thus we arrive at the solution 


10 9 
( ) Pa neg a (1—p,)(1—p,)pi4p5 
The marginal probability of n, units 


in the first subsystem irrespecti 
i “J Spectiv 
the number in subsystem 2 is pective of 


co 


- x nz a (1—p,)p,™ 


No=0 


(11) P 


Ni,e 
Similarly, the marginal probability of my in subsystem 2 is 


(22) P == Py, og = Pee 


12 
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From (10)-(12), we observe the following : 


(a) The numbers of units in the two subsystems are statistically 


in- 
dependent. 


(b) The first subsystem is an M/M/1 system with traffic intensity p,. 
(c) The second subsystem is also an M/M/1 system with traffic intensity 
/», and since the mean service rate is 4g, the mean arrival rate should 
be A. 

Now, the observation (b) above is not at all surprising, since subsystem 1 
could be studied separately, without caring as to what happens to its output, 
as M/M/1. On the other hand, (¢) is noteworthy, since subsystem 2 does not 
experience free input and A is given as the mean arrival rate at subsystem 1. 


The observation (c) implies that the inter-departure times at subsystem 1 


are exponentially distributed with mean 1/A. This important property 
about the output of a Poisson queue will be discussed later. 


From (11) and (12) we see that the mean number in the first subsystem is 
L, = p,/(1—p,) 

and the mean number in the second subsystem is 

(14) L, = p2/(1—pe) 


Therefore, the mean number in the entire system is 


(13) 


(15) Bi Dy fLig er et Ee 


ij, te 
Eq. (15) can alternatively be derived by noting that with probability 
Po ,, there are n,+, units in the entire system. Thus, 
1 "2 


(16) C= 5 s (ny +M2)P,, 


Ne 
nj=0 No=0 » 2 


Now, using (10) and evaluating the double summation we have (15). 


In fact, the observation (c) above implies that the two subsystems can 
be studied separately as M/M/1 systems each with mean arrival rate A. 
Further, if there are any more subsystems added after the second which have 


exponential service times, they will all experience Poisson arrivals with mean 
rate A. 


Ex. 1. In an M/M/1 queue, show that the time between successive 


departures has the same probability distribution as the inter-arrival time 
distribution. 
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mean arrival and 


Now, with probability > 
+ Gee siew cs ne a departing customer leaves behind 
service time. Also, with probability Pes 
leaves an empty station, and in this ae 
consists of the time until the 


more, by the Markovi 


91 
Let A, w be the 


service rates, and let; p= Alp. 


to next departure is the same as a 
= 1—p, the departing customer 


the time to the next departure 
next arrival follow 


ed by a service time. Further- 
4 property of the arriy 


3 al distribution, the time to the 

Arriv 0 F > 7 = - ‘ zi ; 
next irrival has the same distribution as an inter-arrival time. Thus, in 
ae sacORE Lon the inter-departure time is distributed as the sum of an 
inter-arrival time and a service time 

lA "2 ry situ ; + 

Cherefore, the density function f(t) of the inter-departure times is 


f(t) = (1—P,)b(t)-+ Py-a(t)sb (2) 


are the density functions of the inter-arrival time and 
service time, respectively, and the asterisk ( 


*) stands for convolution. Thus, 
we have 


given by 


where a(t) and b(t) 


t 
i} AeWAu « pe Bit-Wday 
0 


A 
f(t) == au west (1— 2) 


= Ae-*t +. A(u—A)e-Ht { e4—Audy 
0 
elt-At_} 
— Ht | | wt ae eae a 
re [ ee | 


= ret 
which is the same as the inter-arrival time density function. 


4.5 Non-Preemptive Priorities 


The concept of priorities was introduced in Section 1.4. We shall now 
consider the waiting time (in queue) of customers in the case of non-preemp- 
tive priority or head-of-the-line priority rule, where a service once started 
is allowed to continue until completion, as explained in Section 1.4. 


Let there be r priority classes, called classes 1, 2, ..., 7, class 1 having the 
highest priority and class r the lowest priority. Let there be a single aah 
and let customers of priority class k arrive according to a Poisson distribution 
with parameter A, (k = 1,2,...,7). Thus the arrivals in different, classe 
have different Poisson distributions. The combined arrival distribution will 


be Poisson with parameter 


(1) ei 


= 
x Ax 
=1 


(4 


[ See. 4.5 toc, 4.5 
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Let the service time distribution for class k have a probability distribution 
r) with mean I/u,, so that the combined service time 


whence, solving for W, we have 


Br(t), (k = 1, 2,..., 


D-1 
eed hiliiea ae 0 gee [ : Fears ] 
distribution is (10) D = PEW x4 K(Ty) \ (l—ay) 


; 
anene 
/ (2) Bt) = crate Ap Br(t) Now, we shall prove by induction that the solution of (10) is 
We assume the queue discipline within a priority class to be first-come firgt- (11) Wee ET) / 
(1—op_,)(1—a,) 


served. 


" py j rs eto yg alaes al i 2 R . . 
The queueing time of a customer of any class will be zero if he finds the For, taking the solution (11) to be true for Wy (k = 1, ..., p—1), we have 
server idle and therefore no customers waiting, on his arrival. To find the : 
mean queueing time of a customer of a priority class p, say Wy, suppose on 


arrival he finds the server busy and nz units of priority class & (k= 1,2,..., p) 


P-1 : p-\ 
(12) x prWe = HT) & ot 
kel kaa (I —FK_1)(1—o) 


waiting in queue. Let 7 be the queueing time of the customer of class p ee 
under consideration, During time 7’, let 7%, units of class k (k = 1, 2, we) P—1) = K(T,) = [ — = ; |, since py = O,—OK_, 
arrive, who will all be served before the unit being considered. Let 7, be ee ts Ares 
the time required to complete the service already in progress, 7; the time to 1 
serve the nj, units already waiting (k = 1, 2,...,r), and 1; the time to serve se | 1—y —1| 
the n, units (tk = 1,...,r). Then, we have ae 
are Substitution from (12) in (10) shows that (11) holds for W,. Also, setting 
(3) f= rt T+. Tr+T, p = | in (10) and noting that an empty sum is zero, we have 
k=l k=l 
' E(T,) E(T,) 
Taking expected values on the two sides, | (13) ae eo (pao) 
2 D-1 ; D : , 
(4) Wy = E(T) = X E(T;)+ ¥ B(T,)+H(T,) | which shows that (11) holds for W,. Thus (11) is established. 
k=1 k=1 | 
| é 
Now. we have Now, (13) can be written 
(5) (Tr) = (1/pee)E (mg) = (1/up)(Ax We) = pW (14) Wy = EKA ps) 
where we use the notation which gives the mean queueing time of a customer of priority class 1. ‘This 
(6) Pe=Arcie, k= 1.2 , } should obviously be the same as the mean queueing time Wg in an ordinary 
— “AIR, O— 1, Ze F: sammnis j irst-c > first- 
Also , M/G/1 queue without priorities, the queue discipline being first-come first 
* ; | served. This queueing time was shown in 3.3(40) to be 
(7) E(T;) = (1/pt)ALE(T) = Wp ; | AB(v?) 
Using the notation (15) a 2(1—p) 
(8) 5 
= ke! Ps %=0;k=1,...,¢ Comparing (14) and (15) we see that 
i= 
na 
we have, from (4), on using (5) and (7) (16) R,)= (A/2)E(v*) = (A/2) J wd Biv) 
‘~ : eit > service time 
(9) Wy = Op Wot & PeWre+E(1,) the expectation H(v*) standing for the second moment of the serv 
oe distribution about the origin. 
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time W, of the priority class p. 

The following direct proof of (16) is interesting. 

Measuring time from the commencement of the service in Progress, th 
= j A ; Z : o“SS, the 
; probability of an arrival in (¢, f+-dé) is Adt and then the residual service tj 

Pe Sar 7 : : ; ime 
up to completion is y—¢, where v is the service time. Also dB(v) 


a is the pro. 
| bility that the service time lies in (v, v-+-dv). Hence, sls 
(17) E(f,) = f Jf w—tAdtdBiv) 
v=0 £=0 


= (A/2) f v°dB(o) 
0 


which is the same as (16). 


Tt may be noted that for (10) to give a positiv 


e value of W 
have Cp <1. 


aan p We should 
The above discussion is in fact a steady-state analysis W 
conclude that the steady-state will exist if age : 
18 
(18) Or = Pit Pat .-.pr < 1 
which is obtained by requiring that Wy for k = 1.2 r are all positi 

- pees sitive, 


PROBLEM SET 4 


1. For the system M~/M/co 


ee. , show that the steady- 
So S 


state equations are 


0=—(Atny)P.1A ¥ L(n- 
(A+ m)P,+-A i ce ar+(n+1)uP, 
0 = —AP,+uzP, 


01 


Hence show that the p.g.f. Piz) = 5 Pz" is 


given b 
n=0 = y 


ae 


P(z) = exp | Aj) 1—A(u) 1u| 


2. From Problem I 


, deduc ; 
re P. educe that, for 


= e-4/4/) / 
e- 7/4) /t)*/n}. 


3. 

ie In Problem 1, by taking the left-hand 
a at the mean number in the System j 
size. Also find the variance of th P 
once again, 


M/M/co, P(e) = expl(A/u)(z—1)), 


derivative of Pi) at z = 1, 
8 Ad/u, where d is the mean group 
€ number in the system by differentiating 


[(A/u)(a-+A"(1)/2] 


Set 4 
On substitution for (7) from (16) in (11), we have the mean queuein 
& 


OP 


Problem Set 4 | 

BULK, TANDEM AN 

AND PRIORITY ead 
4. Ina bulk arrival RITY QUEVES ¢ 

by 4.1(1), let p,,(t) 


(0, ¢]; show that 


; dt Patt) —— —Ap,(t)-+-A z Pn—r\t)ay, (n = 0,1, 2 ) 
bet pay My wee 


Hence prove that the p.of = 5 i 
p-gt. Pe, t) — Ps, P, At)" is given by 


Pe, t) = exp[(A(z)—1)Ag] 
Deduce that the mean number of arriv 


als is Adt, 


| a being the mean group size 
Show also that the variance of the number of arr ; he | 


ivals is At(a+-A"(1)). 

[It may be noted that the results of Problem 4 
service time distribution or the number of 
for the system M*/G/c.] 


do not depend upon the 
servers and will, therefore, hold 


5. In Problem 4, show that from P(, t) we can obtain the probabilities 


p,(t) as 
Polt) = e-*# 
. At , 21 
Pall) = (5) 2 (nia, spilt), (n > 1) 


differentiate partially w.r.t. z both sides of the equation giving P(z,t) and 
equate coefficients of z* on both sides.] 


[po(t) can be obtained by setting z= 0 in P(z,t). To obtain p,(t), n> 1, 


NOTE. The equation P(z, t) = exp[(A(z)—1)At] of Problem 4 is an 
important result in the theory of Stochastic Processes. It may be noted that 
P(z, t) is the probability generating function of the counting distribution of 
events which occur at time points (epochs) distributed according to a Poisson 
distribution such that at each occurrence epoch a random number X of events 
occur with P(X =r)=a,. Such a process is called a Compound Poisson 
Process. In the special case when a, = 4,,, the Kronecker delta, so that at 
each epoch only one event occurs, we get the Poisson Process. 


6. For the system M/M/1, where arrivals are in groups of a fixed size 
R, write down the steady-state transition equations and hence show that 


_ R+ Ra 
the mean number in the system is reap n—BX* 


—_ 
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7. Show that the number of customers in M4/M/1 behaves like the 
number of phases in M/E,/1. Hence i aha ie mean number of Phases in 
the latter (Erlangian service) system is => ‘ er , Where p is the traffic 
intensity of the system. 


[Note that in order to have the mean service time = 1/u in fe) second 
system, each exponential phase must be completed in an average time 1/Ru. 
Hence the « of Problem 6 should be changed to Ry, and the result follows 
by setting A/u =p. (See Eq. 3.8(31).)] 


8. In the queueing system M*/M/I, let the arrival group size have the 
geometric distribution 


ay = gpt!, gq = 1-2, (r=1, 2. ss.) 
Show that 
P, = 1—p 
n = P(1—p)(p+9p)*", mB 1 
Hence obtain the mean and variance of the number in the system. 


re _P __, P2p+(q—p)p} 
q(1—p) ’ g(1—p) g(1—p)? 


9. From the results of Problem 8 deduce the corresponding results (for 
P,,, mean and variance) for the system M/M/1. 


[Take g = 1. (Why ?)] 


10. Consider the following 


queueing system, which may be termed 
M/M®-</1 (Cf. Kambo and Chaudhry 


(1985)) and is applicable to queues for a 
bus at intermediate stations (different from starting station and destination) 


The server is always busy and on the completion of each service he starts 
# new service run, taking b—X customers from those waiting, or the whole 
queue, whichever is less 


, X being a random variable distributed as 
P(X = 1) = & (r = Ly eee 
Thus, a service run may start taking ze 


to from the queue if none are waiting 
or if X=b. (In the latter case the server is full to capacity.) Define the 
steady-state probabilities 


P,, = P{n units in queue (at the bus station)} 
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thes - S . . 
Show that the following equations hold : 


' b 
(A+p)P,, Rc AP, atu z Cb_2P ni (n 2 1) 
r=0 


where 


Using generating function and Rouché’s theorem show that when 


b 
p= Aluib—e) <1, (¢= re, \ 


r=0 
we have 


P,, = (l—a)a" (n > 0) 


where « is the unique root, inside the unit circle |z| = 1, of the equation 
b 
pal Y oz7—(A+pye+A = 0 
r=0 
11. Consider and interpret the special case of Problem 10 when b = 1. 
12. For the queueing system M¥/D/1 with constant service time b, show 
that P,, the steady-state probabilities for the number in the system, are 
given by 


n+1 
oe a (Po+P,)p,(6)+ 2 Pip,i+1(), (n 2 0). 


[Proceed as in Section 3.7.] 


Note that the above equation can be used recursively to evaluate P, 
(n = 1, 2,...) in terms of P, and p,(b), (n = 0, 1, 2,...). Now, the p,(b) can 
be evaluated using results of Problem 5 with ¢ replaced by b. And, P, can 
be evaluated from intuitive reasoning by equating the average input and 


output rates of the system (see, for example, Chaudhry and Templeton (1983), 
p. 72) as follows : 


Ad = (1—P,)(1/b) 
whence 


P, = 1—Adb = 1—p 


*Note that (b— ) is the average intake by the server. 


13 
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, ' ™ ) Dy , ; 
13. Using the equation for P, given in Problem 12 and the Benorating 
} function P(z) : P,2", and proceeding as in Section 3.3, show that foy the 
; nee 0 CHAPTER 
; M*/D/1 system 5 
A) (1—p)(1—z) APPLICATIONS 
(2 1—z exp[Ab(1—A(z))] 
5.1 Introduction 
: ; 3 , ri NO Syste 2 
: Hence find the mean and variance of the number in the system. In this chapter we gh; j Sande ; 
} - c > Shall discuss applications of queueing models in 
; 9 m1) Le oblem solv a eee : ; 
: = 2 pA"(1)/a 3 : ; solving and decision making in various types of human activitie; 
; Mean 4 aT age though apper ance . “4s ‘ ane : nis ate 
| 2(1—p) a Pig sheds ep of queues or waiting lines (physical or, less discernible 
Vantanos = VartN} eo P"(0) ene ce: ee 18 & universal phenomenon. thus making queueing iene 
aris r€ l - 4 4 y f , 7 r hc ee . y 307 * 4 ; j ad . = 
app oot oe 10 almost every walk of life, ranging from daily routine jobs to 
here complicated production processes, yet there are some areas which ha 
+ys . oe ia ~ t au AVe 
6pA"(1) 3p%A"(1))* (1) traditionally been taken ag particularly suited to such an application. It is 
7 » 2 3_1. pt-} 2 — -20 - : hese areas like tele > traffic. machine ; 
P"(1—0)= 61 —p)2 { 6p 10° -|-p*- a } i -2p(1—p) 3 th reas like telephone traffic, machine interference, 


inventories, scheduling, 
18-d-vis application of queueing theory. 
are apparent from the examples and problems 


ete., which we shall he discussing » 


Cf. 3.3(54), et seq. Note that for M/D/1, A’(1) = A’”’(1) = 0 In fact, several applications 


and hence we have in the foregoing chapters. 
Pr 0) pp? 4p |-6) 5.2 Two Types of Models 
6p)? Models, whi > defi i 
Models, which may be defined as abstractions of real-life phenomena 
q., & set of differential-difference equations is a 


used to facilitate analysis (e 
whieh tallies with 3.3(60). | mathematical model of a queueing situation) 


, are classified in several ways. 


One such classification divides models into (a) descriptive, and (b) normative 


(or prescriptive) models. The former type of models describe the state of 
affairs of a system on an ‘ 


14. Notice the similarity between the relation given in Problem 18 and 
Eq. 3.3(28). How would you account for it ? ‘as is’ basis, while the latter tend to seek a norm 
Thus, all “optimization”? models are normative. As 
should be evident from earlier chapters, most of the models used in queueing 


theory are descriptive in nature inasmuch as most of the time we seek to 


5 


or “as it should be’’. 


study the queue length, the waiting time or the busy period through pro- 
bability distributions. However, consideration of business and industrial 
applications has brought in the concept of optimization in queues. 
Optimization. In optimization studies in queueing theory, generally 

arrivals are taken as uncontrollable or uncontrolled, and we try to find the 
optimal service rate or the optimal number of service channels. In such 
studies we consider two types of costs—(i) the cost of providing service 
including the cost of the server’s idle time, and (ii) the cost of the customer's 
waiting time. As we increase the service rate (say 2), the first type of cost 
increases and the second decreases; (see the figure on the next page). We are 
interested in the lowest point in the total cost curve, which corresponds to 
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: value of w, denoted by “*. This is the general format of the 
the optimal value of /, tach, entjentnto slight variations depend. 
We now proceed to describe applications of 
though such a description can never 


optimization studies in queueing 
ing upon the problem studied. 
the theory of queues in various areas, 


be exhaustive. 


cost 


Service 
cost 


Total cost 


Waiting time cost 


5.3 Applications to Telephone Traffic 

As pointed out in Chapter 1, telephone congestion problems accounted 
for some of the initial applications of queueing theory. We start with some 
basic concepts in telephony. 

A full availability system is a set of devices (telephones, trunk lines, etc.), 

which are all accessible to a group of sources of calls (telephone subscribers), 
Otherwise, we have a limited availability system, where different groups of 
sources have only a limited number of devices available to them. A hybrid 
between the two is a link system, where linkages exist between groups of 
sources and destinations with limited access available between other combina- 
tions. A common control system is a loss system where all calls arrive at a 
central control before going into service. Thus, a call is lost when the control 
is busy, even if a service channel is free. In a loss system there is no arrange- 
ment for a call to wait; it is served on arrival or is lost. Alternatively, we 
have a delay system, in which it is possible for a call to-wait until a free channel 
is available for its service. Most systems studied in queueing theory are 
thus delay systems. 


We shall begin this description of queueing theory applications in tele- 
phone systems with : 


(a) Loss Systems. Let calls arrive in a Poisson fashion with mean rate 
A and let the service times (usually called holding times) be exponentially 
distributed with mean rate #. Further, let there be ¢ telephone channels 


Sec. 5.3 | 


es APPLICATIONS 191 
with no provision for waiting, 80 that a call is lost 


busy on its arrival. Thus, we lines f it finds all the channels 


the queueing system M/M/e/e, which has 
1.6(e) with the 
may 


already been discussed in Section 
parameters. Alternatively, we 
equations as follows : 


help of state-dependent 
proceed by writing the steady-state 


(1) (A--np)P.,, — AP, +(n -4 uP, a at 
(2) ARS = BP, 
(3) cuP, = AE 5 


where Py, = P{n channels are busy}. 


As before, the solution is 


eae ape iges 
i") ia 1+p+py2\+.4+pjer? P=AlM, n=0,1,...,¢ 
(4) is called Erlang’s loss distribution*. For n= c. we have Erlang’s loss 
formula. When c— o, (4) becomes 
(5) iE = ePprln ! 


giving a Poisson distribution for the number of busy channels, as we observed 
for the system M/M/co (Section 1.6(c)), 


If we take the mean arrival rate to be A; when i channels are busy 


(i = 0, 1, ..., ¢) and write A4/~ = py, then it can be shown as in Section 1.6 that 


— . | ! 
(6) P,= Py PoP «++ P/N } 
where P, is obtained from the normalizing condition 
c 
(7) ZP,=1 
n=0 


Now, the number of lost calls per unit time is 


(8) AcPc = bPckc = UPoko(PiPs --- Pele !) 
And, the total number of calls arriving per unit time is 
e e APY sax Pai 
(9) EAP, =~ B py Py af = 
n=0 n=0 
PiPa-++ Pn 
re Hpk o + n 7 
n=0 
PiP2 , PrPa --- Pe 
= UpyPo [1 +e+ a7 et 


he ae being independent of the service 


time distribution holds for the more general system M/G/o/c. 


1) 
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= ; APPLICATIONS é 
From (8) and (9), we obtain the proportion of lost calls or the probability Proceeding as in See. 1.6, we obtain the Howi ei is 
; d ? TOhowing ¥. a ier 0 - 
that incoming calls are lost as (1980)] ar ro a (Cohen 
- 
; ae Pa PiP9 «+: Pele ! (15) P C—p, (p | n| | 
: 0) oss) = e 3 * b Aa. 7 ‘ iy P2) “]M | 
. | i ; 1+, +(Prp2/2 !)-+...+(PrP2 +-- Pele !) °—PitpyEe(p,+-p,) No(p, +p.) aoe eee 
j ‘ a 
which is independent of p,. If we set pp=—p (i= 1,2,...,€), (10) yields = pi Bolp, Pa) 
the Erlang loss formula. VSP EMER MP fa) ia Phe 


E ; where 
Assuming calls to originate from a finite set of N sources, and choosing 


the mean service time as the unit of time so that ~—1, we can write (16) N.(p) = s (p'/i!); Ey(p) = (ole WINslo\: 0, an 

&% =(N—i)A. This yields the O’Dell distribution if A is constant and the ee CETL EOS Bathe, Se Vg 
Engset distribution if A depends on loss. Substituting for fi in (10), it is easy The steady-state exists when fee 
to verify that the denominator can be expressed in terms of the incomplete 


beta function ratio : Fhesprobability of a: call being lost is 
: Ss > = cH, ( | 
Ase (17) P( > ya a S __ il Py+- Po) 
J x(1—a)N-¢-* dx n=e C—Pi + pPHe(p,+-py) 
(11) (1-FAyw, 2aty : 
I rp ee ; yep oe : 
f ae(1—x)N-c-2 dy The delay (queueing) time density function faw) is given by 
0 
3 (18 (w)dw = z P e-caw . (cuw)n-e . ; 
Letting N +o, so that AN =X’, a constant. (and consequently A > 0), we ; sas “0 i (n—c) ! ie 


shall see that in this case (10) will yield the Erlang formula : ; : ee 
which simplifies to 


pe x \ \ ‘ ; 
(12) Pfloss} = [(A’)¢/e ] / | y (A’)ifi (19) Sa) as P( 0)-wle—py)e (e—p ;)uw 
i=0 ; 

The average delay for all calls from the first source is 

; (b) toae-Delay Systems. Let there be c channels servicing calls origi- (20) Wa = P( > O)/[(e—p,)) 
nating from two input sources with Poisson arrival rates Aj, Ag, respectively, h Anlaet 
The channels have identical exponential seryice times with parameter LL. Seer ee Ses Semel i 

We assume that the calls from first source can wait for service, but the second (21) Wa/P( > 0) = 1/[(e—py)ye} 


source calls are lost if all the channels are busy, hich is independent of p 
. whic e 0 2: 


Let P, be the steady-state probability that there are n units in the 


system, 7.e., n channels are busy if » ag 1 aes ee Fo avennicn 


and n—c calls are waiting if n > c. 


Then, using state-dependent parameters as in Section 1.6, we have Management of inventories is an area where the concepts of queueing 
(13) he? ot Bad U8 tc, ne theory have an obvious application. A possible model is that in which the 
: ee ae me units in inventory are regarded as forming a queue, and the arrival of a 
= A, aBe customer order is likened to a service completion, since such an arrival results in 
‘ aoa eae 2 he nee 4 
de the gq >. e ‘units’ may arrive in inventory singly or in 
(14) ‘Uy’ of Sec. 1.6 = ny, ees depletion of the queue. Th y 4 g 


certain lots of fixed or variable size. In the latter case we have a bulk- 
eked x, NZ C, arrival queueing system. The time between the arrivals of customer orders 


[ See. 5.4 
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The time of ‘serving’ the customer 
> 


: as » service time. 
may be regarded as the s¢ . 
nay from inventory may be considereg 


the time spent in issuing the item 


18s, Sax iy 
‘service time’’ defined above. 


negligible or included in the 
assumptions can be considered valid, a better corres. 
pondence between an inventory system and a ee oi “* obtained 
by regarding arrivals of customers or customer orders in oe inventory ay stem 
y reg 2 The reason ig that 


as the ‘arrivals’ of customers of the queueing alta 
The inventory replenishment 


Though the above 


both of these are more or less uncontrollable. | 
time or the lead time will then be considered as the service time of the queue- 


ine system, and both of these are, in general, controllable by the management 
of the system. We give below some interesting applications. 


(a) The Multi-Channel Queue Model. Let us keep a certain number ¢ 
of units in stock and order for replenishment each time an item is sold from 


stock. The time taken in effecting a sale from available stock is negligible, 
Let the customer orders for sale arrive at a mean rate A and let the average 
replenishment time be 1/z. Further, assume the inter-arrival time of customer 
orders to be exponential and the replenishment time distribution to be h- 
Erlang, and that any customer orders arriving when there are no items in 
stock are lost. The system may thus be regarded as a queueing system with 
e channels. The arrival of a customer order results in the depletion of the 
stock by one unit and the placing of a replenishment order. Thus we may 
regard the customer order as the arrival of a unit to the queueing system, 
an empty space in inventory as a busy channel, and the replenishment time 
as the service time. The inventory system under consideration may, there- 
fore, be regarded as an M/E;/c/c system with no queues formed. Note that 
in this system the number in the system is equal to the number of busy 
channels, which is the number of replenishment orders outstanding in the 


inventory system. So, the (steady-state) probability that n replenishment 


orders are outstanding is given by 
€ 
(1) P,, = (e"/n !) {| 2 (94/7 !), n = 0, 1, 2, ..., ¢ 
j=0 


where p = A/u. The mean number of outstanding replenishment orders is 


ée 
(2) L= & mP,, = pS(p, c—1)/9(p, c) 
n=l 
where 
c 
(3) S(p, 0) = & [ePpifj 1] 
j~0 
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The functions S(p, ¢) can be evaluated with the help of Poisson distribution 


tables. 
The mean number of sales per unit time is 
: c—1 c 
(4) Ly = A(L—Py) = AE (pill) | E (pif !) 
j=0 j=0 


— AS(p, c—1)/S(p, Cc) = pL 


Since an order is placed every time a sale is made from stock, the mean 
number of orders per unit time is equal to the mean number of sales per unit 


time, t.¢., L, or wL. Also average inventory is equal to the average number 
of idle channels = c—L. 

Let G be the gain or profit made per sale, i.e. the selling price minus the 
purchase price, Cp the cost of procurement or the ordering cost, and OC, the 
inventory carrying cost per unit time. Then, the average net profit per 
unit time is 
(5) P(c) = G-~L—Op-pL—O,(e—L) 

= (gu+C)L—Cye 
= (g4-+Ci)-pS(p, c—1)/S(p, ¢)—Cyc 
where g = G—Cp. If c* is the optimal value of the maximum stock c, we 
should have 
(6) 


Therefore, c* is given by 


P(c*) > P(et—1) and P(c*) > P(c*+1) 


(7) T(p, ¢*)—T(p, c*—1) > Cy{(gu+Cy)p} > Tip, c*+1)—T\p, c*) 
where 
(8) T(p, c) = S(p, c—1)/S(p, c) 


An easy approach to find c* is to look for the smallest c* which satisfies 


the inequality 
(9) 
which gives 


(10) 


P(ct+1) < P(e’) 


T(p, °+1)—T(p, ¢*) < Ci](gu+Ci)p 


It is worth noting that instead of considering the profit as we did in 
writing Eq. (5), we could have written the total cost equation as is usual in 
Thus, the average inventory carrymg cose per unit time 


inventory studies. 
is C, times the average number of items in stock, t.e. C(e—L) and the ordering 


Also, average number of lost sales per unit time is 


cost per unit time is Cpl. 


14 
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Therefore, the shortage cost per unit time, 


given 
Hence, the total relevant cost pe 


YT unit 


: | AP, = A—pL, from (4). 
by the loss of profit, is G(A—p#L). 
time is 
(11) 


where P(e) is given by (5). 


TC. = O(e—L) + Cp eb + G(A—KL) = GA—P(c) 


Since GA is constant (independent of c), it i aes - hia ey | of 
the profit P(c) and minimization of the cost 7’.C. are oe ent. 

(b) Re-order on Depletion of Stock. Next, we conmeley te model where 
we start with a stock of ¢ items and order for i ay when the inven. 
tory is reduced to zero, é.e. when all the units in stock are sold out. Let 
¢ Sai: orders arrive in a Poisson fashion at mean rate A and let the replenish- 


= en a pepe ng, 


ment time have a general distribution with mean I/v. The c units in inventory 
will last an average duration c/A and the average duration of the out-of-stock 
tiation “will hg l/u, the average replenishment time. So, in the steady 
- : 5 i, > n . F . ; } 
state the system will be out of stock for a fraction of time (1/m)/ [(¢/A)+-(1/p)). 
Sou a fa: = . a By : - ; “7. 
It will be convenient for this model to define P,, e be the probability that 
there are n units in stock. So the probability of being out of stock will be 
(12) P, = pi(e+p), where p = Alu 
For the remaining fraction of time c/(c-+p), the system behaves like a pure 
death process with the transition equations : 


‘ AD, SAP Tags) fea ls eee 
(13) } leg ay Neh 
| AP, = #P, 


From Eqs. (13), we have 
(14) P, =P, =... = P, = Wle+p) 


Therefore, the mean number of units in stock is 


ce 
(15) xX aP,, = e(e+1)/2(e+p) 
=1 


Mean number of sales per unit time is 


(16) A(1—P,) = Ac/(e+p) 


Also, with probability P, only one replenishment order is issued during 


one average replenishment time 1/. Therefore, the mean number of orders 
issued per unit time is 


(17) #Py = Al(c+p) 
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Therefore, the average profit per unit time as a function of ¢ is 
Ane) P(e) = [@Ac—OpA—Cye(e-+1)/2\/(e-+p) 
From (18) we ean find the optim 


mizing a function of a discrete 
For 


Regard c¢ 


al value of ¢ by using methoc 


ls for maxi- 
variable, 


[Cf. Eqs. (5) through (7)]. 
an approximate solution, 
a8 & continuous variable 
justified if c is considered | 
zero we see that the 


however, we can proceed as follows. 


and replace e(c+1) by c?. This will be 
Then, setting the derivative P'(c) 
optimal value c* should satisfy 


arge. equal to 


(c + p)(GA—Cie)—[GAc—Cy, - A—Cic2/2] = 0 


or, 


(19) ct 2pe—2A(Gp+Op)/C; ——) 4 || 
which yields 

(20) o* = VP T2NGp+O,)1C,—p 
It is interesting to note that if we take 


the replenishment time to be 
zero, so that p = 0, Eq. (20) yields 


(21) c* = V 2005/0, 
which is the same as the Wilson formula for the economic 
(E.0.Q.). This simplification of the above 
in the above approximate solution we are reg 


order quantity 
result could be expected, since 


arding ¢ as a continuous variable, 
the demand rate is uniform and no shortage is allowed. 


Also, in this case 
the term containing G in (18) 


, Which represents profit made on sales, will not 
contribute anything to the solution, since in this case no sales 


are lost. Thus 
the model discussed reduces to the usual E.0.Q. model. 


It may be noted that the policy of ordering for a single unit for each 
item sold will be suitable in case of bulky and expensive items, so that inven- 
tory holding cost is high compared to ordering cost, whereas the policy of 
ordering for all the ¢ items stored on (or close to) depletion of stock may be 
considered advisable where ordering and handling costs are high in comparison 
to inventory carrying and shortage costs. 


(c) Use of a Re-Order Level. We now consider an improvement of the 
model discussed in (b) above, where a replenishment order for Q items is sent 
whenever the inventory level comes down to D (the re-order level), so that the 
maximum inventory level for the system is Q+D= M, say. As before, 
customer orders are assumed to arrive in a Poisson fashion with mean rate 
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t times are taken to be exponentially distributed With 
ymen ee 
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A and the replenis! 


te equations for this system are 
The steady-state © 


mean I/#. 
APu — “PM_@ = uPp 

(Q<n< M) 

WP, = APaw (D<n <Q) 


9? - 5 
sii (A+2)Pr = AP ys iG er 


pePo = AP, 


AP, == fh P36 + -AP;, +1 


The above equations yield the following solutions : 

1e § é a 

(23) P, = W/AYATM/AY Fo (0<" <D) 
: = (w/A){(A + u)/A}PPo (D <n < Q) 


= (4 TA w)/APP— AFH) /AV OB, (Q<n< M) 
where 


» oper SC] ee ee 
(2 ag 2 ; A / if J 


Therefore, the mean number of units m stock is 


(25) I= S nP = AEM) PE(Q-+2D-+1)e— 2A} 2+ AP QAP + Qua] 


Mean number of sales per unit time is 
(26) S = A(I1—P,) = QAu(A-L pe)? /[AP + Qu(A+ fH)? ] 
And. as in Eq. (17), the mean number of replenishment orders issued 


per unit time is 


D 
(27) Raped P, = AwAtpP/PAPr + QHarey | 
n=0 
Hence, with the notation of Eq. (18), the expected profit per unit time is 
(28) P(Q, D) = GS—Cl—Cpk 


which is a function of two discrete variables Q and D, with S, J, and R given 
by (25)-(27). For an approximate solution for optimal values of Q and D, 
we can differentiate w.r.t. Q and D and equate the derivatives to zero. This 
gives two equations to determine optimal Q and D. 

However, an alternative method of finding the value of D, which includes 
a “buffer” stock, is to fix P,, the probability of being out of stock at a pre- 


assigned value p, which can be chosen by management and then obtain, 
from (24), 


(29) D = log{(A/Qu)(1/p—1)}/log(1+-p/A) 


. 
a 
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It may be noted that (29) expre 
so that the base of the logarithms 


<presses D as the ratio of two logarithms, 
18 Immaterial: 


natural logarithms or common 
enient, 


logarithms may be used as conv 


ey ys 7 rm ¢ - } 
For another approximate analysis. w 


‘fe Mav 1G ac wins 7 
small values of P Y proceed as follows. For 


, Wes ave l 

Then, in (25) (27), 1 s ca have {A/(A+p)}> = h, say, a small quantity 
, 5)- » 4, 8, and F can be expressed i j he 
; ssed in powers of h ] to the 
first power of h. These v; ia 
| I oe es These values can then be used in (28), thus simplifying 
the analysis considerably. As a first approximation, if we set h — 0 it can 


> a rifie 0 0 7) 5 
be verified that the equation aQ P(Q, D) = 0 yields 


i 
(30) Q —— / 2200/0, if 
¥ : F ; : 
showing that even for this complicated model the Wilson formula gives an 
approximate value for the optimal order quantity Q*. However 8 better 


approximation to Q* can be e 
a multiple of h. 


asily obtained containing « correction term as 


5.5 Applications to Machine Interference 


Machine interference problems are typical of textile industry, where 


a 
number of automatic machines, like thread-cutting or wire-dr 


awing machines, 
are put under the charge of a single operator. From time to time a machine 


may stop and require a corrective action by the operator, after which it again 
starts working. Another example of this type is that of workers, who from 
time to time have to stop their work in order to grind their tools at a grind- 
ing machine. If at any time more than one machine need the operator's 
attention or more than one workman report for grinding tools, a queue 
develops. One may like to find the optimal number of machines to be 
assigned to an operator or the optimal number of workmen assigned a grinding 
machine. 


The queueing model corresponding to the above situations, where the 
arrivals are generated by a finite source, has already been discussed in 
Section 1.6(f). Let & machines be under the supervision of a single operator, 
whose service time is exponential with mean 1/m, and let the stoppages of 
any machine occur in a Poisson fashion with mean rate A per unit time; then 
it was shown in 1.6(42)-(44) that the steady-state probability of finding n 
machines failed is given by 


(1) P,, = {k '/(k—n) YAl/n)"P, 


Yn, 
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where 
. &k , 
(2) P t/| E th i/(k—j) YAlp) |. 
= 0 ju 


From (1) and (2) we can write 


usa / k / I ! : ’ 
(3) P, al Ci eed (aks) One We 


5 eA fA)R-I(k—j) | 


jmO 

which shows that k—n, the number of machines running, has a truncated 

i istri i i : ster ys/ fy » mean number of machines 
Poisson distribution with parameter /A, ‘Thus the mean 1 P 1 


running is 


4 ( 
‘ ea quaye Hk 1p) 
(4) x (k—n)P, = (u/A) | : 
n=0 we € HAW] A)* 4/(k J) 
jad 


P(n/A; k—1) 
P(u/A; k) 


the last expression being in terms of the cumulative Poisson distribution 


(ue /A ) 


given by 
r . 
(5) P(m; 7) = & [e-™m4/j 1]. 
i) 
Dividing the expression obtained in (4) by k, we get the expected proportion 
of machines running, which can be maximized as a function of k to give the 
optimal value of k, when A and p are known. 
Next, consider the case of / machines under the supervision of ¢ operators. 
Then we have the model described in Problem 9 of Problem Set 1 (with NV 
replaced by k), i.e. 


(6) A, =(k—n)A, Og n<k 
16). n>k 

(7) Mn = NY, lin<e 
= CH, n>¢ 


We assume that c < kh. Then, proceeding as in Section 1.6, we have 


For 0g n<e 


(8) p= PAE IA... (h—n-+1)A Pp 
fb. 2ft 6. ML 


n 


(NG) 


and, forc<q nck 


KA.(—1)A ... (k—n-+-1)A _ k n | A \” 
( Mees es 
(9) = f.2p ...(€ —T)u.cu.(cy)n-e 9 c | ene Po 
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where 
(10) Po [= (IG eS i a 
nO \ WL fh ) nec ( n } Kn ile ie ¥é, G ) | 
= 1/F(k), say | 
As pe, the mean number of machines running is given by 
1 
(1) S (km, = bP —1y/ (Ly 
17 
c~] k =) | \ k- Pe! P 
(ula) | D (-1\(# Pg ee (TY, 5 ited 
4, ies j Holer\5) | 
C— Ie \ knJ k \ , 
TG are via eme 
BAGS) BA sletel 4) 
(ey! Bk ec ; emt} 4 SU ee IA; k 
= thei A+p 6 | eked P(cw/A; —c—1) 


A+p\k me 

( =) Alt gi \g nO 
, / ; cl ck-e 

where P(. ; .) is given by (5) and B denotes the cumul 

tion given by 


r LY 
(12) B(n, p;r) = = he ) pi —pyn-4 
q=0 ' wy] 
It can be seen that for ¢ = 1, 
Also, for the case c = k, (1 1) yields the 


These results follow also from the fact that for c = k, 
above becomes the brand-share model of Section 1.6(g) 


Further, from Kq. (8)-(10), we find that the 


channels or the mean number of machines being repaired is 


c—1 k 
(13) 2 nP,+ & cP, 


n=0 n=c 


ee eae ney ioe (aim 


n=e+1 B 


ls e-1 ,k—1 \n  k-1 | k—] | A\™ 
= (FPL EC, May +2, let (AY) 


n=c 
ka 
= Uy Si ue) a 
(7) Fea) TT 
fi . & - a \ 
< LIBRARY 


bbc 
Hy ow 

Sd par rt a® > uw 
HA OUT <a 


4 


PlcplA; k—e) 


ative binomial distribu- 


the final result in (11) simplifies to (4). 


mean number of machines running 
= ku/(A+), and hence the mean number of failed machines = kA/(A 


jt). 


the model discussed 


mean number of busy 
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Using (11) and (13), we have the - 
ge i - waiting, given Dy 
system, being repaired or W ri ers.” 
alae L = k—kF(k—1)/F (e) = KO — (PRD 
(14) 2 i. or repair, given by 
and the mean number waiting for repal’, & 


I,g=L —(kA/p) FE (k- 1)/F(k) 
a r Ee 
oe = A —(1-- APR D/F EH 
Also, the mean arrival rate to the system 1s 
a : I L 
. > ; / since 2 nb, = L 
(16) v= Ss h—n)AP, = (b= L)A, sine Ee ae. 
n=0 
Therefore, the mean time spent in system Is 
I W = Liv’ = Li[(k— L)A| 
(17) 
and mean time in queue 1s | 
(18) Wa =L,|v = “qi[(kA—L)A] 
It can be easily verified that 
: ee ear 
(19) W—Wa = l/h 


as expected. 


5.6 Applications to Scheduling 


| Sec. 5.6 


, number of machines in the service 


x j ‘A : é é é shine or ma shin 2 
2 scne i 4 ts Or rocess Ing vt a& mé Ac e 
[ he pr oblem of 5 sh eduling jobs fo p =] LC 3 


€ € j i j S 7x > © t | 7 si if roblem 
tr out f a& waiting line of jobs 18 ust an. example of 1 queuell £ Pp 
centre ut 0 TS 


f f ry 3 rhe era ¢ s ap licable to 
a id th re ore moO) st of the results 0 pre vious ¢ hapters are ay Pp I 3 
al ’ 1e7r “/s Sob ¢ e i 


be } + i S servic ar J } jobs and 
8 en f ju wwe service are the i] 

b ‘e the customers OI unis rec 

S¢ heduling. Here the , ) wing 5 re 


: > j at the shine (centre), 
the service operation is the processing of the jobs at the machine ( ) 


i I asi is snera 4 scheduling 
I the theory of scheduling, the main em yhas 1s 18 ge nerally on = | 
n ne F P : 4 o . . * : a their ef fee 
? ules (corresponding to queue disciplines nn que ue ng theory) and { 
| 5 i 5 shi 4 1e Means 
e u 7 Ne © f job ‘oug h De 1achnine OF nachines . low t n 
on the flow time 0 8 thr (24 th oD hi ~) I 28 | 1 } B, 


: Bats ; si ra 
the waiting (queueing) time plus the processing time which, in a queueing 
me eco ~ | =) 


system is called the total time in the system. 


Scheduling Rules 


Some of the scheduling rules used in practice Byyend upon the ee 
time, like the shortest processing time (SPT) rule according to wpich he P Ds 
are processed in ascending order of the processing nro eer 
As against this, there are three important ae which paid in $i Bae 
processing times: these are the first-come first-served (HORS), the last ai 
first-served (LCFS) and the random selection or service in random o 
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(SIRO) rules. 


Of these, the YOCRS discipline 
in. earlier chapt 


ers. We have 
usual notation the 


(1) 


has been discussed at length 
seen that in an M/M/1/¥ CRS 


system, with 
mean time in the queue, Wh 


@ 18 given by [Cf. 1.5(28)] 
Wa = Alu(u—A)} : PIM 1—p)} 
Proceeding as in 1,5(28) 


, We have the second 
about the orig 


moment of the 
in given by 


queueing time 


o 9 
(2) E(w?) = Jw. py pje-U-ewe diy = 


04 AC p)* 
Let us now consider the M/M/1 
when the server is 


duration of the 


/LCFS system. Suppose a job arrives 


this job’s queueing time 
residual processing time of the 
processing times of 


busy. Then, will consist of the 


job being processed plus the 


all jobs that arrive during 


of all subsequent arrivals unt 
arrivals, 
the 


this processing time and 
il the server first becomes free of 
Thus, referring to the 


queueing time of a job arriyv 
ordinary busy period whose de 


Now, since the queue length d 


all subsequent 
paragraph just after 2.3(11), we 
ing during a busy period is the 
nsity function has the 


see that 
same as an 
Laplace transform Oy. 
tion does not depend on the 
LCES discipline, 


istribu 
cipline, it follows that under the 


an idle period with probability | 
zero. Hence the queueing time d 


(3) 


quene dis- 
too, a job arrives during 
p, in which case his queueing time will be 
ensity function is given by 

fa() = (1—p)8(w)-+p.-\(a,) 

so that the Laplace transform of the 


(4) fords) = 1—p+-pa, 
where 6(-) is the Dirac delta function and a, is give 


(5) 


density function is 


n by (see 2.3(10)) 
ay = [9+ A+ w— (0+ AFM) — aap] (22) 
and the subscript Z signifies the LCFS discipline. 


From the above it follows that the mean 


queueing time under LOFS 
discipline is 


_ Df 2(A-+4 ft) 
7 = — ( = — |— 1. es = 
(6) Wat Jutl) 2a | V(A+p)2—4Au | 
anes i [1- mh } assuming > A 


= A/[p(u—A)] 


12 ac na 


g uelpul 3° 


- piepuels 


oun 
a asog 2Hand UY: f 


a 


man ew! 


& 


~ kenwer Lz 


uged Z 


agg Alet 


oe 
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, oricin of the queueing time is 
sec rent about the orig 

Also. the second mon 


= vy, Vie 2p 
(7) Ewin) = TeTRO) (u—ay = “ae 
Comparing (6), (7) with (1), (2), we see that 
(8) Wot = Wa 
1 aa, 
(9) B(wz) = -—,* Elvi) 


that is, the mean queueing time under LCES is the same as that uneaK FCES, 
but the second moment in the former case is greater ees _ in the latter 
by a factor 1/(1—p). It follows that the queueing time in LCFS has the 
same mean as in FCFS but a larger variance. 

For random selection for service, the reader may refer to Saaty (1961), 
p- 243 ef seq., or to Riordan (1962), pp. 103-106. It can be shown that for 


this system, M/M/1/SIRO, the first two moments are as follows : 


(10) War = Wo 
2 u ; 2 
(11) E(wir) = i pp » E(w?) 


where & stands for random service. Relations (10) and (11) hold for the 
multichannel system, too. 

Non-preemptive priorities were discussed in Section 4.5. For an analysis 
of various forms of preemptive priority scheduling rules, one may refer to 
Jaiswal (1968). 


5.7 Applications to Transportation 
Transportation systems, being basically service systems, are full of 
application areas for queueing theory. In what follows we consider two 


examples of such applications. 


(a) The Double-Ended Queue 

Consider the double-ended queue of customers (passengers) and servers 
(taxis) at a taxi-stand [Kashyap (1965a, b; 1966), Srivastava and Kashyap 
(1982)]. We assume that passengers arrive at the taxi stand in a Poisson 
fashion with mean rate A and form a queue if no taxis are available. Taxis 
also arrive at the stand in a Poisson stream with mean rate pf, and pick up 


one customer each from the queue, if any, or else form a queue. Thus at any 
time there is either a queue of passengers or one of taxis or of neither. Let 
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P,, denote tl 
i No steady-state probability th: iti i 
fe ady-state probability that n passengers are waiting in the 
jueue. We allow n to assume negative values, so that 
(i) when n> 0. it denotes the numbers of passengers waiting 
(1) when » — 0, neither passengers nor taxis are waiting 
g 


(iii) when n < 0, it signifies that —n taxis are waiting 


Ai suming « limited waiting space for N passengers or M taxis, we allow 
n to vary from —M to N. The steady-state transition equations are 


(1) = —AP_y+yuP_m,, 
(2) 0 = —(A+p)P,,+ AP -atHeP na (—-Mene< N) 
(3) 0. = —pPy+APy_, 


Solving Eq. (1) through (3) recursively as in Section 1.5, we have 
(4) P,, = pM*P_y (n = —M, —M+1,... N) 


where p= Aj. Now, using the normalizing condition 


N 
(5) PP eg 
—M 
we have 
. 1—p 
(6) P_y = 1— pMtNA »pF il 
so that (4) yields 
- 1— p)\pMtn 
(7) P= are »pAI,(—-M cn N) 


n ee pM +N+1 


Proceeding to limits as p— 1 we have 


1 


< MiNzi’ when p= 1,(—-M <n < N) 


(8) P 


The probability that the taxi-stand is empty is 


(1—p)pM : 
(9) Py a 1—pM+Nt1 ak a 


1 


ee 
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i sue of passengers (customers) is 
The probability that there is a queu¢ of p 8 ( 


N\pjM+1 
N P nee (l—p )p f p = 1 
1 


(10) Po = n = “| pM+NH 

N any 

= MEN+1' f 

The probability that there is a queue of taxis is 
=1 1—p™ 

og Po = 2 Fy | pment» P ~ 1 
M a 

~ RELI p 


f , * aagq ATC niet ars is 
The mean queue length of passengers (customers) 1 
nn 


N 
y Lo = XnP 
(12) Cc : 7 


pu +41—(N-+1)p¥-+NpN +H} 


eo ea 3 a | 
(l—p)(1—p¥+N+1) Ws 
_ NW pat 
~ 2(M@+N-+1)’ 
The mean queue length of taxis is 
—i 
(13) Ip = X (—n)P,, 
-M 


_ M—(M+1)p+pun 

~ (—p(i—puway P #1 
_ M(M+1) 

~ 2M+N+I)’ 

One may be interested in finding the probabilities that passengers or 

taxis are ‘lost’ to the system due to the limited waiting space. Clearly, we 


have 

(14) P{passengers are lost} = Py 

(15) P{taxis are lost} = P_y 

Also, 
(16) the mean number of passengers lost per unit time = APy 
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and 


(17) the mean number of taxis lost per unit time — LP_y 
where Py and P_y are given by (7), (8). 


Waiting Times 


Clearly, a customer does not have to wait if he finds 
(which occurs with probability Pp), 
present, he has to wait until ( 
function of the customer 


a queue of taxis 
Otherwise, if he finds n customers 
n-+-1) taxis have arrived. Hence. the density 
8’ waiting time is given by 


N , pntlyne—ue 
(18) we(x) = Ppd(z)+ © P, —— —— 

n=0 nt 
where d(x) is the Dirac delta function, and P,,, Pa are given by (7) and (11), 
respectively. Thus, the mean waiting time for customers is 


Le] 
(19) We= J xwe(x)dx 
0 
N (1—p)pM+n (par)n+1e—ue 
a H |—pMtNti al da 


N es 1) sje 
% (ppm nt 
and 1— pM+N* H 
where we have interchanged the order of summation and integration, which 


can be shown to be justified. Eq. (19) simplifies to 


M+1/]—(N-+1)pN-+ Nont1 
aay Wo a. EPL + ph 4 Wows} 


~ ap (1—p)(l—pMtNA) 


1 
=~ Le 
or, 
(21) Le =AWe 


which is Little’s equation. Proceeding similarly, we can obtain the density 
function of the taxis’ waiting time and show that 


(22) Ig = 2pWr 


where Wr is the mean waiting time for taxis, This is Little’s equation for 
taxis, 
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(b) Queues with Variable Capacity ; 
Let us consider queues of passengers waiting for buses at an Rican a: 
station like Kingston, Ontario when buses ply between. Montesil,/ uals 
and Toronto, Ontario. The intake capacity of the bus will be o—X, where 
4 is its maximum capacity (total number of paeceiger seats) ou X is random 
variable denoting the number of seats already occupied on eho at Kingston 
(after the passengers with Kingston laa feos disembarked), Let 
the customer and bus arrivals have Poisson Deena wank mean rates 
A and 4, respectively; so the system operates like an M/M&-*/1 model (Cf. 
Kambo and Chaudhry (1985)). 

The mathematical formulation and solution of the queue length process 
for the above model are given in Chapter 4, Problem 10, where it is shown 
that the steady state number of passengers queueing up for service has a 
geometric distribution. 

5.8 Other Applications 

It is a common experience to find queues or waiting lines developing 
in all spheres of human activity, so that an enumeration of the aieta of 
application of queueing theory can never be exhaustive. For specific applica- 
tion areas the reader may refer to more detailed application-oriented books, 
For example, one may refer to Lee (1966) for a very good account of applica- 
tions in air traffic, to Panico (1969) for applications in hospital administration, 
to Morse (1958) for applications to inventories and maintenance, to Prabhu 
(1965a) for applications to inventories and storage problems, to Hillier and 
Lieberman (1986) for applications to industrial queues, to Kleinrock (1976) 
and Allen (1978) for applications to computer science, to Newell (1982) for 
applications to transportation queues, and to various books on the subject 
like Saaty (1961), Cox and Smith (1961), and others for miscellaneous 

applications, 


Already in its eighties, queueing theory continues to remain a growing 
field. ‘‘May its tribe increase |” 


PROBLEM SET 5 


1. In a large community, people have been observed to go on foreign 

travel in a Poisson fashion and stay abroad for a duration having an approxi- 
mately exponential distribution, On an average, 1000 persons leave on 
any single day and the average duration of their stay abroad is 5 days. 
Regarding persons on foreign travel at any time as forming an M/M/co system, 
find the steady-state number of persons abroad at any time, 


[5000] 
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2. The telephone 


Poisson fashion at the 
of the duration of 


calls in a community appear to be generated in a 
average rate of 5 calls per minute, If the distribution 
@ call is observed to be 
mean of 2 minutes, 
provided in order to ensure a steady state 
the probability that a call is lost 
the average number 


approximately exponential with a 
umber of channels that should be 
- If 15 channels are provided, find 
due to all channels being busy. Also find 
of busy channels. 

{11 channels; 0:037; 9-63] 

3. Machines arrive at a repair shop in a Poisson fashion 
rate A. The machine downtime cost is 
the cost of providing 
Taking the service 


wis AL-VAC JO, 


at an average 
C, per machine per unit time and 
repair service at a rate pe is Cou per unit time. 


times to be exponential, show that the optimal value of 


[ Hine. The total cost of waiting and providing service per unit time is 


A 
O,L+C,p = 0, — +O, 
p—A 
4. A store dealing in costly items follows the inventory policy of carrying 
a few items in stock and ordering for replenishment as soon as an item is 
sold from stock. The cost of the item is 


$ 10,000 and the carrying cost is 
reckoned to be 20% per year. 


The management believe that the shortage 
cost is $ 2,000 per unit short per month. The demand is seen to be approxi- 
mately Poisson with an average rate of 2 units per month and the replenish- 
ment lead time can be considered to be exponential with a mean of one 
month. Considering the system as an M/M/c queueing system, find the 
optimal stock to be carried, 


[5 items] 


[If c items are carried in stock, the total expected cost for the lead time 
duration of 1 month is 


PO(c) = $10,000 x ¢ x (0-20/12)-+ $2,000 x Ly 


ae : e 2(2)e se 9i 9c ws —1 
= $500¢/3-+- $2,000 x “(c—)c— 22 [ = I SSP x | 
(using Eq. 1-6(24) for L,) 


Note that for existence of steady state we should have Aj/cw = 2/e <1, so 
that ¢ > 2. 
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Calculating 7'O(c) for a few values, we see that 
TO(3) = $2,277:77 
TC(4) = $1,013-83 
TC(5) = $ 912-60 
TO(6) = $1,018-02 
Thus optimal stock to carry = 5 items.| 


5. Assume the data as given in Problem 4 except that the management 
- ASS SS ~e ot: rl ‘ 
} 10 precise estimate of the shortage cost. Given that the company 
lave no precise es ssi te 5 i ie 
licy is to carry a stock of 6 items, compute the shortage cost, C, imputed’ 
policy is ary as 


7 the G any licy. 
ee eres Pes [$5,412 < CO < $23,487] 


[Write 7C(c) with $2,000 replaced by C. Now, we require that 
TC(6) < TC(5), and TC(6) < TC(7).] 


6. In a community of 100 individuals, persons switching into a brand 
and those switching out of it are observed to have approximately Poisson 
; is 1 . . 
distributions with mean rates of 30 and 20 per month. Considering the steady 


state of the switching behaviour, find 
(a) the mean and variance of the brand share 


(b) the probability of finding exactly half the community using the 
brand 


(c) the probability of extinction of the brand 


(d) the probability of the whole community being saturated with the 
brand. 


100 
. b 0-24)50 ~ 1-05 10-#2, 
[(a) 60; 24, (b) ( 2 \c ) 


(c) (0-4)! = 1-607 x 10-49, (d) (0-6)10 ~ 6-533 x 10-23] 


[Hint. Use the brand share model of Section 1.6.] 


*The cost that will make the company policy optimal. 


SEARLE TET aR 


APPENDIX A 
SIMULATION 


ory L. e 
A.l ‘To Simulate’ means ‘to make similar’. Thus, 


Sense, is syn 3 with y ; ildi 4 

aE YnNonymous with model-building. However, very often as we are 
doing here, the word simulation is used in its narrow sen 
precisely should be called Monte 


simulation, in its wider 


se to mean what more 
Carlo simulation. The words Monte Carlo, 
derived from the name of a place known for its casinos, have come to signify 
use of probability, and Monte Carlo simulation means the method of gerinest. 
ing values of a random variable from its probability distribution, as explained 
below. 


A.2 Monte Carlo Simulation 


Consider a random variable X with 


(P.D.F.) F(x), so that 


Probability Distribution Function 


(1) P(X < x) = Piz) 

and. 

(2) Bie = < t+ Ag) = (aw) Aa 
d 


where f(x) = a F(z) is the corresponding probability density function ( p.d.f.). 

Take a random fraction fe [0,1) by affixing a decimal point before a 
2 or 3-digit random number, e.g. a number 23 may be changed to the fraction 
R = -23 or 0-23. The fraction R is said to correspond to the value x of the 
variable X, given by 


(3) F(z) =R 
or, 
(4) “= F-1(R) 


where F- is the inverse function of F. 

The justification for the above method lies in the fact that the method 
implies that 
(5) Pe<X< a-+Ax) = F(x+Azx)—F(x) 
to the first order terms in Ax, which agrees with (2) and thus the value x is 
generated according to the given probability distribution. 

Thus, by taking a number of random numbers we can have random 
fractions R,, Ry, ... and generate as many values 2,, 2%, ... of X, 
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I method of simulation applies equally well] to the 
The above 


iables, the only difference being that in the latte 
and discrete variables, ; 


; pach value > variable 
function and, therefore, to each value of the variable a, 
is a step func 


86t o 


. : t abi random 
be associated. Thus, if the discrete variable bakes ii Bh: 
numbers can asso ; 
iti }= 1, 2,..., ”, then 
4, with probabilities pj, + = ae 
Re (Fea Py] > X = XY 
(6) 
i 
: = 
where P; = & pr. 
r=l1 
A.3 Need for Simulation 
pe . desert f TA an e me 
The method of simulation described above : p used to nota 
; Ber ue 1g 4 the real ac vantage 
problem concerning random variables, but al advs ees 


“mulation 
is obtained for solving problems where an analytical solution either dies 
Ss é 2 
not exist or is too involved to be practicable. Whereas an 


alytical Methods 
enable us to express the outcome of a decision as a function of the decision 


mal decision, for 


ear Programming, 
simulation helps us to work out the effect of each feasible decision and then 


to choose the best decision out of those tried. ation helps us in 


As mentioned jp 
any values of a random 
andom numbers. 
random vari 


variables and, therefore, help us in arriving at the opti 
example, through the use of calculus or a technique like lin 


Thus, simul 
arriving at a satisficing* rather than the optimal decision. 
Section A.2, we can use simulation to generate as m 
variable as we like by using sufficiently many r. Evidently, 
by generating a number of values of the able or variables, we 
obtain a sample of values of the variables in question from the population 
of values corresponding to the given probability distribution, As is true for 
8 with the size of the sample. 
“how large the simulated sample should be” 
and the answer lies in striking a balance be 


all sampling, the precision increase 


However, 
the question 


is an involved one, 


tween the cost of sampling and the 


cost of making a wrong decision. Besides, the stati 


stical principles pertaining 
to sampling theory, like variance-reduction techniques,** are applicable. 

Ex. 1. Simulate, taking 20 consecutive arrivals, the 
system with mean arrival rate of 10 per hour and mean 
minutes, and obtain the 


estimates for Ly, Wo, 
that the server ig idle, 


M/M/1 queueing 
service time of 4 
and the proportion of time 


For generating y 


alues of interarriy 
digit random numbey 


al times we take 20 ¢ 
"S, 


onsecutive two- 
Choosing the numbers 


from the 26th row of the 
(1981), Chap. 10, 3 esis in eo ae 
r (1963), p. 93. 


*Seo, for examplo, Grassmann 


**8oe, for example, Toche 
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‘Table of Random Digits’ in the book, we have the numbers 56. 49, 07, ..., 28 
as given in Column, 1 of Table 1 below. 


Let us, 


for example, take the 
we have the 


second number 49. Corr: sponding to this 
value of interarriva] time, 


“, given by 


l—e-** — 0.49 
whence 


« = —(In 0-61)/A = —6In 0-51 min, ~ 4:04 min. 


This value is shown in Column 2 of Table l 


Similarly for other values of 
interarrival time shown in the table, 


TABLE 1 
R.N.* Inter- Arrival R.N.* Service Service Service Waiting Idle 
Arrival Time Time Begins Ends Time Time 
Time (min.) (min.) (min.) 
(min.) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
56 4.93 4.93 92 10.10 4.93 15.03 ~ 4.93 
49 4,04 8.97 60 3.67 15.03 18.70 6.06 _— 
07 0.44 9.41 67 4.43 18.70 23.13 9.29 — 
56 4.93 14.34 66 4.32 23.13 27.45 8.79 — 
89 13,24 27.58 67 4.43 27.58 32.01 —_— 0.13 
70 7.22 34.80 87 8.16 34.80 42.96 — 2.79 
23 1.57 36.37 09 0.38 42.96 43.34 6.59 a 
35 2.58 38.95 65 4.20 43.34 47,54 4.39 — 
11 0.70 39.65 49 2.69 47.54 50.23 7.89 — 
20 1.34 40.99 31 1.48 50,23 51.71 9.24 _ 
32 2.31 43.30 11 0.47 51.71 52.18 8.41 — 
26 1.81 45.11 89 8.83 52.18 61.01 7.07 —~ 
98 23.47 68.58 67 4.43 68.58 a1 - 7.57 
01 0.06 68.64 49 2.69 73.01 75.70 4.37 — 
75.70 77.8) 4.83 ~- 
31 2.23 70,87 41 2.11 
2 10.29 81.16 19 0.84 81.16 82.00 — 8.35 
2 82.00 82.29 0.53 —_ 
81.47 07 0.29 2. 2.28 
i. pie : $8.02 —- 4.68 
60 5.50 86.97 23 1.05 86.97 8.02 a 
3.99 90.89 81 6.64 90.89 97.53 — 2.8 
+ 97.53 98.00 4.67 -- 
28 1.97 92.86 ll 0.47 5S 
82.13 26.32 


* Random Number, 
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A gal vd 20 two-digit random numbers for service time, Say from 
gain, we reé = © ‘ : “fez “4 9R 
| I r lumns of the random digit table, starting with the 26th row,* 
the last two columns i nay: i thes re ian 
1 obtain : 92, 60 11 as shown in Col. 4. From these we obtain the 
and obte : 92, potas Ba as Me Sen 
ice time values as shown in Col. 5. The remaining procedure ig gejg_ 
service PA: a Ne . syste : 
explanatory. Taking the time of the starting of the system as the time Zero, 


we work out the arrival times and the times that service begins and ends. 
. »* « é 


Since the server is idle for 26-32 minutes out of the total time duration 
\ Ss . 
of 98-00 minutes covered above, we have 


26-32 


Proportion of time the server is idle = 18-00 


O27 


. ror 29. 
Again, the 20 customers wait for a total of 82-13 minutes. Hence, 
Wa = 82-13/20 = 4-11 min. 


Now, the total 82-13 minutes of the column for waiting time can be 
considered to be in customer-minute units. So, by dividing this by the total] 


time of 98:00 minutes, we have the average number of customers Waiting 
(in queue), that is, 


2. 
Ig = 8213 ~ 0.84 customer 
98 


Note. For the problem discussed above, analytical results exist and are 
well-known. The problem is given just for the sake of illustration. 


The 
analytical values are : 


Sere e 10 
P{server is idle} = P, = 1—p = 1— ig = 0:33 


2 
ie p 4/9 


Ey a 1/3 = 133 


_ 


We = 1 jA = 7 x6 min. = 8 min. 
Ex. 2. Work out Ex. | taking the mean service time as 2 minutes. 


The calculation work is shown in Table 2, where the s 


ame sets of random 
numbers have bee 


n chosen for comparison’s sake, though we could have 


* The choice of the sa, 


ne me 26-th row for genearing arrival and service times is just 
coincidental. The two starting points can be tak 


en anywhere and do not have to he related 
to each other in any way. 


Sec. 


chosen different sets of 


simulation are seen to he 
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andom numbers. Now, the results obtained through 


feerver iaddigy .. oebe 
P{server is idle} — Q4ag ~ 0°62 
13-99 
Ma aug ~ OS 
13-99 
W,= = V- i 
q 20 0-70 min. 
In the present case, p = 10/30 = 1/3; so the analytical results are 
Py = 1—p = 2/3 = 0-67 i 
Pp Vicky a aes 
Ty fea 33 = 4 = 0-167 customer 


Wa = L4/A = (1/6) x 6 min. = 1 min. 


TABLE 2 | 
R.N.* Inter- Arrival R.N.* Service Service Service Waiting Idle 
Arrival Time Time Begins Ends Time Time 
Time (min.) (min.) (rin.) 
(min.) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
56 4.93 4.93 92 5.05 4.93 9.98 = 4.93 
49 4.04 8.97 60 1.83 9.98 11.81 1.01 — 
07 0.44 9.41 67 2.22 11.81 14.03 2.40 _ 
56 4.93 14.34 66 2.16 14.34 16.50 — 0.31 
89 13.24 27.58 67 2.22 27.58 29.80 — 11.08 
70 7.22 34.80 87 4.08 34.80 38.88 — 5,00 
23 1.57 36.37 09 0.19 38.88 39.07 2.51 — 
35 2.58 38.95 65 2.10 39.07 41.17 0,12 — 
1] 0.70 39.65 49 1.35 41,17 42.52 1,52 — 
20 1.34 40.99 31 0.74 42.52 43.26 1.53 = 
32 2.31 43,30 ll 0.23 43.30 43.53 oa 0.04 
26 1.81 45.11 89 4.41 45.11 49.52 — 1.58 
98 23.47 68.58 67 2,22 68.58 70.80 —_ 19.06 
01 0.06 68.64 49 1.35 70.80 72.15 2.16 — 
3] 2.23 70.87 41 1.06 72.15 73.21 1.28 — 
82 10,29 81.16 19 0.42 81.16 81.58 -- 7.95 
05 0.31 81.47 07 0.15 81.58 81.73 0.11 ~ 
60 5.50 86.97 23 0.52 86.97 87.49 — 5,24 
48 3.92 90.89 81 3.32 90.89 94,21 — 3.40 
28 1.97 92.86 il 0,23 94.21 94.44 1.35 — 
13.99 58.59 


ep 
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We observe that the simulation results are quite close to the 
e opse < F 
in Ex. 2, while they varied considerably in Ex. 1. Generally 
ones 1n ILxX. 2, : id 


j Spe 
we can expect better match for small values of p. We know 


that for 
to 1 we Should 
a long time to elapse before the steady state is obtained. For bette 
simulation 
simulation 
ast 10 of these, 


“ating, 
P>1, 
©Xpect 
r resy Its, 
stimates 


the steady state does not exist; consequently for p close 


it is advisable that several initial values are discarded and 
are based on the last few. For example, one may do the 


eXerciga 
for 30 arrivals and base the calculations on, say, the ] 
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PROBABILITY AND TRANSFORMS 


B.1 Expectations and Moments 


Let X be a continuous random variable (C.V.) with ( cumulative) Probability 
Distribution Fune 


tion (P.D.F.) F(z), so that* 


(1) Pix < x} = F(x) 
The corresponding probability density function (p-d.f.) is f(x), given by 
(2) 


Pee Xe x-+-da} = fi 


x\dx 
from caleulus that F(x 


It is easy to see ) and f(z) are related as follows : 
; d 

3 a) =" p 
(3) fw) = Fa) 


F(z) = f f(u)du 
2 
We define the expectation, or exp 


(5) 


ected value, of a function, (X), of X to be 
E{g(X)] af g(x)d F(x) = ! Go) f (x)dar 


the integrals being taken over the range R of X. 


by us are non-negative 
of the two equiv 


Most of the r.v.’s considered 
and extend over the interval [0, oo). 
alent forms given by the 
(5)"*, we shall usually prefe 


Moreover, out 
second and the third members of 
r the latter, so that we can write 


(6) Blg(X)] = f ga)flx)det 


Let a be any real number (which is usually 


aken to belong to the range 
R of X) and r an integer > 1; then 


(7) B{(X—a)r] = f (—a)*f(a)de 


**P(A)’ stands for ‘ probability that the event A occurs’. 


** Tho second momber expresses expectation as a Stieltjes or Riemann-Stieltjes integral; see 
for oxample, Chaudhry and Templeton (1983), pp. 345-348. 


| Tho Steiltjes integral form is more general, since it includes both the casos when X is & 


continuous or a discrete r.v. (in the latter case F(a) is a step function and so, at some points, f(x) 
cannot be defined). Still, for the sake of simplicity we use the other (Riemann integral) form, 
with which the reader is expected to be more familiar, 
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is called the r-th moment, about the value a, of X (or of the distribution of 
X). If a= 0, we get the r-th moment about the origin, usually denoted by 
#,; thus 


(8) 4, = B[X)] = i x" f(a)dx 


is the r-th moment of X about the origin, In particular, the first moment 
of X about the origin : 


ie) 
(9) fy = E[X) = J xf (a)dx 


is called the mean of X and is often denoted by . From the foregoing it is 
clear that the mean is the expected value of the r.v, X, and corresponds to the 
notion of an ordinary (arithmetic) average. Sometimes, though less techni- 
cally, the words mean and average are used synonymously. Moments about 
the mean, %, play an important part in statistics and probability. The r-th 
moment about the mean is usually denoted by y,; thus 


(10) pr = BYX— By] = fF (eB) fla)de 


The second moment about the mean, /4z, is called the variance of X. The 
variance, and its positive square-root known as the standard deviation (S.D.) 
and denoted by cy, are measures of dispersion or spread of the values of the 
r.v. X around the central value denoted by the mean %. The variance is 


given by 

(11) Var[X] = nw, = E[(X—Z)?| = ox 
It is easy to show that 

(12) My = fa—fy? = E[X*]—# 


Further, it can be seen from the definition ( 11) that the variance, and hence 
the standard deviation, will vanish iff (if and only if) we have a deterministic 
case when the variable can have only one value, say, 


(13) fle) = (wz) 
where 6(:) is the Dirac delta function. 


The mean and the variance (or 8.D.) together provide information which 
is usually sufficient for most of application studies. The mean alone gives a 
central value and, without the knowledge of variability around this value, 
often provides too little information. Moments higher than the second 
provide still more information about the distribution and shape of the p.d.f. 
curve. For example, skewness, a measure of asymmetry, is given by 


oo  _ 


Sec. B.2 | 
(14) 


skewnegs* — |H8? = p,|03 


and kurtosis, a measure of peakedness of the 


~ ito r4 
kurtosis ~ Mal p = plo 


Discrete Variables 
mL. * : . 

The foregoing discussion pertains primarily 

i slic ‘ "Tn i < 

With slight change of notation there are 

for discrete random variables. I 


et N be a discrete rv. 
0), 2)... with 


(16) P{N = 


=} = p,, n = 0,1, 25, awe 
Then, the mean, or the expected value, of NV igs 


(17) E{(N] = D> Np, = N, say 
n=0 


The r-th moment about the mean is 

(18) My = E((N—H)] = F (n—ay 
n= 

and the r-th moment about the origin is 


(19) 
And, in particular, the variance is 


(20) My = Var[N] = E[(N—ay] = 


B.2 The Laplace Transform 
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p.c.f. curve, is 
9 a continuous variable, 


corresponding results and theory 


, taking on values 


x 


rr a. Rin- x 


Given any function f(t), of a real variable t, the Laplace transform (L.T.) 


of f(t) is defined as 


(1) LU fl] = f e-stpeyat 


where s is a complex variable such that, for the sake of convergence. Re(s) > Ot, 
where Re(s) signifies ‘the real part of s’. Thus £[f(t)] is a function of s and 


* Alternatively, some authors define skewness — ftg?} 428. 


T In most cases discussed by us, the integral converges for Re(s)> Ox (Cf, (6)-(7), following.) 


wo 
17 ye Tie, 


£ 


‘ > 2 
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‘3 usually denoted by f(s), f*(s), and sometimes by corresponding lower case 
is usual vs = : : 


letters, as 


& [F()] = f(s) 


From (1) and B.1(6) it is clear that if f(z) be the p.d.f. of a random 


variable X ¢[0, co), we have 
(2) L[f(t)] = Ele**] 
Sometimes, the Laplace-Stieltjes transforms (L.8.T.) are used. The 
L.S.T. of a function, say g(t), may be defined as 
(3) LS{(g(t)] = ‘i e“dg(t), Re(s) > 
There seems to be no standard notation for L.8.T., though some authors use 
(**) as 
(4) LS[9(t)] = 9**(s) 


Very often we use the L.T. of a probability density function (p.d.f.) or 
the L.S.T. of a Probability Distribution ae Let F(-) be the P.D.F. 
of a random variable X, defined by B.1(1), and f(:) the corresponding p.df,, 
defined by B.1(2); then by definition it is clear that » 


(5) L[f(t)] = CS[F()] 


and we shall denote both of them by f(s). Thus, 


(6) jis) = fer Mf(bdt = [ ed F(t), Re(s) > 0 
0 


It may be noted that 


(7) f0) = 
since the left-hand member of (7) stands for the total probability of the r.v. X. 


Some Properties of the Laplace Transform 

Let f(t) be any function, not necessarily a p.d.f. (though in most applica- 
tions discussed by us, it is a p.d.f.) and f(s) its L.T.; then, since the L.T. is 
defined as an integral, all properties of integrals become applicable to it. 
In particular, we have 


(8) £ [of(t)] = cL[f(t)] = ef(s) 


where c is constant, 


(9) Li flt)+9(t)] = LEf)]+ L[9(é)] = fls)+9(s) 
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using the definition of the L.T 


Also, 


. and integrating by 8 sc 
show that : OS ee ae 


a d@ a 
(10) “| Gf0| = £7) = sfe)—fo) 


provided lim [e-stf(t)] = 0. 


mas [This condition is satisfied in most applications.] 
Similarly, using the definition and changi 


the resulting double integral, we have 


f t 
(11) £| Sflujdu | = flsy/s 


From (11), it follows that if F(:)isa PDF 


-, we have 
(12) LF (t) = fi(s)/s = LI f(t)\/s 
Also, if ais a constant, we have 

(13) Let f(t)] = fls+a) 


which is evident from definition. 


And, from (6) on differentiating n times w.r.t. s, we have 


(14) Life) =(—1)» & (Fay, | 


den n & positive integer) 


An important Laplace transform very often used is 


co 
(15) L[t"] = ff erstendy 
0 
l co 
= = | e“u"du, on using the substitution ¢ = “ 
8 0 8 
T(n+1) : nies . 
== a bing the definition of the Gamma function 
n! 


==> if n is a positive integer 
8 E 


From (14), on setting s = 0 on both the sides, we have the important result 
that for a random variable with p.df. f(-), the n-th moment about the origin 


is given by 


(16) = (—1)"fim(0) 


ning the order of integration in 


i 
i 
| 


ALLS 
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where / (0) signifies the n-th derivative of f(s) ats=0. Asa special case 
of (16) the mean is given by 
(17) # = —f'(0) 
B.3 Convolutions (Continuous Variables) 

Let f(t) and g(t) be two functions of a (non-negative) varibable ¢; then 
the convolution of f(t) and g(t), denoted by f(t) * g(t), is defined as follows : 


t 
(1) S(t) + g(t) = J flu)g(t—u)du 


t 
= J f(t—u)g(u)du, as is evident on substitution 
0 


= g(t) + fit). 
The convolution of three functions fi(t), ¢ = 1, 2, 3, is defined as 
(2) Ai) * fa(0) + f,(t) = Aid) (fo(t) « fg(t)) 


From the definition, it is clear that convolution (*) regarded as a binary 
operation on the set of functions of ¢. possesses both the commutative and 
associative properties, so that in the right-hand member of (2) the functions 
Silt), i = 1, 2, 3, can be permuted in any of the 6 different ways. 


The definition (2) can obviously be extended to the case of 7 functions 


Ait), ¢=1,..., n, for any integer n. Consider the case when 

(3) A =f) =... =f.) =f) 

Then, the convolution of the functions f,(t), i = 1, ...,, becomes 
(4) S(t) « f(t) «... (repeated n times) = f(t)*(™) 


which is called the n-fold convolution of f(t) with itself. In particular, the 
(2-fold) convolution of f(t) with itself is, from (13), 


t 
(5) fit) «f) = f fluyft—n)du 
0 
the epithet ‘2-fold,’ in fact, being redundant. 


The Convolution Theorem of Laplace Transforms 

Let h(t) be the convolution of two functions, f(t) and g(t), of t, which for 
our application purposes may be taken to be a non-negative real variable 
in [0, 00); then the L.T. of h(t) equals the product of the L.T.’s of f(t) and g(t), 
that is, with the notation introduced above, 


(6) LL) + 9) = LIU f)1.CIg(t)] 
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or, written alternatively, 


(7) h(s) = Fle)g(s) . 
The proof follows from the definitions of LT. and 


convolution on chang- 
Ing the order of integration in the resulting : 


double integral. 


(6) has obvious generalization to the case of n functions (for any positive 
integer n). In particular, we have ; 


(8) LL f(E)"™)] = (LE fty}n — (fle))” 


Convolution and Probability 

Let X, Y be two independent (non-negative) random variables having 
probability density functions (p.d.f.) f(x) and gly) and let Z = Xi y. Then, 
using the multiplication theorem of probabilities for independent r.v.’s, we 
see that the density function for Z, say h(z), is given by 


z 
(9) he) = § fe)ge—zx)dx = fe) s g(z) 
0 


and hence, using (6), the L.T. of h(z) is given by 
(10) h(s) = fis)gis) 
where f(-), g(‘) are the L.T.’s of f-), 9), respectively. 


Like (6), (10) is also capable of generalization to the case of any number 
of independent random variables. In particular, if X,, X,,..., X,, be n 
i.id.r.v. (independent identically distributed random variables) having the 
same p.d.f., f(x), the L.T. of the p-d.f. h(z), say, of their sum 


(11) Z= S X, 

is given by 

(12) h(s) = (fis) 
where 

(13) fis) = ‘f e*2f(x)de 


is the L.T. of f(z). 
(12) is very often used in queueing theory; e.g. Cf. 1.4(19), 2.3(11). For 
an example of the use of (10) one may refer to 3.3(34). 
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It may be noted that (10) and (12) will hold if we replace the p.d.f. by 
P.D.F. (probability distribution function) and the L.T. by L.S.T. and use 
B.2(5). Note further that the analogue of (9) in terms of the distribution 
functions will be 


z 
(14) H(z) = { Oe—a)d F(x) = F@) * Ge) 


where F(-), G(-) and H(-) denote P.D.F. of X, Y, Z, respectively, and I(s), 


etc, in (10), (12) are defined as in B,2(6). 


B.4 Convolutions (Discrete Variables) 
Let p,, J, (n = 0,1, ...) be two functions of a discrete variable; then 
n 
(1) Pn*In = % PmIn-m 
m=0 
is called the convolution of p, and q,. 
As in the case of functions of continuous variables, the above definition 
admits of similar extensions to the definition of convolutions of more than 
two functions and convolution of a function with itself. (See B.3 (2), (4).) 


Generating Functions 


Given a function p, (n = 0,1,...), the function 
(2) P() = 2 P»2” 
0 


is called the generating function’ of p,. The variable z may be regarded as 
a complex variable with a condition imposed on z to ensure convergence of 
Pz). 

If p, be the probability mass function (p.m.f.) of a discrete random 
variable N, say, i.e. 
(3) P(N = n) = Dy 
we shall require |z| < 1 for convergence. Also, in this case P(z), being the 
generating function of the p.m.f. of WN, is called the p.g.f. of N. 


Using the discrete analogue of B.1(6), we see that 
(4) Piz) = Bz") 


*The generating function is often called the Z-transform, and may be taken as tho discrote 
analogue of the Laplace transform, However, some authors define the Z-transform of Pn 8 


J 
2s Dnt, 
n=0 


COTES 


sap meee 
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Let p,, q, be the probability 


mass functio ) ] 
Fa ig ca awe! ns (p.m.f.) of two independent 
~ bb 2» 


respectively, i.e., 

* 

(5) P(N, = n) = P,(n = 0,1,...) 
(6) P(N, = n) = q, (n = 0,1, ...). 


Then, using the multiplicat 


ion theore  mrohahilite. <: 
of the sum N,+N, is rem Of probability we see that the p-m.f. ; 


re 


(7) PINALN, = 9) 


n 
» he P seg ae 
~ Pm 4In—m = Pn* In = Ty, BAY 


— 


Also, the generati i i 
g ing function of Tr, 1.€., the p.gf. of N,+N, is 


(8) R= Eres fF 
m=t) n=0 moo siti egal 
= st 2 PmIn—me (on changing the order of summation) 
Z Pm” ZL dy_mzr—m 
m=0 n=m 
= P@)Qe@) 


(8) is the discrete analogue of B.3(10) and can likewise he extended to the sum 
of any number of discrete random variables. In particular, if N,, N Ny 
es . 5 Mig Mea sak 
ne iid. r.v.’s having the same p-gf. P(), the p.gf. R(z) of N,4+N,+4+...4N, 

will be given by | 


(9) R@) = (P@)yr 
Some Properties of Z-T ransforms 


Let us write the Z-transform (or generating function) of p, as Z[p,], ice., 


(10) Zip.) = 3 pz 


0 


Thus, Z[p,] = P(), as given by (2). Then, the following relations are evident 
from definition 


(11) ZAlep,| = cZ[p,] 
(12) AP +9n) = ZPal+ Zan) 
And, (8) can be written as 

Z[p,, * Yn] — Z[PrlZ[dn] 


il 


— ee 2 mm~ x 


) 
) 
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Also, we have : : tia : 
= 2 Pr 4re” 
( 14) Z[Pn+r] “a a : 
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a > Pay" 
n=0 


Cl Rouche’s Theorem 


A function Sz) of a complex variable z 

ane aed i act ‘ if it is Single-valued and differentiable 
moO If two functions f@) 

C on the complex z pl 

the same number of 


is said to be analytic in a region R 
at every point of R. 

and g(z) are 
ane and | g(z)| 
zeros inside QC, 


analytic inside and on a closed contour 


< |f(z)| on G, then f(z) and f(z) + g(z) have 


C.2 Lagrange’s Differential Equation 5 
The general solution of the partial differential equation 


¢ | 
| (1) ed: Q ms =R 
where P, Q, R are functions of x, y, z, is 
(2) d(u, v) = 0 
or, equivalently, 
(3) » = flu) 
where ¢, f are arbitrary functions, determined from 


initial conditions. and 
u=da,v=b,aandb being 


arbitrary constants, 
of the corresponding 


are two independent solutions 
simultane 


ous ordinary differential equations 
dz dy dz 

4 — = 

(4) = ge 


PE —————— 


(4) are known as the Lagrange system of equations, 


C.3 Difference Equations 


Just as a consideration of increments in the argument(s) of a function 
of continuous variable(s) gives rise to differential equations, a simil 
deration in the case of a function of discrete variable(s) yields equations 
known as difference equations. An example of a homogeneous linear difference 
equation of order n with constant coefficients is 


ar consi- 


(5) WY rt nt WY etn—1+ 1s TOnYe =v 
Where dy, d, dg, ..., @, are constants. 


In order to solve (5), we use an operator B, defined by 
18 
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9 
(6) Bry, = Yor t= 9, 1,2, .-- 


7 , an . ; , 5 
so that operation with # results in enhancing the argument x by unity. Using 
EB, (5) can be written as 


(7) apH"y, +a,E"y,+.---+onY2 = 9 


or, separating the operator, 

(8) (a,H®+a,H"+...+an)¥, = 9 

The equation 

(9) a,H"+a,H"1-+...+a, = 0 

regarded as an equation (of degree n) in £ is called the characteristic equation, 


or determinantal equation, of (5). Let (9) have the n distinct roots 


E = My, Mo, «1+, My 
Then, the general solution of (5) is 
(10) Yo = GyM{+Coms+... +l, My 


where ¢, c, are arbitrary constants, to be determined from initial 
Piatt Be 


conditions. 

If the roots of (9) are not all distinct, let m, = m,, say, the other roots 
being all distinct and not equal to m,; then the first two terms of (10), corres- 
ponding to the equal roots m,, m,, will be replaced by (c,+c,x)m{, This 
result has an obvious generalization to a root repeated any number of times. 


The reader will have observed the analogy between the above and the 
corresponding theory regarding differential equations. 


For a more detailed but introductory reading on difference equations, 
the reader may refer to Yamane (1968). 


C.4 Differential-Difference Equations 

A consideration of variations in a continuous variable often results in a 
differential equation and, in the same way, a consideration of variations in a 
discrete variable leads to a difference equation. In queueing theory, when 
we deal with functions like P,,(t), where n is a discrete variable and t a conti- 
nuous one, and consider changes in both nm and #, the result is a differential- 
difference equation (sometimes also called a difference-differential equation). 
Examples of this are 1.5(4),(6). By using the Z-transform or generating 
function, the ‘difference’ aspect of the equation is done away with, thus 
yielding a ‘pure’ differential equation, see, for example, 2,2(4). However, 
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if we use Laplace transforms, the ‘differential’ 


left with a diffe ss beet is G00 rid of andwe are 


Asa ea ee a, for example, in Srivastava and K 
. ae : AK ) me plait: differential or difference equations may 

ved by usual methods*. However, usually both the Z and the Laplace 
mranstorms Gre used thus reducing the equation to a simple algebraic salon 
as 1s done in 2.2(8), 2.3(7). It may be noted that consideration of the sinaks 
state of the system makes the derivatives w.r.t. the time ¢ vanish tin dot 
away with the ‘differential’ aspect of the equation, as in 1.5(7) (8 | As ; 
further illustration. refer to Section 3.6, whe eek 
state of the system from the 


equation in P, (2), given 


ashyap 


re, since we consider the steady 
very beginning, we obtain a differential-difference 
ek by 3-6(9), (10). Otherwise, consideration of the 

re neve resulted in what may be called a partial-differential- 
difference equation, in P,,(a, t) as, for example, in Srivastava and Kashyap 
(1982), Eq. 4.4(3)-(5). For more details on differential-difference equations 
one may refer to Bellman and Cooke (1963) | 


Thus, the use of the Laplace and/or the Z-transforms results in a consi- 
derable simplification of the queue equations, but as is true for all ‘benefits’, 
this simplification, too, can be had at a certain ‘cost’. The cost in this case 
is that we obtain solutions in terms of transforms, which should be inverted 
to arrive at a transform-free, or explicit, solution. This necessitates considera- 
tion of the problem of inversion of the transforms. 


C.5 Inversion of Laplace Transforms 
If, as heretofore, 
(1) L[f()] = fis) 
we say that f(t) is the inverse Laplace transform or Laplace-inverse of I), 
and write 
(2) £Tfls)] = f(t) 
A simple way of finding inverse Laplace transforms is by using the results 


for Laplace transforms in the reverse order, Thus since. for example, we 
have 


(3) L[e-a#] = 1/(s+-a) 
we can write 
(4) 11 /(s+-a)] = e-at 


*The method of solution of a differential equation is simple and is covered in 
undergraduate calculus. Tho method for a difference equation is briefly discussed in Section 0.3, 


[ See, 0.6 4 
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The general method of finding the inverae Laplace transform is through Pir | 
. lace transforms which states 3 ifferentiation of a D fini 
the use of the complex inversion formula of Laplace trans r ‘tates ete 
; = C ider ay 
that if f(t) and f(a) ar¢ related as in (1), we have wnsider the definite integra} 
(1) i ‘. 
p(x) f fix, 2\dz 

a 


tite 7 1d 

hi i] oh f(a)d4, | Bes I 
(5) I ont ¢ o 

‘(t) here the limi ‘ 

(6) fij=0, t<0 where the limits of integration a,b are functior ‘ 

i ; Z mis of Zz, 

where the integration in (5 c in the complex It can be shown by using the d« finiti ; 

: > ; c : ‘ ; 

iy), the real number ¢ ition of the derivative that 

¢ lies to the right of all singularities* of f(s). 


the complex inversion formula, one may 
In the special c: 

yf al case when a. bh are o ; 

are constants, the last two terms in (2) will vanish 


) is performed along the line 
being chosen arbitrarily subject to the 

@) “He) f Llw.2) 4, db da , 

ae j qzt@: 1) 


z-plane (z w+ 
condition that z 

For more details regarding 
refer to Kleinrock (1975) or to hooks on 


(1946) and Doetsch (1974). 


Laplace transforme like Widder 


C46 Inversion of Z-Transforms 
This, too, can be achieved by using the Z-transform results in the reverse 
a Z-transform 


order as in the case of the L.T, However, the inversion of 


(1) P,] = & Pam 
Ned 


when the function Z[-] in (1) is given. 


implies that we have to evaluate P,, 
Thus the problem is simply that of writing down the coefficient of 2” in a fune- 


tion of z. for which several methods are known to the average reader. The 
basic problem is that of expansion of Z{-| in powers of z. In many cases, 
use of known expansions like the hinomial theorem is found useful; see, for 
example, Eq. 1.7(10), (11), Sometimes, Z[-] may have to be resolved into 
partial fractions before the use of the binomial theorem. Also, in a com- 
plicated equation, inversion may be effected by equating coefficients of like 
powers of z on both sides; for an example see Problem Set 4, Problem 5. 
For # discussion of the general method of inversion of Z[-] via the use 


of complex variable theory, the reader may refer, for example, to Kleinrock 


(1975), and Chaudhry and Templeton (1983). ; 
For an interesting example of inversion of generating functions, one F 

; . ; f 

may refer to Problem Set 4, Problem 12; one may also refer to Chaudhry 
i 


et al. (1987) and Briere (1987) for some more complex examples of inversion 


of generating functions using roots of characteristic equations in connection 


with bulk queues. 


A singular point, of a function of a complox variable, is a point whore the function fails 
to be analytic, ‘ Singularitios’ include singular points, and essential singular pointe. (Cf. Churchill 


(1972), 


om Tae 


APPENDIX D 


LIST OF SYMBOLS 


The following is the list of symbols and notation most commonly used 
in the text. These and other symbols and notation are, however, also 


explained in the text when they are introduced. 
D.1 Notation for Queueing Systems 
A queueing system is generally denoted by 
AX/BY/c/N 


where A and B, respectively, stand for the type of inter-arrival mane and 

service time distributions, c is the number of channels (servers) and JV is aie 

maximum number of customers admissible to the system. The superscripts 

X and JY, respectively, indicate that a variable number of eM D may be 

arriving or may be served together in a group. When N is missing in the 

| above notation, it means that N is infinite. A and B may be replaced by the 
following symbols with their meanings shown against each. 


D constant inter-arrival or service times 
Ey, k-Erlang distribution, i.e. the k-fold convolution of the exponen- 


tial distribution 


G general (arbitrary) distribution, commonly used for service 
times (replacing B) 


GI general (arbitrary) independent input distribution, used for 
° - ; . pee , oy 
inter-arrival times (replacing A) in preference over ‘G 


HE hyper-exponential distribution 
M (negative) exponential distribution 


Replacement of X by R indicates that the arriving group size is constant and 
equal to R. 


D.2 Symbols from Roman Alphabet 


A(u) P.D.F. of inter-arrival times 
A(z) p.g-f. of arrival-group size (= © a,2’) 
r 


i 


Te ce a ST TEIN RS ROSE NR oe 
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a 


ay 


a(u), as) 
Biv) 
b(v), 5(s) 
EH 
E[-], E(-) 
exp(x) 
FCFS 
FIFO 

Fw), f(s) 
Saw), Fa(s) 
I,(-) 
Lid. 
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mean arrival-group size (= ¥ Tar) 

r 
probability (mags) 
P(X = r) 


function of arrival-group size Pe 8 


"> 


p.d.f, of inter-arrival times and its L.T. 
P.D.F. of service times 


P-d.f. of service times and its L.T. 
finite difference Ope 


trator as in Hr Ee Boos 
expect 


ation, or expected value or mean value, of 


e” where e = 2.71898 --. 1s the base of natural logarithms 
first-come, first-served 


first-in, first-out 


p.d.f. of system time (waiting time in the system) and its L.T. 
p.d.f. of queueing time (waiting time in the queue) and its L.T, 


modified Bessel function of the first kind 
independently and identic 


generating function of ky 


and order py 
ally distributed 


probability of r arrivals (se 


rvice completions) during a service 
time (an inter-arriva] time 


) of arbitrary duration 
mean number of customers in system 
mean number of customers in queue 
mean number of custome 
Laplace transform of 
inverse Laplace transform of 
last-come, first-served 


last-in, first-out 


Ts in queue just after a departure 


Laplace-Stieltjes transform 
Laplace transform 
number of customers in system 


number in system just after the n-th departure ( 
the n-th arrival) 


number in queue 
number being served 


number in system at time ¢ 


just before 


(4 
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: ‘he order of x as % — @, SO that $04 
ci ee ee | ; times spent, in syst MIST OF SYMBOLS 145 
2 Stem ; . 
| O(a) | < Ka z stieretiee m and in queue 
. eMng Tuncti ee 
for x close to a, K being a eonstant |z| ncaa 10 Pparametey 
8 of z 
o(z) tterms negligible in comparison with « as v7 a, so that / Ba” Ge: dipkaie | 
; may! bel and Misee : 
é o() T(-) : Wwe ellaneougs Symbols 
lim fete ti) gamma function iT 
z—a aC efined as 
P(A) probability of the event A aye (pace | 
A 7 7 aes = : 2 
P(A |B) (conditional) probability of the event A, given that the event - | 
ag re. Y = (n—})! AP in ' 
B has occurred Omn Kronecker delta en ee yi nis an integer 
Pp. (steady-state) probability that there are n customers in system 1 (x) Dirac delta fy i Tia Ih = 96, ceno: Gtheewiga 
eo 7 pax : x . F é ; = a NnCti0on (— 0 ax se 
tap (steady-state) probability of finding 7 in system just after a i (= 0 except when x = 0, such that 
E b 
= F 
departur | J 9(x)de i bap i 
P,,(é) probability that there are n in system at time ¢ S ° < 5; zero otherwise 
b 
. : . f Als 
P(z enerating function of P,, (= & P,,z” SO, 8la\fin\a. ; 
(z) g g a aon ) : J (®)f(x)de =f(9), ifa<o < b) 
+ ati sti + 
Pte) generating function of P} A mean arrival rate 
; -st ransiti ability of going from state i to state j 4 : 
Py one-step transition probability of going : j fe mean service rate 
pY n-step transition probability of going from state 7 to state 9 f Ans Me mean arrival and i 
F rea x : « F ' and service r yom : 
PDK: (cumulative) probability distribution function be: vas ates when there are n in system 
: tA Lig 3 ~~ Moment of a rand ; os 
mae a be ; : i om Vari a ¢ 
p-d.f. probability density function the origin py atin ariable about the mean and about 
as : } v1 = Mean, fig = variance 
m.f, robability mass function 5 ‘ete 2 ) 
‘| é .f I ber i (relative) traffic intensity, also called wtilinats 
Re(:) real part of (a complex number) : as (mean arrival rate is Ge Utilization factor, defined 
: > system)/(mean sarc: 
r.v. random variable : system). It is a bet nig m)/(mean service rate of the 
; er, but is taken ee : 
8 Laplace transform parameter : Erlangs, ken to be measured in 
ae [ o2 3 
SIRO service in random order variance 
alg 
CEO; an : 
a> Fy variance of arrival-or, 4 
of arriy al-group size and service time respectively 


f time | 
. Ss if 
behaves as; J(%) ~ g(x) as «> b if and only if 


aS Sa times of successive arrivals (departures) : 
Uy, U,...  inter-arrival times : lim [fle)/g()) = 1 
Var{-] variance of i cae on 

ee a ae . primes (dashes) denoting differentiation 
W mean (waiting) time in system 

Wa mean (waiting) time in queue or mean queueing time 


Tt Usually 6 is zero or infinite. In the text, 6 is infinite 
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APPENDIX E 
SHORT TABLE OF LAPLACE TRANSFORMS 


We give below some important results in Laplace transforms, most of 


which are used in the text. This table of Laplace transforms can also be 


used as one for inverse-Laplace transforms by reading entries from the second 


column and their inverses from the first column. 


The fu netion f(t) 


(1) f(t) 


(2) Tt” estfis)ds 


2771 c—ia 
(3) af(t)+bg(t), a, 6 constants 


(4) f(t) 
(5) f(t), the n-th derivative of f(t) 


t 
(6) I flaou 


t 

(7) fit) + g(t) = S (w)g(t—w)du 
0 

(8) etf(t) 


(9) f(t), n > 0 


(10) f(t)/t 

(11) 1 

(12) ¢ 

(13) , n> 0 

(14) et 
(15) a(t), Dirac delta function 
(16) I,(at), Re(v) > —1 


(17) tI,(at), Re(v) > —2 


| afis)+b9(8) 

| sf(s)—f(0) ° 

| s” f(s) —s"1f(0) —s"-*f’(0)—... —f (0) 
| F(s)/s 


F(s)g(s) 


fte—a) 
(—1)afimys) = (—£)" Fis) 


f flu)du 
& 

1/s 

1/s? 


T'(n+1)/sm*1 (= n!/s"44, if nis a 
positive integer) 
1/(s—a) 


1 
a(s*—a)-V%(5-+ Vs*@—a?)-”, 
Re(s) > | Re(a)| 
a"(s-+vVv 6*—a*)(s?@a?)-9/2, 
(s-+V's*—a*)-’, Re(s) > | Re(a)| 


yta"(s+-V8—a?), Re(s) > | Re(a) | 
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